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ABSTRACT 


This  experimental  study  was  conducted  at  the  N’aval  Post¬ 
graduate  School  to  investigate  the  heat  transfer  character¬ 
istics  of  flat  vertical  plates  in  a  rectangular  gas - f luidi zed 
bed.  The  primary  objective  was  to  determine  what  effect 
variations  in  the  bed  width-to-height  ratio  had  on  heat  trans¬ 
fer  to  the  vertical  flat  plates  forming  the  container  walls. 

The  experiment  was  conducted  using  a  specially  heated  and  in¬ 
strumented  fluidized  bed  equipped  with  a  movable  side  wall 
which  permitted  modification  of  the  bed  geometry.  As  the 
width  of  the  bed  was  adjusted,  the  settled  bed  height  was  main¬ 
tained  at  a  constant  level  by  the  addition  or  removal  of  bed 
material . 

A  secondary  objective  of  the  study  was  to  determine  the 
effect  of  variations  in  the  fluidization  gas  flow  rate  on  heat 
transfer  to  the  bed  walls.  Flow  rates  ranging  from  fixed  bed 
fluidization  to  pneumatic  conveying  were  studied. 

Pressure  drop  measurements  as  well  as  visual  observations 
were  used  to  determine  minimum  fluidization  flow  rates  for 
each  configuration.  Heat  transfer  coefficients  were  calculated 
for  each  flow  rate  and  bed  geometry  using  temperature  data  ob¬ 
tained  from  a  computer  controlled  thermocouple  network.  In 
addition,  a  study  was  made  of  variations  in  wall  temperature 
with  changes  in  vertical  position. 
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I.  INTRODUCTION 


A  fluidization  bed  in  its  simplest  form  consists  of  a  con¬ 
tainer  filled  with  small  particles  which  are  subjected  to  an 
upward  flow  of  fluid  emerging  from  the  bottom  of  the  container. 
When  the  upward  viscous  force  (drag)  exerted  on  the  particles 
by  the  moving  fluid  balances  or  exceeds  the  downward  force 
of  gravity,  the  particles  become  suspended  and  begin  to  move 
about  the  container  in  a  fluid-like  manner,  hence  the  term 
"fluidization."  This  churning  mixture  of  particles  and  gas 
is  often  referred  to  as  the  "emulsion  phase."  When  the 
velocity  of  the  fluid  is  increased  to  the  point  that  the 
particles  are  carried  out  of  the  container,  a  transition  is 
made  to  a  different  form  of  fluidization  known  variously  as 
pneumatic  conveying,  entrained  flow  fluidization,  or  carryover. 
When  the  velocity  of  the  fluid  is  too  low  to  counteract  the 
gravitational  force  and  cause  fluidization,  the  bed  is  referred 
to  as  being  "packed"  or  "fixed"  and  the  process  is  known  as 
flow  through  a  porous  medium  or  fixed  bed  fluidization. 

Distinctly  different  modes  of  behavior  are  observed  depend¬ 
ing  on  whether  the  fluid  passing  through  the  container  is  a 
liquid  or  a  gas.  Liquid  fluidized  beds  are  often  referred  to 
as  homogeneous  or  particulately  fluidized  beds.  This  type  of 
bed  tends  to  display  a  very  smooth  transition  through  the 
various  stages  of  fluidization  from  fixed  to  fluidized  to 


carryover.  Concurrently  a  very  smooth  transition  in  heat 
transfer  characteristics  and  hydrodynamic  behavior  is  observed. 
On  the  other  hand,  gas  fluidized  beds  tend  to  display  a  more 
erratic  performance.  Most  researchers  attribute  the  difference 
in  behavior  to  the  differing  particle  to  fluid  density  ratios. 
Gas  fluidized  beds  are  known  as  "aggregate"  or  "bubbling 
fluidized  beds."  In  this  type  of  bed,  once  fluidization  has 
begun,  the  gas  usually  forms  into  bubbles  which  displace  the 
bed  particles  as  the  bubbles  rise.  The  appearance  of  a  bubbling 
bed  is  remarkably  similar  to  that  of  a  boiling  liquid,  as  shown 
in  Figure  1.  During  conditions  of  very  active  fluidization 
the  bubbles  sometimes  form  together  or  coalesce  to  create  a 
large  bubble  whose  diameter  approaches  the  width  of  the  bed, 
as  shown  in  Figure  2.  This  condition  is  known  as  "slugging" 
or  "heaving"  and  can  significantly  effect  bed  performance. 
Another  phenomenon  of  gas  fluidized  beds  is  shown  in  Figure  3. 

It  is  a  condition  known  as  "channelling."  When  this  occurs, 
the  gas  travels  through  a  fairly  confined  area  forming  a  gas 
filled  tube  within  the  bed  rather  than  passing  evenly  through¬ 
out  the  cross  section.  This  behavior  is  believed  to  be 
caused  by  poor  distributor  design  and,  like  slugging,  can 
strongly  effect  the  heat  transfer  characteristics  of  the  bed. 

Fluidized  beds  offer  a  wide  variety  of  beneficial  char¬ 
acteristics  such  as  high  rates  of  heat  transfer,  high  thermal 
inertia,  large  solid  to  fluid  contact  area,  and  isothermal 
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bed  temperatures.  In  addition,  once  fluidized,  the  pressure 
drop  across  a  bed  remains  relatively  constant  despite  varia¬ 
tions  in  fluid  flow  rate.  The  pressure  drop  is  a  function 
of  bed  height  and  particle  density  and  is  independent  of 
particle  size.  [Ref.  1]  This  permits  increasing  the  surface 
area  by  using  smaller  particles  without  paying  a  penalty  in 
increased  pressure  loss. 

Fluidized  beds  also  suffer  from  certain  limitations.  The 
size  of  particles  which  can  be  fluidized  is  limited  to  a 
specific  range  extending  from  approximately  0.001  inch  in 
diameter  to  about  0.2S  inch  diameter,  depending  on  the  density 
of  the  particles  and  the  fluid.  [Ref.  2]  In  addition,  some 
materials  cannot  be  fluidized  regardless  of  size  because  of 
their  tendency  to  clump  together  and  form  particles  too  large 
to  fluidize,  or  because  of  their  tendency  to  break  down  into 
particles  which  are  so  small  that  they  become  entrained  in  the 
fluidizing  medium  and  are  carried  from  the  bed.  Another  dis¬ 
advantage  is  that  fluidization  requires  the  expenditure  of 
energy  to  force  the  fluidfzing  medium  through  the  bed  at  the 
required  velocity. 

In  general  the  beneficial  aspects  of  fluidized  beds  far 
exceed  the  limitations  and  as  a  result,  fluidized  bed  prin¬ 
ciples  have  been  utilized  in  many  diverse  applications 
over  the  years.  Early  mining  engineers  used  fluidization 
to  separate  solids  of  different  densities  by  suspending 
them  in  rising  streams  of  water.  [Ref.  1]  Prior  to  World 


IVar  II,  chemical  engineers  discovered  that  fluidized  beds 
made  outstanding  chemical  reactors  because  of  the  extremely 
large  contact  area  between  the  particles  and  the  fluid. 

For  example,  a  container  measuring  one  foot  by  one  foot  by 
one  foot,  filled  with  0.01  inch  diameter  spherical  particles 
yields  approximately  4,000  squares  feet  of  surface  area. 

This  characteristic  of  fluidized  beds  has  been  utilized  to 
great  advantage  in  a  wide  number  of  chemical  processes,  most 
notable  the  catalytic  cracking  of  petroleum.  [Ref.  5] 

Because  the  particles  within  a  fluidized  bed  are  constantly- 
churning  about  and  mixing,  the  temperature  at  one  point  in 
the  bed  is  usually  in  very  close  agreement  with  the  tempera¬ 
ture  at  any  other  point.  Industry  has  capitalized  on  this 
isothermal  property  by  using  fluidized  beds  to  improve  the 
heat  treatment  of  metals  in  such  processes  as  quenching, 
tempering,  and  carburizing.  As  a  combined  result  of  the 
large  surface  area  and  the  thorough  mixing  of  the  particles, 
fluidized  beds  possess  impressive  heat  transfer  capabilities 
Consequently,  the  power  generation  industry  has  begun  utiliz 
ing  this  technology  to  create  a  new  generation  of  heat 
transfer  devices  such  as  steam  generators  which  efficiently 
burn  poor  quality  fuels  while  simultaneously  yielding  very 
low  levels  of  harmful  emissions.  These  devices  are  capable 
of  burning  high  sulfur  coal,  pelletized  garbage,  and  even 
raw  crushed  shale,  just  as  it  comes  from  the  ground.  Anothe 
application  is  seen  in  waste  heat  recovery  units  which  are 


fluidized  by  the  exhaust  gases  from  conventional  boilers  or 
gas  turbines.  The  thermal  energy  which  would  otherwise 
have  been  wasted  is  transferred  to  a  secondary  fluid  cir¬ 
culating  through  tubes  emersed  within  the  bed.  This  re¬ 
covered  energy  is  then  utilized  to  increase  overall  plant 
efficiency.  [Ref.  4] 

The  use  of  fluidized  bed  technology  has  expanded  greatly 
in  recent  years  as  hundreds  of  new  applications  have  been 
discovered,  ranging  from  advanced  ship  propulsion  plants  to 
sophisticated  solar  energy  conversion  and  storage  devices. 

As  the  utilization  of  these  devices  increases  so  does  the 
quest  for  a  better  understanding  of  the  many  complex  phenom¬ 
enon  associated  with  fluidization. 

Fluidized  bed  research  has  been  conducted  throughout  the 
world  and  more  has  been  written  about  fluidization  than  can 
be  read  in  a  lifetime,  however  much  remains  to  be  learned 
in  order  to  unlock  the  full  potential  of  this  fascinating 
process.  This  thesis  represents  a  modest  attempt  to  add  to 
the  knowledge  of  fluidization  in  the  area  of  heat  transfer 
to  vertical  flat  plates  in  gas  fluidized  beds. 
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OBJECTIVES  AND  METHOD  OF  APPROACH 


A.  GENERAL  STATEMENT  OF  THE  PROBLEM 

The  bulk  of  fluidized  bed  heat  transfer  research  conducted 
to  date  has  been  related  to  the  study  of  heat  transfer  to  or 
from  objects  immersed  within  the  bed,  particularly  tubes  and 
tube  bundles.  Comparatively  little  research  has  been  done 
on  the  study  of  heat  transfer  to  the  walls  of  the  containing 
vessel  itself.  In  addition,  a  large  majority  of  the  test 
apparatus  utilized  thus  far  to  collect  exper imenta'l  data 
have  been  cylindrical  in  configuration.  Very  few  studies 
have  been  conducted  on  heat  transfer  within  rectangular 
fluidized  beds;  a  geometry  which  is  being  used  with  increasing 
frequency  in  a  wide  variety  of  applications. 

The  principle  objective  of  this  thesis  was  to  gain  a 
greater  understanding  of  the  heat  transfer  to  the  flat  ver¬ 
tical  plates  which  form  the  walls  of  a  rectangular  gas 
fluidized  bed  and  to  determine  what  effects,  if  any,  varia¬ 
tions  in  bed  geometry  had  on  heat  transfer  performance. 

A  secondary  objective  was  to  study  the  variations  in  heat 
transfer  which  occurred  as  the  fluidization  gas  flow  rate 
was  adjusted  over  a  wide  range. 


B.  METHOD  OF  APPROACH 


This  investigation  was  primarily  experimental  in  nature 
with  limited  utilization  of  analytical  methods  for  comparison 


purposes.  To  acquire  the  necessary  data  for  this  study,  a 
computerized  data  acquisition  system  was  linked  to  a  network 
of  up  to  75  copper-constantan  thermocouples  installed  in  a 
specially  instrumented  and  heated  fluidized  bed.  This  system 
was  used  to  obtain  a  thermal  "snapshot"  of  the  test  apparatus 
which  showed  the  temperature  at  numerous  key  locations  at 
a  specific  point  in  time.  This  data  was  then  combined  with 
measured  values  of  input  power  and  incoming  air  flow  rates 
to  produce  a  calculated  heat  transfer  coefficient.  Experi¬ 
mental  data  was  obtained  for  a  number  of  different  air  flow¬ 
rates  and  bed  conditions  beginning  with  fixed  bed  fluidiza¬ 
tion  and  progressing  through  to  the  onset  of  entrained  flow 
fluidization.  Once  data  had  been  collected  for  the  entire 
range  of  flow. rates  ,  the  bed  geometry  was  altered  by  moving 
an  adjustable  side  wall  within  the  test  apparatus  so  that 
the  width  of  the  rectangular  bed  was  increased.  Additional 
particles  were  then  added  to  the  bed  so  that  the  bed  depth 
remained  at  the  same  level  as  in  the  previous  configuration. 

The  data  collection  process  was  then  repeated  in  its 
entirety  for  the  new  configuration.  A  total  of  four  differ¬ 
ent  bed  geometries  were  studied  in  the  course  of  this  work. 
Over  5,700  temperature  readings  were  recorded  for  eighty  dif¬ 
ferent  test  runs.  All  experimental  data,  including  narrative 
comments  regarding  observed  bed  behavior  were  recorded  on 
magnetic  disks.  The  recorded  data  was  automatically  pro¬ 
cessed  by  computer  and  then  printed  out  in  such  a  manner 
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that  bed  performance  for  a  particular  test  run  was  readily 
apparent.  The  calculated  results  for  the  individual  runs 
were  then  analyzed  to  determine  overall  trends. 

In  addition  to  the  heat  transfer  investigation,  an  ex¬ 
perimental  study  of  bed  pressure  drop  characteristics  was 
conducted  for  each  geometric  configuration  in  order  to 
establish  the  point  of  fluidization  onset  and  to  confirm 
proper  operation  of  the  fluidizing  apparatus. 


III.  EXPERIMENTAL  APPARATUS  AND  PROCEDURES 


A.  DESCRIPTION  OF  FLUIDIZATION  APPARATUS 

The  fluidization  apparatus  consisted  of  a  rectangular 
plexiglas  box  filled  with  0.012  inch  diameter  silica  glass 
beads  as  shown  in  Figure  4.  The  box  was  equipped  with  a 
movable  side  wall  which  permitted  the  inside  dimensions  of 
the  bed  to  be  varied  from  6.25  inches  wide  by  6  inches  deep 
by  18  inches  high  to  24  inches  wide  by  6  inches  deep  by  18 
inches  high.  The  movable  wall  was  equipped  with  felt  gasket 
which  prevented  the  escape  of  beads  through  the  gap  between 
the  movable  side  wall  and  the  fixed  front  and  rear  walls. 

0.5  inch  thick  clear  plexiglas  (grade  GM)  was  used  through¬ 
out  the  construction  and  all  components  were  assembled  using 
plastic  laminating  cement  with  wood  screw  reinforcements. 

The  floor  of  the  apparatus  also  served  as  the  distributor 
plate  whose  function  was  to  evenly  distribute  the  flow  of 
incoming  air  over  the  entire  cross  section  of  the  bed.  The 
distributor  plate  was  constructed  from  a  plexiglas  sheet 
with  0.125  inch  diameter  holes  drilled  0.25  inches  on  center 
Secured  over  this  plate  was  a  140  mesh  stainless  steel  wire 
screen  to  prevent  the  beads  from  falling  down  into  the  holes 
and  to  provide  an  even  finer  distribution  of  air. 

-  Located  beneath  the  distributor  plate  was  an  air  chamber 
or  plenum  which  served  to  ensure  that  the  air  pressure  on 


the  underside  of  the  distributor  plate  was  the  same  at  any 
point.  This  was  accomplished  by  first  directing  the  incoming 
air  through  a  2.5  inch  inside  diameter  plexiglas  distribution 
tube  centered  lengthwise  within  the  24  inch  by  6  inch  by  6 
inch  plenum.  This  tube  had  0.5  inch  diameter  holes  drilled 
on  1  inch  centers  on  its  underside  and  served  to  evenly  direct 
the  incoming  air  downward  over  the  entire  length  of  the  plenum 
In  order  to  prevent  the  escape  of  air  when  the  bed  was  con¬ 
figured  to  a  width  less  than  the  24  inch  width  of  the  distri¬ 
butor  plate,  the  movable  side  wall  was  equipped  with  a 
horizontal  baffle  plate  which  laid  flat  over  the  top  of  the 
unused  portion  of  the  distributor  plate.  Once  the  desired 
position  of  the  movable  wall  was  set,  the  baffle  plate  was 
securely  clamped  to  the  distributor  plate  and  the  edges  were 
sealed  with  silicon  caulking.  After  passing  through  the 
distributor  plate  and  finally  the  bed  itself,  the  fluidizing 
air  escaped  from  the  apparatus  via  one  or  more  of  three 
plexiglas  exhaust  tubes  positioned  on  the  top  of  the  apparatus 
These  tubes  were  quipped  with  wire  mesh  covers  to  prevent  the 
loss  of  glass  beads  due  to  carryover. 

In  order  to  measure  the  air  pressure  at  various  depths 
the  apparatus  was  equipped  with  static  pressure  taps  posi¬ 
tioned  along  the  rear  wall  with  2  inch  vertical  spacing 
between  centers.  Additionally  taps  were  installed  just 
below  the  distributor  plate  and  immediately  above  it.  These 
taps  were  connected  by  flexible  plastic  tubing  to  a  manifold 


which  permitted  selective  connection  of  the  taps  to  a 
Meriam  Model  33KA35  .30  inch  water  manometer  as  shown  in 
Figure  5. 

The  fluidizing  apparatus  was  supplied  with  air  from  a 
Spencer  axial  flow  turbo  compressor  rated  at  500  SCFM. 

The  fluidizing  air  was  not  filtered  or  dehydrated  prior 
to  entry.  The  volumetric  flow  rate  of  air  entering  the 
apparatus  was  measured  by  a  calibrated  Fisher  and  Porter 
flowmeter  rated  at  34  CFM. 

B.  DESCRIPTION  OF  HEAT  TRANSFER  APPARATUS 

t 

The  test  apparatus  was  equipped  with  two  instrumented 
flat  plate  heater  assemblies  each  measuring  5  inches  wide 
by  10  inches  high.  These  assemblies  were  mounted  flush  with 
the  inside  surface  of  both  the  movable  side  wall  and  the 
fixed  side  wall  opposite  it.  Each  heater  assembly  was 
powered  by  a  Watlow  strip  heater  rated  at  250  watts.  The 
heater  consisted  of  a  grid  of  nickel  alloy  resistance  wire 
embedded  in  a  thin  sheet  of  rubberized  supporting  material 
and  backed  by  a  0.5  inch  thick  sheet  of  silicone  sponge 
rubber  insulation.  This  insulation  was  in  turn  bonded  to 
a  0.5  inch  thick  clear  plexiglas  mounting  block.  The  face 
of  the  assembly  consisted  of  a  0.25  inch  thick  copper  plate 
mounted  flush  with  the  heating  surface  of  the  strip  heater. 
In  order  to  avoid  the  use  of  screws  or  other  conventional 
attachment  methods  which  may  have  altered  or  impeded  the 


heat  flow  from  the  strip  heater  to  the  copper  plate,  the 
two  were  held  together  using  a  strongback  arrangement  as 
shown  in  Figures  6  and  7.  This  mounting  method  ensured 
that  the  strip  heater  was  held  firmly  against  the  copper 
plate,  thus  helping  to  minimize  thermal  contact  resistance 
between  the  two.  To  further  reduce  this  resistance,  a  ther¬ 
mally  conductive  paste  was  applied  to  the  heating  surface 
of  the  strip  heater  prior  to  positioning  the  copper  plate 
over  it.  Power  was  supplied  to  the  left  hand  heater  by  a 
Lambda  regulated  DC  power  supply  model  LK34SA,  while  the 
right  hand  heater  was  supplied  by  a  Hewlett  Packard  DC 
power  supply  Model  6296A.  Both  power  supplies  were  rated 
at  60  Volts  and  2  Amps  and  were  each  equipped  with  an  in¬ 
dividual  voltmeter  and  ammeter. 

All  thermocouples  used  in  this  apparatus  were  10  gauge, 

ANSI  Type  T,  copper-constantan ,  exposed  junction  type.  Each 

heater  assembly  was  equipped  with  twenty  thermocouples 

arranged  in  five  rows  of  four  thermocouples  each  as  shown 

in  Figures  8  and  9.  These  thermocouple's  were  mounted  to  the 

copper  plate  in  such  a  manner  so  as  to  be  least  disruptive 

of  the  heat  distribution  patterns  within  the  plate.  The 

thermocouple  junctions  were  soldered  into  small  shallow  holes 

drilled  into  the  back  of  the  plates  so  that  the  top  of  the 

junction  was  flush  with  the  plate  surface.  The  thermocouple 

leads  were  then  routed  straight  back  through  holes  drilled 
«  * 

into  the-  strip  heater  and  its  mounting  block.  Figure  10 


Figure  6.  Heater  Assembly  Mounting  Diagram;  Side  View 
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Figure  7.  Heater  Assembly  Mounting  Diagram;  Rear  View 
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NOTE:  Arrangement  is  as  seen  from  within  the 
bed  (front  view). 


Figure  8.  Right  Hand  Heater  Thermocouple  Arrangement 
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shows  the  thermocouple  mounting  procedure.  X-ray  photographs 
were  used  to  determine  the  exact  location  of  the  resistance 
wires  embedded  within  the  heaters  so  as  to  avoid  striking 
them  while  drilling  the  holes.  This  thermocouple  arrangement 
made  it  possible  to  read  the  plate  inside  surface  temperature 
immediately  adjacent  to  the  surface  of  the  strip  heater  while 
causing  minimum  disruption  of  the  heat  flux  patterns  within 
the  plate  itself.  A  determination  of  the  plate  outside 
surface  temperature  as  seen  by  the  bed  was  then  accomplished 
by  means  of  simple  calculations  as  shown  in  section  1 1 1 H . 

The  temperatures  within  the  bed  itself  were  monitored  by 
15  retractable,  24  inch,  stainless  steel  sheathed,  exposed 
junction,  Omega  Engineering  Type  T  thermocouple  probes.  The 
position  of  these  probes  was  held  fixed  at  a  bed  depth  of 
7  inches  above  the  distributor  plate  for  all  data  runs.  The 
number  of  bed  probes  actually  utilized  depended  on  the  geo¬ 
metric  configuration  being  studied,  as  shown  in  Figure  11. 

Incoming  air  temperature  was  monitored  by  3  retractable, 

4  inch,  stainless  steel,  sheathed  thermocouple  probes  posi¬ 
tioned  within  the  air  plenum  as  shown  in  Figure  12.  Out¬ 
going  air  temperature  was  monitored  by  3  similar  probes 
mounted  one  each  in  the  center  of  the  3  air  exhaust  tubes, 
as  shown  in  Figure  13.  Because  of  the  positioning  of  these 
tubes  with  respect  to  the  movable  side  wall  only  one  thermo¬ 
couple  was  used  to  read  the  air  outlet  temperature  for  all 
data  runs. 
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NOTE:  Arrangement  is  as  viewed  from  top  front 


Figure  12.  Air  Inlet  Thermocouple  Probe  Arrangement 


NOTE:  Arrangement  is  as  viewed  from  top  front. 


Figure  13.  Air  Outlet  Thermocouple  Probe  Arrangement 


Figure  14.  Right  Hand  Heater  Mounting  Block  Thermocouple 
Arrangment 


Figure  IS.  Left  Hand  Heater  Mounting  Block  Thermocouple 
Arrangment 
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Three  insulated  thermocouples  were  mounted  to  the  outside 
surface  of  each  of  the  plexiglas  heater  mounting  blocks,  as 
shown  in  Figures  14  and  15,  in  order  to  determine  the  mounting 
block  outer  surface  temperatures.  These  temperatures  were 
used  to  estimate  the  amount  of  heat  lost  through  the  back  of 
the  heater  assemblies  as  explained  in  section  1 1 1 E .  Thermo¬ 
couples  were  also  mounted  on  the  inside  and  outside  surfaces 
of  the  plexiglas  sidewalls,  as  shown  in  Figures  16  and  17, 
in  order  to  estimate  the  losses  to  the  atmosphere  and  to  the 
bed  itself  from  the  thin  plexiglas  wall  strips  adjacent  to 
the  edges  of  the  heater  assemblies.  A  final  set  of  thermo¬ 
couples  was  mounted  to  the  inside  and  outside  surfaces  of 
the  container  rear  wall,  as  shown  in  Figure  18,  in  order  to 
obtain  an  estimate  of  losses  through  the  rear  wall.  Because 
of  symmetry  it  was  assumed  that  the  losses  through  the  front 
wall  were  identical  to  the  rear  wall  losses. 

All  thermocouples  used  in  the  test  apparatus  were  sampled 
by  an  HP  3497A  data  acquisition  system  controlled  by  an  HP  85 
desk  top  computer.  This  system  was  capable  of  sampling  and 
storing  the  output  of  seventy  five  thermocouples  within  seven 
seconds.  Because  of  the  steady  state  conditions  in  existance 
during  each  of  data  runs,  the  seven  second  sampling  time  was 
considered  brief  enough  to  treat  the  thermocouple  readings 
as  through  they  had  all  been  taken  simultaneously. 

Once  all  of  the  thermocouples  had  been  sampled,  the  HP  85 
retrieved  the  stored  thermocouple  voltages  and  converted  them 
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to  temperature  values  using  equations  previously  developed. 
These  temperature  values,  along  with  other  stored  parameters, 
were  then  recorded  on  magnetic  disks  using  a  Hewlett  Packard 
Model  HP  82901M  dual  flexible  disk  drive.  This  recorded 
data  was  later  processed  and  permanently  displayed  using  an 
HP  82905B  printer  and  an  HP  7225A  plotter. 

Computer  programs  used  to  control  data  acquisition  and 
data  processing  are  shown  in  Appendixes  J  and  K,  respectively. 
The  arrangement  of  the  heat  transfer  data  collection  and  pro¬ 
cessing  system  is  shown  in  Figure  19  while  Figures  20  and  21 
show  photographic  views  of  the  entire  experimental  system. 

C.  EXPERIMENTAL  APPARATUS  CALIBRATION 

The  inlet  air  flowmeter  was  calibrated  using  an  ASME 
Herschel- type  venturi  with  an  inlet  diameter  of  4.26  inches 
and  a  throat  diameter  of  2.13  inches.  Experimental  data 
(shown  in  Appendix  B)  was  collected  for  nine  different  data 
points  and  from  this  data,  volumetric  air  flow  rates  were 
calculated  using  equations  shown  in  Appendix  H.  A  plot  of 
air  flowrate  vs  flowmeter  reading  is  shown  in  Figure  22.  A 
least  squares  curve  fit  was  applied  to  the  data  points  to 
obtain  a  linear  equation  which  converts  flowmeter  readings 
to  air  flow  rates  in  CFM.  This  equation  was  then  incorporated 
into  the  data  reduction  computer  program  shown  in  Appendix  K. 
The  rotameter  calibration  data  shows  that  a  1001  flowmeter 
reading  equates  to  34.0  CFM. 


Figure  20.  Photograph  of  experimental  Apparatus 


Lgure  21.  Photograph  of  Experimental  Apparatus 


Thermocouples  number  1  through  67  were  calibrated  using 
a  Rosemount  Commutating  Bridge  Model  920A  and  a  Rosemount 
Calibration  Bath  Model  913A  using  liquid  nitrogen  and  water. 
The  HP  3497A  and  the  HP  85  were  included  in  the  system  during 
calibration,  as  shown  in  Figure  23.  Data  was  collected  for 
twenty  different  bath  temperatures  and  a  least  squares  curve 
fit  was  applied  to  the  data  to  obtain  a  fourth  order  poly¬ 
nomial  function  for  each  of  the  thermocouples.  The  constants 
from  these  functions  were  then  incorporated  into  the  data 
acquisition  computer  program  shown  in  Appendix  J  to  allow 
automatic  conversion  of  thermocouple  output  voltage  to  tem¬ 
perature,  in  degrees  Fahrenheit.  A  sample  of  the  calibration 
data  in  shown  in  Appendix  D. 

Thermocouples  number  68  through  79  were  not  calibrated 
but  were  instead  assigned  constants  from  the  standard  power 
expansion  for  Type  T  thermocouples.  [Ref.  5]  They  were  then 
compared  in  a  water  bath  of  known  temperature  to  confirm 
proper  operation. 

The  water  manometer  used  to  determine  static  pressure 
drop  within  the  bed  was  regularly  rechecked  for  zero  readings 
and  was  readjusted  as  necessary. 

D.  EXPERIMENTAL  PROCEDURE 

Preparation  for  the  taking  of  heat  .transfer  data  was 
accomplished  by  adjusting  the  movable  side  wall  to  the  proper 
position  for  the  geometry  to  be  studied.  The  baffle  plate 
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was  then  clamped  into  position  and  the  joints  were  sealed 
with  silicone  rubber  caulking.  This  caulking  was  allowed 
to  harden  for  a  minimum  of  twenty  four  hours  to  ensure  a 
leak  tight  seal.  Carefully  weighed  quantities  of  glass 
beads  were  then  added  to  the  bed  through  the  right  hand 
air  exhaust  tube  until  the  depth  of  the  bed  reached  12.5 
inches  above  the  distributor  plate.  The  air  supply  turbo 
compressor  was  then  lined  up  and  energized.  Following  this, 
the  air  inlet  control  valve  was  opened  and  adjusted  until 
the  desired  air  flow  rate  registered  on  the  flowmeter.  The 
heater  power  supplies  were  then  energized  and  adjusted  to 
the  desired  level. 

The  power  supply  controls  were  set  so  that  the  voltage 
going  to  each  of  the  heaters  was  the  same.  Current  into 
the  heaters  was  determined  by  the  resistance  of  the  nickel¬ 
heating  element  which  was  in  turn  determined  by  the  temper¬ 
ature  of  the  heaters.  Once  the  voltage  was  equilized,  the 
bed  was  allowed  to  come  to  steady  state  conditions.  Prior 
to  the  taking  of  data,  the  system  was  checked  for  thermal 
equilibrium  by  monitoring  the  temperature  at  specific  points 
using  a  special  computer  program  which  sampled  thermocouples 
at  five-second  intervals.  Once  it  was  confirmed  that  the 
heater  assemblies  and  the  bed  itself  were  at  a  steady  tem¬ 
perature,  the  data  acquisition  program  was  initiated.  At 
the  beginning  of  the  program  the  computer  requested  that 
the  operator  provide  certain  information  regarding  the  data 
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run  to  be  recorded.  The  first  item  of  information  requested 
was  the  run  number.  This  number  was  structured  so  as  to  be 
as  informative  as  possible.  The  first  two  digits  represented 
the  bed  width  in  inches.  These  two  digits  were  followed  by 
a  letter  indicating  which  of  the  heaters  was  energized;  "L" 
designated  left,  "R"  stood  for  right,  and  "B"  designated  both. 
This  letter  was  followed  by  two  digits  indicating  the  air 
flowmeter  reading.  Finally,  a  single  letter  was  used  for 
special  purpose  designations  such  as  a  repeated  run  without 
a  change  in  configuration  or  flow  rate.  The  next  item  of 
information  requested  by  the  program  was  the  date  and  time. 
This  was  followed  by  requests  for  bed  width  and  height,  air 
flow  rate,  ambient  temperature  (read  on  a  mercury  in  glass 
thermometer  positioned  near  the  test  apparatus),  voltage 
and  current  readings  from  the  two  power  supplies  to  each  of 
the  heaters,  the  number  of  bed  thermocouple  probes  in  use, 
and  the  amount  of  bed  expansion,  if  any.  Finally,  general 
comments  concerning  observed  bed  behavior  were  entered. 

As  the  push  of  a  button,  all  of  the  thermocouples  were 
then  sampled  and  the  corresponding  temperatures  were  then 
recorded  on  disc  along  with  the  other  data  entered  from  the 
keyboard.  Once  the  recording  process  was  completed,  the 
system  was  ready  to  begin  another  data  run.  If  the  next 
run  to  be  studied  was  for  the  same  bed  configuration,  all 
that  ha4  to  be  done  was  adjust  the  air  flow  rate  to  the  new 
reading  and  wait  for  the  temperatures  to  stabilize.  In 
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some  cases,  the  power  supplied  to  the  heaters  was  also  ad¬ 
justed  in  order  to  maintain  as  large  a  temperature  differ¬ 
ential  as  possible  between  the  incoming  air  and  the  heaters. 

In  order  to  avoid  excessive  softening  of  the  plexiglas 
adjacent  to  the  heaters,  the  maximum  heater  temperature  was 
limited  to  approximately  one  hundred  eighty  degrees 
Fahrenheit . 

After  data  was  collected  over  the  entire  range  of  flow 
rates  for  a  given  geometry,  the  caulking  was  removed  and  the 
movable  wall  was  repositioned  for  the  next  configuration. 

The  entire  data  collection  process  was  then  repeated.  In 
all,  four  different  geometric  configurations  were  studied. 

Bed  widths  of  6.25  inches,  8  inches,  10  inches,  and  12  inches 
were  used.  In  addition,  data  runs  were  taken  in  the  6.25 
inch  configuration  with  the  bed  empty.  Data  was  also  col¬ 
lected  with  the  underside  of  the  distributor  to  the  left  of 
the  movable  wall  masked  off  as  a  check  to  see  if  the  baffle 
plate  was  performing  effectively.  A  total  of  eighty  data  runs 
were  completed  during  this  study.  Experimental  heat  transfer 
data  as  well  as  calculated  results  are  shown  in  Appendix  C. 

Pressure  drop  data  was  collected  for  each  bed  configura¬ 
tion  by  aligning  the  valve  manifold  (shown  in  Figure  S)  so 
that  the  static  pressure  tap  immediately  above  the  distri¬ 
butor  plate  was  connected  to  the  30  inch  water  manometer. 

The  fluidizing  air  flow  rate  was  then  adjusted  to  the  desired 
level  and  the  pressure  reading  at  the  base  of  the  bed  was  read 


on  the  manometer  and  recorded.  The  manifold  was  then  re¬ 
aligned  so  that  the  static  pressure  tap  near  the  top  of 
the  bed  housing  was  connected  to  the  manometer.  This  read¬ 
ing  was  then  read  and  recorded.  By  substracting  this  read¬ 
ing  from  that  obtained  at  the  base  of  the  bed,  the  bed 
pressure  drop  was  determined.  Pressure  drop  data  was  ob¬ 
tained  in  this  fashion  over  the  entire  flow  range  for  each 
of  the  bed  configurations.  During  the  pressure  drop  data 
runs,  the  power  to  each  of  the  heaters  was  secured. 

E.'  EXPERIMENTAL  DATA  PROCESSING 

All  raw  data  was  processed  using  an  HP  85  desk  top 
computer.  The  final  results  consisted  of  six  pages  of 
printed  data  for  each  run.  The  first  page  was  a  listing  of 
experimental  conditions  and  comments  recorded  during  the 
data  acquisition  phase  of  the  study. 

The  second  page  was  a  listing  of  all  thermocouple  out¬ 
puts  and  their  corresponding  temperatures. 

The  third  page  showed  a  temperature  profile  of  the  right 
hand  heater  surface  as  seen  from  within  the  bed.  In  addition, 
an  average  value  of  the  thermocouple  readings  in  each  hori¬ 
zontal  row  was  calculated  and  displayed.  Thermocouple  number 
four  was  found  to  be  malfunctioning therefore  it’s  output 
was  disregarded.  In  a  separate  display,  the  average  row 
temperatures  were  plotted  against  bed  depth  to  show  how  the 
heater  temperatures  varied  from  top  to  bottom. 
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The  fourth  page  of  the  results  printout  showed  the  left 
hand  heater  temperature  profile  in  a  format  identical  to 
that  of  page  three. 

The  fifth  page  showed  the  temperature  profile  of  the  bed 
itself,  as  measured  by  the  thermocouple  probes  immersed 
within  it.  The  profile  is  arranged  as  seen  from  above  the 
bed . 

The  sixth  and  final  page  contained  the  calculated  results. 
The  first  value  displayed  on  this  page  was  the  average  heater 
temperature  (T^) .  This  value  was  a  simple  arithmetic  average 
of  the  five  horizontal  row  temperatures  previously  calculated. 

Th  =  {R(l)  +  R(2)  *  R(3)  i-  R  (4 )  +  R(5)}/5  (3.1) 

The  average  bed  temperature  (T^)  was  similarly  calculated 
using  the  temperatures  displayed  on  page  five  of  the  printout. 

Tb  =  {T  (50)  +  T  (51 )  +  T ( 5  2 )  +  T(53)  +  T (54 )  +  T(56) 

+  .  +  T  (X)  }/  (X-  50)  (3.2) 

X  *  Highest  number  bed  probe 

The  electrical  energy  into  the  heater  (q^g)  was  calculated 
from  the  voltage  and  current  measurements  read  off.  The  heater 
power  supplies  as  shown  below: 


lie 


I  x  E  x  0.0569  ( 


Btu 


Watt -Min 


) 


(3.3) 


An  estimate  of  the  heat  less  from  the  heater  to  the  atmos¬ 
phere,  through  the  heater  mounting  block  (q^)  was  obtained 
as  follows:  First  an  average  temperature  (T  ^)  was  calculated 
for  the  three  thermocouples  attached  to  the  outer  surface  of 
the  heater  mounting  block. 


Tmb  =  {T(41)  +  T(42)  +  T (4 3) }/ 3 


(3.4) 


Next  it  was  assumed  that  the  thermal  contact  resistance 
between  the  copper  plate  and  the  strip  heater  was  negligible, 
therefore  the  temperature  of  the  strip  heater  surface  was 
equal  to  T^.  The  contact  resistance  between  the  strip  heater 
and  the  insulation  was  also  assumed  to  be  negligible. 

The  heater  insulation  and  mounting  block  was  then  con¬ 
sidered  to  be  a  composite  wall  as  shown  in  Figure  24.  Fourier’s 
Law  of  Conduction  was  then  utilized  to  calculate  the  heat  flux 
through  the  mounting  block. 


Th  '  Tmb 
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(3.5) 
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Figure  24.  Heater  Assembly  Temperature  Profile 
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0.5  inches 


Then  an  estimate  was  made  of  the  heat  lost  to  the  atmos¬ 
phere  from  the  hot  plexiglas  side  wall  sections  adjacent  to 
the  side  edges  of  the  heater  (q ^ ) -  This  estimate  was  based 
on  the  side  wall  outer  surface  temperature  and  used  Newtons 
Law  of  Cooling  and  the  convection  heat  transfer  coefficient 
previously  calculated  for  the  heater  mounting  block.  Since 
the  side  wall  sections  and  the  mounting  block  were  in  such 
close  proximity  it  was  assumed  that  the  coefficients  would 
be  in  close  agreement  with  one  another. 

<■„  ■  hmba  *  As  *  (T(72)  '  V  <3 

1  Ft 2 

Ag  =  0.5  in.  x  10  in.  x  ^44  ~in2  x  2 

(one  strip  on  each  side  of  heater) 
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A  similar  loss  estimate  was  then  calculated  for  the  heat 
flux  to  the  atmosphere  from  the  side  wall  sections  adjacent  to 
the  top  edge  of  the  heater  (q  ) . 

L  cL 

"ta  *  hmba  x  At  *  <T<68>  '  V  (3- 

1  pt  ^ 

Af  *  0.5  in.  x  6  in.  x  — : - T 

144  in" 

Yet  another  estimate  was  then  calculated  for  the  loss  to 
the  atmosphere  from  the  side  wall  section  adjacent  to  the 
bottom  edge  of  the  heater  (q^a)  • 


%a  =  hmba  x  Ab  x  {T(71)  '  T a> 


(3.9) 


1  Ft 

A.  =  0.5  in.  x  6  in.  x  - =- 

b  144  in2 


All  atmospheric  losses  were  then  totaled, 


<a  =  ^mb  +  ^sa  +  ^ta  +  ^ba  .  [3-10 

The  difference  between  the  electrical  power  input,  qei, 
and  losses  from  the  heater  to  the  atmosphere,  q  ,  is  the  heat 

a 

flux  transferred  to  the  bed.  In  order  to  identify  that  portion 
transferred  directly  from  the  heated  copper  plate  to  the  bed, 
an  estimate  of  the  heat  flux  from  the  heated  plastic  strip 
surrounding  the  copper  plate  to  the  bed  must  be  made.  This 
is  done  by  assuming  that  same  heat  transfer  coefficient  applies 
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to  both  the  copper  plate  and  to  the  plastic  strip.  The  method 
selected  was  an  iterative;  wherein  it  was  first  assumed  that 
the  heat  flux  passed  through  the  copper  plate,  from  which  the 
heat  transfer  coefficient  was  calculated.  Using  this  coeffi¬ 
cient,  the  heat  flux  from  the  heated  plastic  strip  was  cal¬ 
culated  which  in  turn  yielded  a  new  heat  flux  from  the  copper 
plate.  The  process  was  repeated  until  the  heat  flux  values 
did  not  change  significantly.  Thus  following  this  procedure 
the  heat  flux  into  the  bed  from  the  heater  q  ^  was  assumed 
to  be  equal  to  the  incoming  electrical  power  qecJ  minus  the 
losses  from  the  heater  to  the  atmosphere  qa> 


^cb  ^ei  '  ^a 


(3.1 


Once  the  heat  flux  through  the  copper  plate  qcb  was  known, 
the  plate  outside  surface  temperature  T^  seen  by  the  bed 
could  be  calculated. 


fqcb  x  Lc, 

lJc — T~K  ‘ 


Lc  -  0.25  in.  x 

kc  =  232  -  (0.032  x  (Th  -  70)}  Btu/Hr-Ft-°F 


(3.1 


A 


1  Ft 

5  in.  x  10  in.  x  (tthi — rx2) 


Once  the  outside  plate  surface  temperature  was  known,  the 


convection  heat  transfer  coefficient  from  the  plate  to  the 
bed  (hc^)  could  be  calculated. 
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(3.1 


Using  the  plate  to  bed  heat  transfer  coefficient  (h  b) 
just  determined,  an  estimate  was  made  of  the  heat  flux  which 
reached  the  bed  from  the  hot  plexiglas  strips  adjacent  to  the 
heater  sides  (qsb)  • 


lsb 


-  hcb  x  As  x  (T (73) 


V 


(3.1 


Similar  estimates  were  then  prepared  for  the  flux  entering 
the  bed  from  the  hot  sidewall  section  adjacent  to  the  top  edge 
of  the  heater  (qtb)  as  well  as  the  section  adjacent  to  the 
bottom  edge  of  the  heater  (q^^ ) . 

itb  ■  hcb  *  \  x  tTl69)  ■  V  f3-1 

%b  •  hcb  *  Ab  x  tT(70)  •  Tb}  (3a 

The  heat  flux  into  the  bed  from  all  sources  other  than  the 


heater  plate  itself  were  then  added  to  obtain  a  total  value 

(<W- 


qob  =  qsb  +  qtb  +  qbb 


(3.17) 


The  original  copper  plate  heat  flux  value  was  then  cor¬ 
rected  to  allow  for  that  amount  (q^)  which  did  not  pass 
through  the  plates. 

<>cb  *  %i  •  <ta  •  ‘lob  t3'ls) 

It  was  then  necessary  to  recalculate  all  of  the  parameters 
dependent  on  this  value.  This  process  was  repeated  in  an  iter¬ 
ative  fashion  until  the  calculated  flux  through  the  plate  did 
not  change  signf iciantly .  Once  the  iterative  process  had  been 
completed,  the  final  values  were  printed  out  and  the  entire 
calculation  process  was  repeated  for  the  left  hand  heater. 

The  left  hand  heater  calculations  were  similar  to  those  used 
for  the  right  hand  with  the  exception  of  variable  names  and 
thermocouple  numbers.  Because  of  a  limitation  on  the  total 
number  of  thermocouples  which  could  be  controlled  by  the 
HP  3497A,  there  were  no  thermocouples  available  to  measure 
the  inside  and  outside  surface  temperature  of  the  plexiglas 
sidewall  section  adjacent  to  the  bottom  edge  of  the  left 
hand  heater.  This  prevented  the  calculation  of  a  loss 
estimate  based  on  a  direct  measurement.  To  compensate  for 
this,  it  was  assumed  that  the  losses  from  the  sidewall 
sections  adjacent  to  the  top  and  bottom  edges  of  the  heater 
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were  identical  to  one  another.  Once  the  left  hand  heater 
computations  were  completed  and  the  results  displayed,  cal¬ 
culations  were  begun  on  an  energy  balance.  The  objective 
was  to  determine  the  total  amount  of  energy  entering  the  bed 
and  then  compare  that  with  the  amount  calculated  to  be  leaving 
in  order  to  validate  the  data  and  the  calculations  used.  The 
first  parameter  to  be  determined  was  the  density  of  the  air 

within  the  bed  o„ . 

a 


1.325  x  Pa 

(T  +459 .69)  (5.19 

cl 

P  =  29.92  in.  Hg 
a  d 

Ta  =  T(S5) 


Since  the  air  entering  a  fluidized  bed  is  known  to  assume 
the  bed  temperature  within  a  very  short  distance  of  entry, 
the  air  outlet  temperature  was  used  as  the  bed  air  temperature 
Using  this  value  of  density,  the  amount  of  energy  carried 
out  of  the  bed  by  the  fluidizing  air  stream  (qao)  was 
calculated . 


q  =  m  C  (T  -  T  ■  )  (3.20 

Hao  p  v  ao  ai'  v 

qao  =  ^  x  pab  x  60  Min/Hr  x  °*241  x  (Ta0  ‘  Tai)  (3.21 


T  •  =  (T (47 )  +  T (4 8 )  +  T(49)}/3 


The  heat  loss  through  the  front  and  rear  walls  (qfra)  was 
then  estimated.  Only  the  rear  wall  was  instrumented,  however, 
because  of  symmetry,  it  was  assumed  that  the  loss  through 
the  front  wall  would  be  identical. 


°*f  ra 


kfr  x  Afr  x  {T  (7  5 ) 


L 


fr 


T (74 )  } 


(3.22) 


kf  =  0.112  Btu/Hr- Ft  - °F 

1  Ft 

=  12  in.  x  W  x  ( j-g-A  -  )  x  2  (to  account  for 
both  front  and  rear  walls) 

W  =  Bed  Width 

L-  =  0.5  in. 


The  total  heat  flux  leaving  the  bed  qQ  was  then  calculated. 
It  was  assumed  that  the  only  flux  paths  out  of  the  bed  were 
via  the  front  and  rear  walls  and  via  the  outgoing  air  stream. 
The  side  walls  were  considered  to  be  at  a  sufficiently  elevated 
temperature  as  to  block  the  flow  of  heat  outward  from  the  bed. 
The  floor  of  the  bed  was  considered  to  be  free  of  losses  since 
any  heat  passing  through  it  would  be  returned  to  the  bed  by 
the  fluidizing  air  stream  as  it  passed  through  the  distributor. 


%  =  %a  +  ^f  ra 

The  total  energy  into  the  bed  from  all  right  hand  sources 
was  then  calculated  (q^)  •  This  included  both  the  flux 
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Following  the  completion  of  this  calculation  the  energy  balance 
section  of  page  six  was  printed  out  and  the  program  was 
terminated . 
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IV.  PRESENTATION  AND  DISCUSSION  OF  RESULTS 


A.  FLOW  PATTERNS 

The  particle  flow  patterns  observed  in  the  course  of  this 
study  were  in  general  agreement  with  those  described  in  the 
literature  for  a  bubbling  bed.  The  particles  moved  upward 
near  the  center  line  of  the  bed  as  did  most  of  the  bubbles. 

Once  the  particles  reached  the  surface  of  the  bed  they  were 
thrown  outward  from  the  center  and  then  proceeded  to  travel 
downward  along  the  bed  walls  as  shown  in  Figure  25.  The 
area  of  greatest  activity  was  roughly  cylindrical  in  shape 
with  the  corner  areas  exhibiting  considerably  less  particle 
motion.  Observed  particle  stagnation  regions  are  shown  in 
Figure  26.  No  indications  of  channeling  were  observed  visually 
at  any  time.  At  high  fluidization  flow  rates  in  the  narrow 
bed  configuration,  slugging  was  occasionally  observed.  This 
was  characterized  by  the  appearance  of  large  bubbles  and  the 
sudden  lifting  of  the  bed  surface  to  a  level  two  to  three 
inches  above  its  previous  position.  Once  the  bubble  reached 
the  top  of  the  bed,  the  surface  would  suddenly  collapse  and 
resume  its  former  level  until  the  next  bubble  appeared. 

The  particles  within  the  bed  were  observed  to  behave  in 
a  manner  much  like  a  viscous  fluid.  The  layer  of  particles 
directly  adjacent  to  the  walls  showed  little,  if  any,  move¬ 
ment.  Because  of  their  transparent  composition  it  was  possible 
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Figure  25.  Observed  Particle  Circulation  Flow  Paths 


Figure  26.  Observed  Particle  Stagnation  Regions 


to  observe  that  the  next  layer  of  particles  was  moving  at  a 
slightly  higher  velocity  and  that  the  third  layer  inboard 
was  moving  even  faster.  This  movement  was  not  steady  but 
tended  to  occur  in  surges.  It  is  possible  that  a  static 
charge  contributed  to  the  "No  slip"  condition  observed  at 
the  boundary.  After  the  fluidizing  air  flow  was  secured  upon 
completion  of  a  test  run  it  was  noted  that  numerous  particles 
were  clinging  to  the  container  walls  in  the  upper  portion  of 
the  bed.  In  an  attempt  to  eliminate  this  condition,  an  anti¬ 
static  film  was  sprayed  on  the  upper  portions  of  the  bed  walls. 
For  a  short  period  of  time,  the  walls  remained  free  of  cling¬ 
ing  particles,  however,  the  effect  did  not  last  for  more  than 
two  or  three  data  runs.  Following  this,  no  further  attempts 
were  made  to  alter  this  phenomenon.  Research  by  Miller  and 
Logwinuk  [Ref.  2]  suggested  that  the  presence  of  such  a  static 
charge  could  contribute  to  erratic  heat  transfer  results,  how¬ 
ever  for  the  purpose  of  this  study,  any  affect  was  considered 
to  be  negligible. 

Closer  examination  of  the  bubble  flow  path  showed  that 
many  bubbles  tended  to  originate  from  the  left  side  of  the 
bed  near  the  bottom  of  the  movable  side  wall.  As  a  result, 
there  was  a  shift  in  overall  particle  activity  to  the  left 
of  center  and  a  corresponding  alteration  of  the  quiescent 
zone  near  the  bottom  of  the  bed  such  that  the  majority  of 
stagnant  particles  were  located  in  the  right  hand  corner. 


70 


This  phenomenon  was  particularly  noticeable  in  the  wide  bed 
configurations.  An  example  of  this  off-center  activity  is 
shown  in  Figure  27. 

Examination  of  the  apparatus  revealed  that  a  stream  of 
air  was  being  injected  into  the  bed  from  the  area  between 
the  baffle  plate  and  the  distributor.  Air  from  the  inlet 
plenum  was  passing  through  the  unused  portion  of  the  dis¬ 
tributor  to  the  left  of  the  movable  wall.  Once  the  air  was 
through  the  distributor,  further  upward  progress  was  blocked 
by  the  baffle  plate  and  the  sealed  joints.  There  was  however, 
a  very  thin  horizontal  passage  between  the  baffle  and  the  top 
of  the  distributor  as  shown  in  Figure  28.  The  air  flowed 
along  this  path  until  it  emptied  into  the  bed  itself.  Because 
of  this  "peripheral  injection"  there  was  an  uneven  distribution 
of  fluidizing  air  and  a  resulting  shi£t  in  the  center  of  activ¬ 
ity  toward  the  left  hand  side.  In  order  to  determine  the  effect 
of  this  phenomenon  an  entire  set  of  data  was  taken  with  the 

unused  portion  of  the  distributor  masked-off  from  the  under- 

* 

side  so  that  no  air  could  enter.  During  these  masked  data 
runs  it  was  observed  that  the  bubble  path  was  very  well 
centered  as  was  the  quiescent  zone  near  the  bottom  of  the  bed. 
All  other  data  runs  were  taken  with  the  distributor  unmasked 
in  order  to  study  the  difference  in  heat  transfer  performance 
between  the  two  sides.  A  detailed  study  of  particle  cell 
flow  patterns  for  this  apparatus  was  presented  by  Morgan 
[Ref.  2]  in  1981. 


B.  PRESSURE  DROP  DATA 


Measurements  were  taken  of  the  pressure  drop  across  the 
bed  for  all  flow  rates  and  bed  configurations.  This  data 
can  be  found  in  Appendix  A.  Plots  were  constructed  of  LOG 
AP  vs.  LOG  G  for  each  bed  configuration  in  order  to  study 
the  variation  in  pressure  drop  with  changes  in  the  super¬ 
ficial  mass  velocity.  These  plots  are  shown  in  Figures  29 
through  32.  An  ideal  pressure  drop  curve  is  shown  in  Figure 
33.  Each  of  the  plots  display  the  same  basic  characteristic 
of  the  ideal  curve  in  that  the  pressure  drop  increases 
steadily  until  the  point  of  minimum  fluidization  (G  j)  is 
reached.  From  that  point  onward,  the  pressure  tends  to 
remain  relatively  constant  despite  further  increases  in  air 
flow  rate.  This  is  because  of  the  fact  that  once  the  bed  is 
fluidized,  any  further  increase  in  the  flow  rate  imparts  an 
increased  lift  to  the  particles  because  of  the  viscous  drag 
force  exerted  on  them  by  the  moving  air  stream.  As  the  lift 
increases,  expansion  of  the  bed  occurs  and  voidage  (e)  in¬ 
creases.  As  the  voidage  increases,  the  interstitial  spacing 
between  the  particles  also  increases.  This  results  in  a 
reduction  of  the  interstitial  gas  velocity,  which  in  turn, 
lowers  the  drag  exerted  on  the  particles  until  the  gravita¬ 
tional  force  and  the  lift  are  once  again  in  balance.  As  a 
result  of  the  increased  voidage,  the  pressure  drop  across 
the  bed  remains  constant  or  may  even  decrease. 
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Examination  of  the  pressure  drop  curve  for  the  6.25  inch 
configuration  (Fig.  29)  shows  a  deviation  from  ideal  behavior 
in  that  the  pressure  drops  sharply  after  reaching  the  point 
of  minimum  fluidization.  A  pressure  drop  of  this  sort  is 
normally  associated  with  channeling.  As  indicated  previously, 
no  visual  indications  of  channeling  were  detected,  however 
it  is  possible  that  "intermittent  channeling"  was  taking 
place.  This  type  of  channeling  would  be  difficult  to  differ¬ 
entiate  from  rapid  bubbling  because  the  channel  does  not 
extend  through  the  entire  depth  of  the  bed  as  it  does  in 
"through  channeling."  An  alternate  explanation  for  the  drop 
in  pressure  is  the  effect  of  particle  interlocking  which  some¬ 
times  manifests  itself  during  the  initial  fluidization  of  a 
bed.  It  is  noted  that  subsequent  data  runs  did  not  demonstrate 
such  behavior. 

Based  on  the  curve  shown  in  Figure  29  the  point  of  minimum 

fluidization  for  the  6.25  inch  configuration  occurred  at  G=159.6 
2 

Lbm/Hr-Ft  (261  flow).  It  was  was  noted  ,the  visual  observation 
of  fluidization  onset  occurred  at  the  same  air  flow  rate  as 
was  predicted  by  this  curve. 

The  pressure  drop  curve  for  the  8  inch  configuration 
(Fig.  30)  is  in  very  close  agreement  with  the  ideal  curve 
with  no  indications  of  channeling  or  slugging.  Based  on  this 

curve  the  point  of  minimum  fluidization  occurred  at  G=132.6 

2 

Lbm/Hr-Ft  (281  flow).  Fluidization  was  visually. observed 

2 

to  commence  at  137  Lbm/Hr-Ft  (29$  flow). 
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The  curve  for  the  10  inch  configuration  (Fig.  31)  was 
smoother  than  the  previous  configurations  and  matched  the 
ideal  curve  closely  except  for  a  gradual  decrease  in  pres¬ 
sure  following  fluidization  onset.  The  reason  for  this 
decrease  was  not  readily  apparent.  The  curve  shows  that 
fluidization  commenced  at  G=148  Lbm/Hr-Ft2  (40%  flow). 

The  observed  starting  point  of  fluidization  was  also  at 
40%  flow. 

The  12  inch  configuration  curve  (Fig.  32)  was  also  very 
smooth  and  similar  to  the  ideal.  Like  the  10  inch  curve, 
a  slight  decrease  in  differential  pressure  was  noted.  The 
starting  point  of  fluidization  (Gm£)  indicated  by  the  curve 
was  G=146.5  Lbm/Hr-Ft^  (48%  flow).  Fluidization  was  visually 
observed  to  commence  at  the  same  point. 

A  combined  graph  showing  the  pressure  drop  curves  for  all 
four  configurations  is  shown  in  Figure  34.  With  the  exception 
of  the  8  inch  configuration,  they  show  a  slight  gradual  de¬ 
crease  in  the  minimum  fluidization  superficial  mass  \ elocity 
(Gmf)  as  the  bed  width  is  increased.  A  similar  trend  was 
observed  by  Morgan  [Ref.  2]  who  attributed  this  to  a  reduction 
in  the  effect  of  the  wall  frictional  drag  as  a  fraction  of 
the  overall  drag  experienced  by  the  bed. 

The  average  bed  pressure  drop  at  minimum  fluidization 
( Pm£ )  for  all  configurations  was  approximately  16.0  in.  H-,0 
or  0.577  psi.  Rese.arch  by  Trivedi  et  al  [Ref.  6]  indicated 
that  the  pressure  drop  across  a  fluidized  bed  should  equal 
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the  weight  of  the  bed  material  per  unit  cross  sectional  area. 

The  settled  bed  density  for  the  glass  beads  used  in  this  ex-  • 

periment  was  96  Lbm/Ft^  which  resulted  in  a  weight  -  to  -  area 

2 

ratio  of  0.681  Lbm/In  .  Thus,  the  pressure  drop  observed 
was  slightly  less  than  the  theoretical  value.  Recalculation 
of  the  weight - to-area  ratio  based  on  the  expanded  bed  density 
yielded  essentially  the  same  results. 

C.  TEMPERATURE  DISTRIBUTION 

The  temperature  distribution  across  each  of  the  heater 
surfaces  was  examined  for  all  data  runs  and  two  distinct 
trends  were  found  to  be  evident.  First;  the  four  thermo¬ 
couple  readings  in  any  horizontal  row  were  found  to  be  in 
very  close  agreement  with  one  another  in  all  cases  as  shown 
in  Appendix  E.  This  even  distribution  of  temperature  hori¬ 
zontally  across  the  plate  suggested  that  the  local  rate  of 
heat  transfer  from  the  plate  was  the  same  for  any  given  point 
along  the  plate  surface  at  a  given  bed  depth.  This  conclusion 
was  based  on  the  assumption  that  the  heat  flux  into  the  plate 
was  also  horizontally  uniform.  The  second  trend  observed  was 
that  there  was  a  variation  in  heater  temperature  from  the  top 
of  the  plate  to  the  bottom  and  that  this  "temperature  profile" 
displayed  a  distinct  and  predictable  change  in  shape  as  the 
air  fluidization  flow  rate  was  increased.  Figures  35  through 
44  illustrate  these  trends.  At  low  air  flow  rates  the  point 
of  greatest  temperature  was  always  at  the  top  of  the  plate. 
Moving  downward,  the  temperature  then  gradually  dropped  in 


a  non-linear  but  smooth  fashion  to  a  level  approximately  three 
to  five  degrees  less  than  that  observed  at  the  top  of  the 
plate.  Examination  of  the  plots  of  temperature  vs  bed  depth 
for  the  low  flow  rate  runs  showed  all  of  the  curves  to  have 
the  same  shape  with  each  curve  displaying  a  positive  slope. 

It  was  noted  that  as  the  flow  rates  increased,  the  slope  of 
these  curves  became  more  positive  until  they  assumed  a  ver¬ 
tical  orientation  and  then  reversed  their  slope  to  become 
negative.  It  was  further  noted  that  this  slope  reversal 
occurred  at  or  near  the  point  of  fluidization  onset.  The 
change  in  slope  of  these  curves  meant  that  the  peak  tempera¬ 
ture  within  the  plate  was  moving  downward  from  the  top  of 
the  plate  toward  the  bottom  as  the  fluidization  air  flow  rate 
was  increased.  A  possible  explanation  of  such  behavior  lies 
in  the  fact  that  prior  to  fluidization  onset,  when  the  par¬ 
ticles  were  not  in  motion,  the  heater  was  primarily  cooled 
by  an  upward  flow  of  air  from  the  distributor.  Because  of 
this  upward  flow  path,  the  temperature  of  the  heater  tended 
to  increase  with  an  upward  movement  along  the  plate.  Once 
the  bed  was  fluidized,  however,  a  significant  amount  of 
cooling  was  then  provided  by  the  moving  particles  sweeping 
along  the  plate  surface.  Since  the  particles  were  moving 
downward  along  the  wall  a  thermal  gradient  was  introduced 
which  was  opposite  to  that  caused  by  the  upward  flow  of  air. 
The  peak  temperature  within  the  plate  would  consequently  be 
observed  at  the  point  where  these  two  gradients  met.  If  the 


cooling  due  to  the  particles  was  greater  than  that  provided 
by  the  upward  air  flow,  the  peak  temperature  would  be  found 
toward  the  bottom  of  the  plate.  Conversely,  if  the  air  flow 
cooling  predominated,  the  peak  temperature  would  be  seen  in 
the  upper  section  of  the  plate.  An  increase  in  the  fluidiza¬ 
tion  air  flow  rate  would  provide  more  upward  flowing  cooling 
air  while  simultaneously  increasing  the  downward  moving 
particle  activity. 

Experimental  results  showed  that  an  increase  in  air  flow 
rate  caused  the  peak  temperature  to  move  downward.  This  sug¬ 
gested  that,  as  the  flow  rate  was  increased,  the  increase  in 
particle  cooling  predominated  over  the  increase  in  air  cooling. 
The  reason  for  this  predominance  may  have  been  that  the  par¬ 
ticles  had  a  much  higher  specific  heat  capacity  than  the  air 
and  a  lower  mean  velocity.  Consequently,  they  were  able  to 
remove  a  larger  quantity  of  heat.  Another  possible  cause  for 
the  particle  cooling  predominance  as  flow  rates  were  increased 
was  that  the  motion  of  the  cooling  air  may  have  been  altered 
once  fluidization  commenced  so  that  more  of  the  upward  moving 
gas  was  concentrated  in  the  center  of  the  bed.  Although  there 
was  no  objective  evidence  that  this  was  occurring,  it  may  have 
been  that  any  increase  in  flow  rate  was  merely  diverted  to 
the  center  of  the  bed  and  consequently  did  not  increase  the 
amount  of  cooling  air  available  near  the  wall. 

Another  trend  which  was  observed  while  studying  the  heater 
temperature  distributions  was  that  as  the  bed  width  increased, 
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the  shift  in  the  orientation  of  the  temperature  curve  became 
less  pronounced.  In  the  6  inch  configuration  the  peak 
temperature  point  moved  from  the  top  of  the  bed  downward  to 
the  very  bottom  as  fluidization  activity  increased.  The 
maximum  temperature  differential  observed  between  the  top 
and  the  bottom  of  the  plate  was  approximately  17  degrees. 

In  the  8  inch  configuration,  the  peak  temperature  point 
also  transited  down  the  entire  length  of  the  plate,  however 
it  did  so  at  a  slower  rate  than  in  the  6  inch  configuration. 
The  maximum  top-to-bottom  differential  was  essentially  the 
same  as  that  seen  for  the  6  inch  bed.  With  the  bed  con¬ 
figured  to  the  10  inch  width,  the  peak  temperature  did  not 
move  down  the  entire  length  of  the  plate,  but  instead  stopped 
at  a  level  approximately  3  inches  above  the  bottom  edge.  The 
shift  in  temperature  orientation  was  even  more  gradual  than 
in  the  8  inch  configuration.  The  maximum  observed  tempera¬ 
ture  differential  was  approximately  11  degrees.  While  in 
the  12  inch  configuration,  very  little  change  was  observed 
in  the  temperature  distribution  as  the  flow  rate  was  in¬ 
creased.  The  peak  temperature  shifted  from  the  top  of  the 
plate  to  just  above  center,  however  the  maximum  temperature 
differential  of  6  degrees  was  significantly  less  than  that 
observed  in  the  other  configurations. 

In  general,  the  right  hand  heater  temperature  distribu¬ 
tion  proved  to  be  more  sensitive  to  changes  in  flow  rate 
than  did  the  left  hand  heater. 
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A  characteristic  unique  to  the  left  hand  heater  was  that 
immediately  prior  to  commencement  of  the  downward  transit  of 


the  peak  temperature  point,  the  temperature  differential 
between  the  top  and  the  bottom  of  the  plate  suddenly  increased 
to  approximately  three  times  its  previous  value.  As  the  flow 
rate  was  further  increased,  this  differential  abruptly  dropped 
back  down  to  its  former  level.  This  trend  is  shown  in  Figures 
40  through  43.  It  is  believed  that  this  fluctuation  could 
be  related  to  a  possible  reorientation  of  air  flow  patterns 
as  fluidization  commences.  This  phenomenon  was  not  evident 
in  the  right  hand  heater  temperatures. 

Examination  of  the  data  in  Appendix  C  shows  that  for  each 
flow  rate  and  bed  configuration,  the  average  temperature  of 
the  left  hand  heater  was  lower  than  for  the  right  hand  heater. 

In  many  cases,  the  difference  in  cooling  of  the  plates 
was  quite  significant.  During  the  data  runs  with  the  distri¬ 
butor  masked-off,  it  was  observed  that  the  left  and  right 
hand  heaters  were  being  cooled  to  approximately  the  same 
degree,  since  the  average  heater  temperatures  were  in  close 
agreement  with  one  another. 

i 

The  temperature  profile  of  the  bed  itself,  as  measured 
by  the  immersed  probes,  showed  a  non-uniform  distribution 
prior  to  the  commencement  of  fluidization  for  all  configura¬ 
tions.  The  temperatures  were  highest  near  the  heater  sur¬ 
faces  where  particle  heating  was  taking  place  by  conduction. 
Because  of  the  relatively  low  rate  of  heat  transfer  the 
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NOTE:  Numerical  values  represent  flowmeter  readings 

percentage  of  maximum  rotometer  capacity. 


Figure  36.  Right  Hand  Heater  Temperature  Profiles;  6 
Masked  Configuration 


Vertical  Position 


NOTE:  Numerical  values  represent  flowmeter  readings  as  a 
percentage  of  maximum  rotometer  capacity. 
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temperature  tended  to  decrease  toward  the  center.  Immediately 
upon  commencement  of  fluidization,  however,  the  temperature 
distribution  within  the  bed  became  extremely  uniform  and 
remained  so  as  fluidization  flow  rates  increased.  This  was 
directly  attributable  to  the  thorough  mixing  induced  by  the 
fluidized  particle  motion.  In  most  cases,  the  temperature 
of  the  outgoing  airstream  was  within  two  or  three  of  the 
temperature  of  the  bed  itself. 

D.  HEAT  TRANSFER  COEFFICIENTS 

In  order  to  establish  a  performance  baseline,  data  was 
collected  for  a  bed  from  which  all  of  the  particles  had  been 
removed  and  an  average  heat  transfer  coefficient  of  3.8  Btu/ 
Hr-Ft“-  F  was  calculated.  All  remaining  runs  were  accom¬ 
plished  with  the  apparatus  filled  to  a  settled  bed  depth  of 
approximately  12  inches.  The  resulting  heat  transfer  data 
for  all  runs  is  tabulated  in  Appendix  C  and  plots  of  heat 
transfer  coefficient  vs.  superficial  mass  velocity  are 
shown  in  Figures  45  through  S3.  Each  of  these  plots  demon¬ 
strates  the  same  basic  trend  in  that  initially,  during  low 
gas  flow  rate  conditions,  the  heat  transfer  coefficient  was 
correspondingly  low  (approximately  4.5  Btu/Hr-Ft“-0F)  and 
remained  low  until  a  flow  rate  approaching  the  minimum 
fluidization  superficial  mass  velocity  (Gm£)  was  attained. 

From  this  point  onward,  the  heat  transfer  coefficient  in¬ 
creased  in  a  roughly  linear  fashion  as  the  gas  flow  rate  was 
increased. 
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Figure  47.  Heat  Transfer  Coefficient  vs  Superficial  Mass 
Velocity;  Right  Hand  Heater,  6.25  in. 
Configuration 
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Figure  48.  Heat  Transfer  Coefficient  vs  Superficial  Mass 
Velocity;  Left  Hand  Heater,  6.25  in. 
Configuration 
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Figure  49.  Heat  Transfer  Coefficient  vs  Superficial  Mass 
Velocity;  8.0  in.  Configuration 
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A  second  significant  trend  which  was  observed  was  the 
fact  that  the  heat  transfer  coefficient  for  the  left  hand 
heater  was  significantly  higher  than  for  the  right  hand 
heater,  particularly  as  the  bed  width  was  increased.  In 
the  wide  bed  configurations  the  heat  transfer  coefficient 
for  the  right  hand  heater  showed  very  little  improvement 
in  response  to  an  increase  in  flow  rate.  With  a  high 
height-to-width  ratio  (narrow  bed),  the  heat  transfer 
coefficient  for  the  left  hand  heater  was  roughly  twice 
the  value  for  the  right  hand  heater  at  high  gas  flow  rates. 
With  a  low  height-to-width  ratio  (wide  bed)  ,  the  left  hand 
heat  transfer  coefficient  exceeded  that  for  the  right  by 
a  factor  of  15  to  1. 

The  exception  to  this  trend  occurred  during  the  data 
runs  accomplished  with  the  distributor  masked-off.  Under 
this  condition,  the  heat  transfer  coefficients  for  the  left 
hand  and  right  hand  heaters  were  in  much  closer  agreement 
with  one  another  than  in  all  other  cases.  The  effect  of 
masking  the  distributor  was  to  increase  the  heat  transfer 
coefficient  for  the  right  hand  heater  by  approximately  10% 
while  simultaneously  reducing  the  heat  transfer  coefficient 
for  the  left  hand  heater  by  approximately  25%.  It  would  be 
expected  that,  with  the  distributor  masked-off,  the  heat 
transfer  coefficients  would  be  the  same  for  both  the  left 
hand  and  right  hand  heaters.  Although  there  is  a  much 
closer  match  than  in  the  unmasked  configuration,  it  is 
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observed  that  the  left  hand  heater  still  achieved  higher  rates 
of  heat  transfer.  There  are  several  possible  explanations 
for  this  difference.  The  first  of  these  is  that  the  unused 
portion  of  the  distributor  may  not  have  been  completely  sealed 
by  the  masking  and  as  a  result,  some  peripheral  injection 
could  have  still  have  been  taking  place.  Secondly,  it  was 
noted  that  there  was  a  small  amount  of  gas  leakage  past  the 
movable  side  wall  gaskets  which  may  have  altered  the  gas  and 
particle  flow  patterns  along  the  left  wall.  This,  in  turn, 
may  have  altered  the  heat  transfer  coefficients.  Thirdly, 
it  is  possible  that  the  movable  side  wall  may  not  have  been 
perfectly  vertical.  Although  care  was  taken  to  accurately 
position  and  secure  the  side  wall,  deviations  from  a  true 
vertical  orientation  of  as  much  as  1°  are  considered  possible. 
Research  conducted  in  the  Soviet  Union  by  Filippovskii  and 
Baskakov  [Ref.  7]  has  indicated  that  plate  angle  has  a  con¬ 
siderable  influence  on  heat  transfer  coefficients.  Because 
of  the  small  angle  involved  in  this  case,  the  contribution 
of  this  effect  is  not  considered  to  be  the  principle  cause 
of  the  difference  in  performance  between  the  two  plates. 

A  final  effect  which  may  have  made  a  minor  contribution  to 
this  difference  is  the  fact  that  the  losses  from  the  heater 
backing  to  the  atmosphere  were  lower  for  the  left  hand  heater 
than  for  the  right.  As  a  result,  more  of  the  electrical 
energy  going  into  the  heater  was  directed  through  the  copper 
plate  and  on  to  the  bed,  thus  resulting  in  a  minor  improvement 


in  the  heat  transfer  coefficient.  Because  of  the  configura- 
tion  of  the  apparatus,  the  left  hand  heater  mounting  block 
was  enclosed  within  the  unused  portions  of  the  apparatus 
walls.  Because  of  this,  there  was  a  stagnation  of  the  air 
in  the  vicinity  of  the  mounting  block  as  well  as  an  elevation 
of  the  local  ambient  temperature.  These  effects  combined  to 
reduce  the  amount  of  heat  lost  through  the  mounting  block. 
Because  of  the  specialized  conditions  which  were  found  to 
exist  near  the  left  hand  wall,  the  performance  of  the  right 
hand  heater  shall  be  considered  as  being  representative  of 
a  typical  flat  plate. 

The  magnitude  of  the  heat  transfer  coefficients  calculated 
as  a  result  of  this  study  were  within  the  range  of  values 
found  in  the  literature,  however  it  should  be  noted  that 
there  is  a  wide  variation  in  the  coefficients  found  by 
other  researchers  due  to  the  effect  of  differing  bed  sizes 
and  geometries,  particle  size  and  composition,  fluidizing 
gas  properties,  superficial  mass  velocities,  and  bed  operat¬ 
ing  temperatures. 

There  is  also  a  wide  variety  of  empirical  correlations 
available  for  the  prediction  of  heat  transfer  coefficients, 
however  each  of  these  serves,  for  the  most  part,  only  to 
describe  the  particular  test  apparatus  and  conditions  from 
which  they  were  derived.  For  this  reason,  no  attempt  was 
made  to  fit  the  data  to  the  correlations  currently  available 


nor  was  an  attempt  made  to  develop  a  new  correlation  to 
describe  the  performance  of  the  apparatus  used  in  this 
experiment . 

A  review  of  the  literature  shows  that  numerous  researchers 
have  observed  heat  transfer  curves  similar  to  that  shown  in 
Figure  54.  This  typical  curve  displays  a  peak  heat  transfer 
coefficient  (h  )  at  some  optimum  superficial  mass  velocity 
(GQ  ) .  Flow  rates  above  this  optimum  value  resulted  in  a 
reduction  in  the  heat  transfer  coefficient.  This  has  been 
attributed  to  the  fact  that  as  the  voidage  is  increased, 
there  is  a  point  at  which  the  beneficial  contribution  of 
increased  particle  activity  is  offset  by  the  negative  effect 
of  reducing  the  solids  concentration  within  the  emulsion. 

In  this  study  such  an  optimum  flow  rate  was  not  achieved. 

It  is  believed  that,  due  to  the  design  of  the  fluidization 
apparatus,  the  range  of  gas  flow  rates  was  restricted  by 
the  relatively  shallow  "freeboard"  region  above  the  level 
of  the  settled  bed.  .  When  the  bed  was  undergoing  active 
fluidization,  particles  within  the  "splash  zone"  reached 
the  top  of  the  container  and  blocked  the  air  exit  ports. 

This  was  believed  to  be  taking  place  at  gas  flow  rates 
below  the  optimum  level.  It  is  possible  that  a  reduction 
in  the  settled  bed  depth  and  a  consequent  increase  in  the 
height  of  the  freeboard  region  would  permit  higher  super¬ 
ficial  mass  velocities  and  result  in  higher  levels  of  heat 


transfer . 


Heat  Transfer  Coefficient  vs  Superficial  Mass 
Velocity;  Typical 
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The  difference  in  the  level  of  heat  transfer  between 
the  left  and  right  hand  walls  is  clearly  related  to  the 
high  level  of  cooling  observed  at  the  left  hand  wall. 

This,  in  turn,  is  related  to  the  large  amount  of  particle 
activity  observed  to  the  left  of  center.  The  exact  mech¬ 
anism  of  heat  transfer  which  takes  place  at  the  wall  cannot 
be  stated  with  certainty.  Although  there  are  many  theories 
as  to  the  mechanism  of  heat  transfer  within  a  fluidized 
bed,  the  bulk  of  the  models  postulated  to  date  fall  into 
one  of  two  broad  categories.  The  first  of  these  models 
proposes  that  the  heat  transfer  is  primarily  through  a  gas 
film  adjacent  to  the  heat  transfer  surface.  The  heat  is 
transferred  by  steady  state  conduction  through  the  film  to 
the  downward  moving  particles.  This  film  is  kept  thin  by 
the  scrubbing  action  of  the  solid  particles  which  come  in 
contact  with  the  surface.  This  reduction  in  the  boundary 
layer  thickness  is  thought  to  be  the  key  to  achieving  high 
rates  of  heat  transfer.  The  second  model  theorizes  that 
heat  transfer  is  primarily  by  unsteady  conduction  to  "packets" 
which  are  made  up  of  several  solid  particles  grouped  together 
for  a  short  period  of  time.  These  packets  are  thought  to 
be  in  contact  with  the  heat  transfer  surface  for  only  brief 
periods  of  time  before  they  move  away  from  the  surface  and 
are  replaced  or  "replenished"  by  fresh  packets.  As  a  result, 
there  is  always  a  high  temperature  differential  between  the 
surface  and  the  packets  and  this,  in  turn,  yields  high  levels 
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of  heat  transfer.  A  significant  factor  shared  by  both 
the  packet  model  and  the  film  conduction  model  is  that 
they  both  rely  on  the  free  motion  of  particles  to  effect 
high  rates  of  heat  transfer.  It  is  believed  that  in  this 
apparatus  a  high  level  of  particle  activity  near  the  wall 
promotes  the  frequent  exchange  of  particles  from  within  the 
bed  core  where  it  is  comparatively  cool.  As  the  frequency 
of  particle  or  packet  replenishment  increases,  the  tempera¬ 
ture  differential  between  the  wall  and  the  particles  also 
increases.  This,  in  turn,  increases  the  overall  heat  transfer 
Conversely,  when  there  is  comparatively  little  particle  activ¬ 
ity  near  the  wall,  the  replenishment  rate  is  low  and  the  same 
particles  remain  in  contact  with  the  hot  wall  for  longer 
periods  of  time  as  they  move  downward.  Because  of  the  long 
residence  time,  the  temperature  differential  is  diminished 
and  the  heat  transfer  coefficient  is  reduced.  In  this  study, 
the  level  of  particle  activity  near  the  right  hand  wall  was 
comparatively  low  and  so  was  the  heat  transfer  coefficient. 

A  study  of  the  temperature  profiles  for  the  right  hand  plate 
shows  a  steady  increase  in  plate  temperature  with  a  downward 
change  in  vertical  position.  This  profile  is  consistent 
with  that  expected  for  long  particle  residence  times  since 
the  downward  moving  particles  would  experience  a  steady  in¬ 
crease  in  temperature  as  they  approached  the  bottom  of  the 
plate.  It  is  noted  that  when  the  configuration  of  the  bed 
is  tall  and  narrow,  the  bubble  flow  path  is  very  close  to 


the  bed  walls.  This  causes  a  high  degree  of  particle  activity 
near  the  wall  and  a  corresponding  increase  in  the  particle 
replenishment  rate.  As  the  bed  height - to-width  ratio  is 
decreased  by  increasing  the  bed  width,  the  bubble  flow  path 
moves  further  away  from  the  wall.  This  causes  a  reduction 
in  the  amount  of  activity  near  the  wall  and  a  corresponding 
decrease  in  the  replenishment  rate.  This  scenario  is  con¬ 
sistent  with  the  observed  reduction  in  heat  transfer  coef¬ 
ficient  as  the  bed  width  was  increased.  Conversely,  the 
amount  of  activity  seen  by  the  left  hand  plate  was  consis¬ 
tently  high  as  a  result  of  the  peripheral  injection  effect 
described  previously.  The  degree  of  activity  and  consequently 
the  heat  transfer  coefficient,  remained  unaffected  by  an  in¬ 
crease  in  bed  width.  In  addition,  the  heater  temperature 
profile  showed  no  sharp  transition  from  top  to  bottom,  thus 
indicating  that  the  particle  residence  time  was  short. 

Based  on  the  findings  of  this  study,  it  may  be  concluded 
that  peripheral  injection  is  an  effective  means  of  promoting 
particle  activity  or  turbulence  along  the  wall  of  the  con¬ 
tainer.  The  idea  of  turbulence  promoters  in  fluidized  beds 
is  not  a  new  one.  Early  research  by  the  Babcock  and  Wilcox 
Company  [Ref.  8]  included  the  use  of  spirally  shaped  strips 
to  create  a  rotational  motion  which  forced  the  solids  toward 
the  walls  by  centrifugal  force.  The  result  was  a  significant 

increase  in  heat  transfer.  Recently  Sokolov  et  al  [Ref.  9] 

% 

in  the  Soviet  Union  studied  the  use  of  angled  louvers 
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positioned  within  the  bed  in  order  to  redirect  the  motion  of 
the  particles  outward  toward  the  walls,  thus  destroying  the 
stagnant  layer  which  was  insulating  the  surface.  The  result 
was  a  marked  improvement  in  heat  transfer  performance. 
Although  more  study  is  obviously  warranted,  the  use  of  pe¬ 
ripheral  injection  would  appear  to  offer  benefits  similar 
to  other  enhancement  techniques  while  eliminating  the  need 
for  in-bed  structures  which  occupy  valuable  space  and  which 
require  periodic  replacement  as  a  consequence  of  particle 
errosion.  In  addition,  the  amount  of  turbulence  created  at 
the  wall  by  the  injection  technique  could  easily  be  adjusted 
to  optimum  levels  by  controlling  the  flow  rate  of  the  in¬ 
jection  gas.  Guch  a  system  could  be  used  to  offset  the  low 
rates  of  heat  transfer  at  the  walls  of  wide  fluidized  beds. 

As  discussed  previously,  an  energy  balance  was  conducted 
in  order  to  confirm  the  validity  of  the  numerical  results. 
The  results  of  this  energy  balance  were  disappointingly 
poor.  In  most  cases,  the  measured  heat  flux  leaving  the 
bed  with  the  outgoing  airstream  was  found  to  be  only  half 
of  that  calculated  to  be  entering  the  bed  via  the  electrical 
heaters.  It  is  believed  that  the  disagreement  was  primarily 
due  to  incorrect  measurement  of  the  fluidizing  gas  flow 
rate.  Despite  the  poor  energy  balance,  the  conclusions 
reached  as  a  result  of  this  study  are  considered  to  be  valid 
since  they  were  based  on  observed  trends  rather  than  on  the 
magnitude  of  specific  numerical  results. 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

Based  on  the  experimental  results  of  this  study  the 
following  conclusions  may  be  drawn: 

1.  The  plate*to-bed  heat  transfer  coefficient  for  a  flat 
vertical  plate  serving  as  a  containing  wall  decreases  with 
a  decrease  in  the  bed  he ight - to -width  ratio. 

2.  The  plate-to-bed  heat  transfer  coefficient  increases 
with  an  increase  in  superficial  mass  velocity  above  the  mini¬ 
mum  required  to  initiate  fluidization. 

3.  The  plate-to-bed  heat  transfer  coefficient  remains  small 
and  relatively  constant  at  superficial  mass  velocities  below 
the  minimum  required  to  initiate  fluidization. 

4.  The  plate-to-bed  heat  transfer  coefficient  increases 
with  an  increase  in  the  particle  activity  or  turbulence  in 
the  vicinity  of  the  wall. 

5.  The  particle  turbulence  near  the  vertical  plate,  and 
consequently  the  heat  transfer  coefficient,  may  be  enhanced 
by  the  injection  of  air  from  beneath  the  plate  in  a  direction 
perpendicular  to  the  primary  fluidizing  air  flow  path. 

6.  The  temperature  profile  of  a  vertical  flat  plate  may  be 
used  to  provide  insights  as  to  the  nature  of  particle  activity 
in  the-vicinity  of  the  plate. 
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7.  The  particle  flow  pattern  within  a  rectangular  fluidized 
bed  is  in  very  close  agreement  with  that  described  in  the 
literature  for  cylindrical  fluidized  beds  with  upward  motion 
of  the  particles  taking  place  along  the  centerline  of  the 
container  and  downward  return  flow  of  the  particles  taking 
place  along  the  periphery. 

8.  Slugging  or  heaving  of  the  bed  was  more  prevalent  in  the 
narrow  bed  configuration  as  opposed  to  the  wide  bed  configura 
tion  where  this  phenomenon  was  not  observed. 

9.  Fluidized  bed  particle  motion  along  a  non-moving  boundary 
matches  the  behavior  of  a  viscous  fluid  in  that  a  "no-slip" 
condition  is  observed. 

10.  The  bed  pressure  drop  characteristics  for  a  rectangular 
fluidized  bed  match  those  found  in  the  literature  for  a  cy¬ 
lindrical  fluidized  bed  in  that  the  pressure  drop  across  the 
bed  increases  linearly  with  an  increase  in  air  flow  rate 
until  the  point  of  fluidization  onset.  Once  fluidization 
has  commenced,  the  pressure  drop  remains  constant. 

11.  The  temperature  of  the  heater  surface  at  a  ’given  depth 
remains  constant  for  any  point  in  a  horizontal  direction. 

12.  The  bed  temperature  will  assume  a  uniform  temperature 
distribution  once  fluidized. 

B.  RECOMMENDATIONS  FOR  FURTHER  STUDIES 

In  order  to  better  pursue  the  ojbectives  of  this  study  at 
sometime  in  the  future,  the  following  recommendations  and/or 
suggestions  are  made: 


1.  Replace  the  three  circular  air  exhaust  tubes  on  the  top 
of  the  fluidization  apparatus  with  a  rectangular  exhaust 
port  which  runs  the  entire  twenty-four  inch  length  of  the 
bed.  This  would  provide  a  non- restrict ive  exit  path  for 
the  fluidizing  air  which  would  minimize  the  amount  of  swirl¬ 
ing  and  mixing  in  the  region  above  the  bed. 

2.  Repair  severed  wire  on  thermocouple  *15. 

3.  Replace  malfunctioning  thermocouples  *4  and  *40. 

4.  Install  a  minimum  of  four  new  thermocouples  to  monitor 
air  exit  temperature.  At  present  only  one  thermocouple  can 
be  used  for  this  purpose  (T(55))  for  bed  configurations  of 
12  inches  or  less. 

5.  Provide  a  more  positive  means  of  securing  the  baffle 
plate  to  the  distributor  plate  so  as  to  prevent  leak-by  into 
the  bed. 

6.  Install  a  bypass  valve  around  the  existing  air  inlet  con¬ 
trol  valve  to  permit  fine  adjustment  of  fluidization  flow  rates. 

7.  Install  two  additional  thermocouples  in  the  left  hand 
movable  wall  1/4  inch  below  the  bottom  of  the  heater  in 
order  to  obtain  a  more  accurate  estimate  of  wall  losses. 

8.  Install  two  additional  thermocouples  in  the  front  wall 
of  the  bed  so  as  to  obtain  a  more  accurate  estimate  of  wall 
losses . 

9.  Substitute  the  existing  power  supplies  with  units,  which 
are  capable  of  providing  additional  current  and  which  are 

t 

equipped  with  controls  to  vary  power  as  well  as  voltage. 
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This  would  permit  equalization  of  input  power  to  the  two 
heaters 'and  allow  the  desired  temperature  differential 
between  the  plate  and  the  bed  to  be  maintained. 

10.  Install  a  movable  wall  in  the  air  plenum  beneath  the 
distributor  plate  so  that  the  flow  path  followed  by  the  in¬ 
coming  air  is  identical  for  the  left  and  right  sides. 

11.  Prior  to  the  commencement  of  future  data  runs,  apply 
caulking  material  to  the  movable  side  wall  so  as  to  prevent 
leakage  of  fluidizing  air  past  the  gaskets. 

12.  Modify  the  fluidization  apparatus  sc  that  a  controlled 
amount  of  injection  air  may  be  applied  to  each  of  the  heater 
plates.  Such  a  modification  would  permit  additional  study 
of  the  effects  of  the  peripheral  injection  phenomenon. 

C.  CLOSING  REMARKS 

Gas  fluidization  is  a  fascinating  process  which  promises 
to  provide  a  multitude  of  beneficial  applications.  Full 
enjoyment  of  the  potential  which  exists  depends  on  a  mere 
complete  understanding  of  the  many  complex  phenomenon 
associated  with  fluidization.  It  is  hoped  that  this  thesis 
has,  in  some  small  way,  contributed  to  that  understanding 
or  helped  to  stimulate  additional  thought  and  research  in 
this  exciting  field  of  study. 
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Table  I.  Bed  Pressure/Gas  Muss  Flow  Data  (6.25  in.  Bed) 
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Table  II.  Bed  Pressure/Gas  Mass  Plow  Data  (8. 
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Table  III.  Bed  Pressure/Gas  Mass  Flow  Data  (10.0  in.  Bed) 
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Table  IV.  Bed  Pressure/Gas  Mass  Flow  Data  (12.0  in.  Bed) 
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APPENDIX  B. 

EXPERIMENTAL  ROTOMETER  CALIBRATION  DATA 


**ROTOMETER  CALIBRATION** 
ASME  HERSCHEL-TYPE  VENTURI 
4.26/2.13  IN  DIA 


20  PERCENT  READING 


**EXPER I MENTAL  DATA** 


ATMOSPHERIC  PRESSURE  »  29.76910  IN  Hg 
VENTURI  INLET  TEMPERATURE  -  72.0  DEG  F 
VENTURI  INLET  STATIC  GAGE  PRESSURE  -  10.70  IN  H20 
VENTURI  DIFFERENTIAL  PRESSURE  *  .010  IN  RED  OIL 
SPECIFIC  GRAVITY  OF  RED  OIL  -  .834 


* t CALCULATED  VALUES** 


ATMOSPHERIC  PRESSURE  « 

VENTURI  INLET  GAGE  PRESSURE  - 
VENTURI  INLET  ABSOLUTE  PRESSURE 
VENTURI  DIFFERENTIAL  PRESSURE  - 
VENTURI  INLET  AIR  DENSITY  - 
VENTURI  PRESSURE  DROP  RATIO  (X) 
EXPANSION  FACTOR  <Y>  - 
ASSUMED  REYNOLDS  NUMBER  - 
DISCHARGE  COEFFICIENT  <C>  - 
VELOCITY  OF  APPROACH  FACTOR  <E> 
FLOW  COEFFICIENT  <K>  - 
AIR  MASS  FLOW  RATE  - 
AIR  VOLUMETRIC  FLOW  RATE  - 
VENTURI  INLET  REYNOLDS  NUMBER  - 


14.6239  PSI 
.38654  PSI 
14.23932  PSI 
.00030  PSI 
.07230  LBM/FT3 
. 00002 
1.00 
3000 
.889 
1 . 03280 
.91816 
.01021  LBM/SEC 
8.469  CFM 
3002 


125 


**ROTOMETER  CALIBRATION** 
ASME  HERSCHEL-TYPE  VENTURI 
4.26/2. 13  IN  DIA 

25  PERCENT  READING 


♦•EXPERIMENTAL  DATA** 


ATMOSPHERIC  PRESSURE  *  29.76910  IN  Hg 
VENTURI  INLET  TEMPERATURE  *  72.0  DEG  F 
VENTURI  INLET  STATIC  GAGE  PRESSURE  -  12.35  IN  H20 
VENTURI  DIFFERENTIAL  PRESSURE  -  .011  IN  RED  OIL 
SPECIFIC  GRAVITY  OF  RED  OIL  -  .834 


♦♦CALCULATED  VALUES** 


ATMOSPHERIC  PRESSURE  - 
VENTURI  INLET  GAGE  PRESSURE  - 
VENTURI  INLET  ABSOLUTE  PRESSURE  * 
VENTURI  DIFFERENTIAL  PRESSURE  - 
VENTURI  INLET  AIR  DENSITY  * 
VENTURI  PRESSURE  DROP  RATIO  <X>  * 
EXPANSION  FACTOR  (Y)  * 

ASSUMED  REYNOLDS  NUMBER  - 
DISCHARGE  COEFFICIENT  <C>  « 
VELOCITY  OF  APPROACH  FACTOR  (E)  - 
FLOW  COEFFICIENT  (K>  =* 

AIR  MASS  FLOW  RATE  ■ 

AIR  VOLUMETRIC  FLOW  RATE  » 

VENTURI  INLET  REYNOLDS  NUMBER  » 


14.6259  PSI 
.44614  PSI 
14. 17972  PSI 
.00033  PSI 
.07200  LBM/FT3 
. 00002 
1.00 
3150 
.891 
1 . 03280 
.91970 
.01070  LBM/SEC 
8.916  CFM 


3147 


**ROTOMETER  CALIBRATION** 
ASME  HERSCHEL-TYPE  VENTURI 
4.26/2. 13  IN  DIA 


30  PERCENT  READING 


**EXPER I MENTAL  DATA** 


ATMOSPHERIC  PRESSURE  »  29.76910  IN  Hg 
VENTURI  INLET  TEMPERATURE  *  72.0  DEG  F 
VENTURI  INLET  STATIC  GAGE  PRESSURE  -  14.75  IN  H20 
VENTURI  DIFFERENTIAL  PRESSURE  -  .015  IN  RED  OIL 
SPECIFIC  GRAVITY  OF  RED  OIL  »  .334 


* *CALCULATED  VALUES** 


ATMOSPHERIC  PRESSURE  - 
VENTURI  INLET  GAGE  PRESSURE  * 
VENTURI  INLET  ABSOLUTE  PRESSURE 
VENTURI  DIFFERENTIAL  PRESSURE  « 
VENTURI  INLET  AIR  DENSITY  * 
VENTURI  PRESSURE  DROP  RATIO  <X> 
EXPANSION  FACTOR  <Y>  - 
ASSUMED  REYNOLDS  NUMBER  » 
DISCHARGE  COEFFICIENT  <C>  - 
VELOCITY  OF  APPROACH  FACTOR  <E> 
FLOW  COEFFICIENT  <K>  * 

AIR  MASS  FLOW  RATE  » 

AIR  VOLUMETRIC  FLOW  RATE  ■ 
VENTURI  INLET  REYNOLDS  NUMBER 


14.6259  PSI 
.53284  PSI 
14.09302  PSI 
.00045  PSI 
.07156  LBM/FT3 
. 00003 
1.00 
3700 
.895 
1 . 03280 
. 92435 
.01252  LBM/SEC 
10.497  CFM 
3682 
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**ROTOMETER  CALIBRATION** 
ASME  HERSCHEL-TYPE  VENTURI 
4.26/2. 13  IN  DIA 

35  PERCENT  READING 


**EXPER I MENTAL  DATA** 


ATMOSPHERIC  PRESSURE  ■  29.76910  IN  Hg 
VENTURI  INLET  TEMPERATURE  *  72.0  DEG  F 
VENTURI  INLET  STATIC  GAGE  PRESSURE  *  17.50  IN  H20 
VENTURI  DIFFERENTIAL  PRESSURE  »  .020  IN  RED  OIL 
SPECIFIC  GRAVITY  OF  RED  OIL  »  .834 


♦♦CALCULATED  VALUES** 


ATMOSPHERIC  PRESSURE  - 
VENTURI  INLET  GAGE  PRESSURE  - 
VENTURI  INLET  ABSOLUTE  PRESSURE 
VENTURI  DIFFERENTIAL  PRESSURE  - 
VENTURI  INLET  AIR  DENSITY  » 
VENTURI  PRESSURE  DROP  RATIO  (X) 
EXPANSION  FACTOR  <Y)  « 

ASSUMED  REYNOLDS  NUMBER  - 
DISCHARGE  COEFFICIENT  <C>  » 
VELOCITY  OF  APPROACH  FACTOR  <E> 
FLOW  COEFFICIENT  <K>  - 
AIR  MASS  FLOW  RATE  - 
AIR  VOLUMETRIC  FLOW  RATE  - 
VENTURI  INLET  REYNOLDS  NUMBER  - 


14.6259  PSI 
.63219  PSI 
13.99367  PSI 
.00060  PSI 
.07105  LBM/FT3 
. 00004 
1.00 
4250 
.899 
1 . 03280 
. 92797 
.01446  LBM/SEC 
12.211  CFM 
4253 


1 
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**ROTOMETER  CALIBRATION** 
ASME  HERSCHEL-TYPE  VENTURI 
4. 26/2. 13  IN  DIA 


40  PERCENT  READING 


**EXPER I MENTAL  DATA** 


ATMOSPHERIC  PRESSURE  -  29.76910  IN  Hg 
VENTURI  INLET  TEMPERATURE  =  72.0  DEG  F 
VENTURI  INLET  STATIC  GAGE  PRESSURE  *  20.30  IN  H20 
VENTURI  DIFFERENTIAL  PRESSURE  »  .025  IN  RED  OIL 
SPECIFIC  GRAVITY  OF  RED  OIL  -  .834 


*  * CALCULATED  VALUES** 


ATMOSPHERIC  PRESSURE  * 

VENTURI  INLET  GAGE  PRESSURE  » 
VENTURI  INLET  ABSOLUTE  PRESSURE 
VENTURI  DIFFERENTIAL  PRESSURE  * 
VENTURI  INLET  AIR  DENSITY  - 
VENTURI  PRESSURE  DROP  RATIO  <X> 
EXPANSION  FACTOR  <Y)  * 

ASSUMED  REYNOLDS  NUMBER  * 
DISCHARGE  COEFFICIENT  (C)  » 
VELOCITY  OF  APPROACH  FACTOR  (E> 
FLOW  COEFFICIENT  <K)  » 

AIR  MASS  FLOW  RATE  - 

AIR  VOLUMETRIC  FLOW  RATE  - 

VENTURI  INLET  REYNOLDS  NUMBER  =» 


14.6259  PSI 
.73334  PSI 
13.89252  PSI 
.00075  PSI 
.07054  LBM/FT3 
. 00005 
1.00 
4750 
.901 


1 . 03280 
. 93086 
.01616  LBM/SEC 
13.744  CFM 
4753 


**ROTOMETER  CALIBRATION** 
ASME  HERSCHEL -TYPE  VENTURI 
4.26/2. 13  IN  DIA 


45  PERCENT  READING 


**EXPER I MENTAL  DATA** 


ATMOSPHERIC  PRESSURE  =  29.76910  IN  Hg 
VENTURI  INLET  TEMPERATURE  *  72.0  DEG  P 
VENTURI  INLET  STATIC  GAGE  PRESSURE  -  24.30  IN  H20 
VENTURI  DIFFERENTIAL  PRESSURE  *  .031  IN  RED  OIL 
SPECIFIC  GRAVITY  OF  RED  OIL  «  .834 


♦♦CALCULATED  VALUES** 


ATMOSPHERIC  PRESSURE  » 

VENTURI  INLET  GAGE  PRESSURE  - 
VENTURI  INLET  ABSOLUTE  PRESSURE  * 
VENTURI  DIFFERENTIAL  PRESSURE  - 
VENTURI  INLET  AIR  DENSITY  * 
VENTURI  PRESSURE  DROP  RATIO  <X>  * 
EXPANSION  FACTOR  (Y)  * 

ASSUMED  REYNOLDS  NUMBER  = 
DISCHARGE  COEFFICIENT  <C>  * 
VELOCITY  OF  APPROACH  FACTOR  <E>  » 
FLOW  COEFFICIENT  <K)  - 
AIR  MASS  FLOW  RATE  * 

AIR  VOLUMETRIC  FLOW  RATE  » 

VENTURI  INLET  REYNOLDS  NUMBER  * 


14.6259  P SI 
.87784  PS I 
13.74802  PS I 
.00093  PS I 
.06981  LBM/FT3 
. 00007 
1.00 
5300 
.  904 
1 . 03280 
. 93365 
.01795  LBM/SEC 
15.431  CFM 
5281 


**ROTOMETER  CALIBRATION** 
ASME  HERSCHEL-TYPE  VENTURI 
4.26/2.13  IN  DIA 

50  PERCENT  READING 


** EXPERIMENTAL  DATA** 


ATMOSPHERIC  PRESSURE  »  29.76910  IN  Hg 
VENTURI  INLET  TEMPERATURE  *  72.0  DEG  F 
VENTURI  INLET  STATIC  GAGE  PRESSURE  »  27.30  IN  H20 
VENTURI  DIFFERENTIAL  PRESSURE  »  .037  IN  RED  OIL 
SPECIFIC  GRAVITY  OF  RED  OIL  »  .834 


* (CALCULATED  VALUES** 


ATMOSPHERIC  PRESSURE  » 

VENTURI  INLET  GAGE  PRESSURE  » 
VENTURI  INLET  ABSOLUTE  PRESSURE 
VENTURI  DIFFERENTIAL  PRESSURE  » 
VENTURI  INLET  AIR  DENSITY  » 
VENTURI  PRESSURE  DROP  RATIO  <X) 
EXPANSION  FACTOR  ( Y)  » 

ASSUMED  REYNOLDS  NUMBER  - 
DISCHARGE  COEFFICIENT  (C)  - 
VELOCITY  OF  APPROACH  FACTOR  <E> 
FLOW  COEFFICIENT  <K>  - 
AIR  MASS  FLOW  RATE  - 
AIR  VOLUMETRIC  FLOW  RATE  « 
VENTURI  INLET  REYNOLDS  NUMBER  « 


14.6259  PSI 
.98621  PSI 
13.63965  PSI 
.00111  PSI 
.06926  LBM/FT3 
. 00008 
1.00 
5750 
.906 


1 . 03280 
. 93592 
.01958  LBM/SEC 
16.967  CFM 
5760 


**ROTOMETER  CALIBRATION** 
ASME  HERSCHEL-TYPE  VENTURI 
4.26/2. 13  IN  DIA 

35  PERCENT  READING 


**EX PER I MENTAL  DATA** 


ATMOSPHERIC  PRESSURE  *  29.76910  IN  Hg 
VENTURI  INLET  TEMPERATURE  *  72.0  DEG  F 
VENTURI  INLET  STATIC  GAGE  PRESSURE  «  31.85  IN  H20 
VENTURI  DIFFERENTIAL  PRESSURE  =  .044  IN  RED  OIL 
SPECIFIC  GRAVITY  OF  RED  OIL  -  .834 


♦♦CALCULATED  VALUES** 


ATMOSPHERIC  PRESSURE  * 

VENTURI  INLET  GAGE  PRESSURE  - 
VENTURI  INLET  ABSOLUTE  PRESSURE  * 
VENTURI  DIFFERENTIAL  PRESSURE  » 
VENTURI  INLET  AIR  DENSITY  * 
VENTURI  PRESSURE  DROP  RATIO  (X)  * 
EXPANSION  FACTOR  <Y)  - 
ASSUMED  REYNOLDS  NUMBER  ■ 
DISCHARGE  COEFFICIENT  <C>  » 
VELOCITY  OF  APPROACH  FACTOR  (E)  - 
FLOW  COEFFICIENT  <K)  - 
AIR  MASS  FLOW  RATE  - 
AIR  VOLUMETRIC  FLOW  RATE  * 

VENTURI  INLET  REYNOLDS  NUMBER  » 


14.6259  PSI 
1. 15058  PSI 
13.47528  PSI 
.00133  PSI 
.06842  LBM/FT3 
.00010 
1.00 
6250 
.909 
1 . 03280 
. 93830 
.02128  LBM/SEC 
18.662  CFM 
6259 


132 


♦♦ROTOMETER  CALIBRATION** 
ASME  HERSCHEL-TYPE  VENTURI 
4.26/2. 13  IN  DIA 


60  PERCENT  READING 


**EXPER I MENTAL  DATA** 


ATMOSPHERIC  PRESSURE  -  29.76910  IN  Hg 
VENTURI  INLET  TEMPERATURE  =  72.0  DEG  F 
VENTURI  INLET  STATIC  GAGE  PRESSURE  =  37.00  IN  H20 
VENTURI  DIFFERENTIAL  PRESSURE  -  .031  IN  RED  OIL 
SPECIFIC  GRAVITY  OF  RED  OIL  -  .834 


♦♦CALCULATED  VALUES** 


ATMOSPHERIC  PRESSURE  = 

VENTURI  INLET  GAGE  PRESSURE  * 
VENTURI  INLET  ABSOLUTE  PRESSURE 
VENTURI  DIFFERENTIAL  PRESSURE  - 
VENTURI  INLET  AIR  DENSITY  - 
VENTURI  PRESSURE  DROP  RATIO  <X> 
EXPANSION  FACTOR  <Y>  - 
ASSUMED  REYNOLDS  NUMBER  - 
DISCHARGE  COEFFICIENT  (C)  - 
VELOCITY  OF  APPROACH  FACTOR  (E) 
FLOW  COEFFICIENT  (K)  - 
AIR  MASS  FLOW  RATE  » 

AIR  VOLUMETRIC  FLOW  RATE  - 
VENTURI  INLET  REYNOLDS  NUMBER  « 


14.6259  PSI 
1.33663  PSI 
13. 28923  PSI 
.00154  PSI 
.06748  LBM/FT3 
.00012 
1.00 
6700 
.911 
1 . 03280 
. 94036 
.02280  LBM/SEC 
20.277  CFM 
6707 


Table  V.  Experimental  Heat  Transfer  Data  (6.25  in.  Bed) 
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Table  VI.  Experimental  Heat  Transfer  Data  (6.25  in.  Masked  Bed) 
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Table  VII.  Experimental  Heat  Transfer  Data  (8. 
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Table  VIII.  Experimental  Heat  Transfer  Data  (10.0  in.  Bed) 
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Table  IX.  Experimental  Heat  Transfer  Data  (12. 


s- 

4) 

O  *-» 
—  03 
JU  O 


ON 

sO 

a 

f-l  nO 

00 

nO 

0 

^4 

00 

as 

CN 

*■4 

0 

A 

O 

00 

ca 

A 

vO 

CN 

O 

00 

A 

•4- 

CN 

CN 

*—• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

ca 

vO 

A 

CN 

a 

A 

A 

'4’ 

A 

CN 

CA 

A 

A 

vC 

>—  I-. 
4-1  O 
O 

26  2 

r-  "-** 

S  a 


ca 

nO 

CN 

A 

A 

*-«  A 

00 

a 

ON 

A 

vO 

ao  no 

as 

CA 

CA 

CN 

CN 

0 

—  0 

0 

NO 

NO 

ao 

ON  O 

CA 

CA 

CA 

CA 

CA 

<n 

CA 

CA  CA 

CA 

CA 

CA 

CA 

CA 

CA 

CA  ^ 

-4 

,<D 

,  H  T3 
.  ca  <a< 
IH  •  ffl 


<v|  ON 

A 

as 

A 

00 

>0  0 

ON 

NO 

NO 

'4* 

as 

•m 

lA 

CN 

CA 

NO 

ON 

NO 

CN 

as 

CN 

sO 

<r 

CN 

O 

CN 

A 

NO 

nO 

O 

O 

ON 

ao 

A 

?NJ  O 

0  0 

<r 

O 

00 

00 

ON 

A* 

00 

on 

Os 

00 

00 

00 

00 

00 

ao 

CO 

ao 

ao 

00 

00 

ao 

00 

00 

00 

<A 

CN 

ao 

ON 

0 

00 

a\ 

A 

ao 

as 

<r 

vO 

A 

O 

CA 

NO 

A 

NO 

0 

CN 

CA 

00 

CN 

ao 

t— < 

vO 

<f  0 

ON 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

A 

A 

ao 

nO 

as 

CN 

CA 

sr 

CA 

\0 

NO 

•41 

A. 

A 

O 

• 

• 

• 

• 

A 

A 

sO 

vO 

sC 

NO 

vO 

vO 

NO 

A 

«■* 

O 

O 

A 

NO 

sO 

A 

*-< 

— 

^•4 

^4 

as 

ON 

ON 

as 

__  O  4-> 
*  r  fl 
,  2  4) 

It-  •  =c 


\OODOC'IOO<>>4'00\0<T'CNr»vOr>'l>-«->— 

n  (s  n  n  n  -  «  »  »  «  rs  «  — •  —  — <oo 

««®aooo®M^rsrsNaoooooeoooa) 


'f,  >stu 

(/)  +J  1 

(/)  r-H 


NO 

A 

nO 

on  ao 

A 

4 

A 

ON 

*4 

CA 

ON 

vO 

ON 

00 

00 

A 

^4 

A 

>y 

A 

© 

A  O 

A 

ON 

A 

0 

00 

CA 

ON 

• 

• 

• 

• 

• 

• 

« 

• 

• 

• 

• 

• 

• 

• 

• 

• 

CN 

co 

CA 

ON 

A 

A 

00 

CA 

00 

A 

ao 

-0 

© 

vO 

nO  O  © 

*-4 

CN 

CA 

CA 

<r 

NT 

A 

A 

NO 

A 

A 

00 

ON 

ON 

*4 

^4 

*■4 

H 

^4 

^4 

3  c 
V*  o  u 

■rtHW 

<  a,  w 


nooi*iooci«MNNr*Nr'Nr>NP'N 
tn  —  a3<T-"r'<ror«<*iOvO<-i<J'so<Na' 

r—*  <—*  —  **  "H  •— «  CN  CN 


0 

A 

4 

nO 

00  0 

A 

sO 

00 

O 

A 

NT 

sC 

00 

© 

A 

CA 

CA 

A 

CA 

CA  nT 

NT 

<r 

A 

A 

A 

A 

A 

nO 

vO 

APPENDIX  D. 

SAMPLE  THERMOCOUPLE  CALIBRATION  DATA 


THERMOCOUPLE  t  14 
i  t  t  i  ************* 


TEMP 

PLAT 

SET 

RESIST 

**** 

****** 

1 

27 . 4555 

27 . 5 5 '3 5 

7 

28 . 1373 

4 

28 . 6638 

5 

23 . 4563 

6 

23 . 3557 

( 

30 . 7738 

31 . 3107 

Q 

w  -  l  y  y » 

10 

32 . 7447 

1  1 

3  3  6  3  6  4 

12 

33 . 3723 

13 

33 . 3723 

1  4 

32 . 4532 

15 

31  .  3132 

16 

31  0431 

17 

30 . 5287 

IS 

23  7070 

13 

23.2117 

20 

28.4275 

TC 

BATH 

OUTPUT 

TEMP 

****** 

**** 

.  733 

1  8 . 53 

.  774 

13.61 

1  .  024 

25  82 

1  .221 

30 . 53 

1 . 543 

38  33 

1  .  753 

43  35 

2.107 

51 . 43 

2 . 340 

56 . 85 

65 . 63 

2. 374 

71  .  20 

3  401 

80 . 65 

3.531 

83  53 

3  257 

77 . 50 

2 . 344 

68 . 33 

2  605 

62  32 

2 . 225 

54.24 

2 . 002 

43 . 05 

1 . 652 
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1  121 
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APPENDIX  E. 

SAMPLE  HEAT  TRANSFER  DATA  REDUCTION  COMPUTER  PRINTOUTS 

Non- fluidized  6 . 2 S  in.  configuration. 

**£XPER I MENTAL  DATA** 

RUN  #  06B20 
06/08/82  -  0425 


RUN  NUMBER  ■ 

06B20 

BED  WIDTH  - 

6.25  INCHES 

STATIC  BED  HEIGHT  « 

12.5  INCHES 

AIR  FLOW  RATE  ■ 

7.29  CFM 

AMBIENT  TEMPERATURE  - 

74.50  DEG  F 

RIGHT  HEATER  VOLTAGE  - 

44.0  VOLTS 

RIGHT  HEATER  CURRENT  » 

.82  AMPS 

LEFT  HEATER  VOLTAGE  ■ 

45.0  VOLTS 

LEFT  HEATER  CURRENT  » 

1.00  AMPS 

BED  EXPANSION  * 

0.0  INCHES 

COMMENTS  t  NON-FLU I D 1 2 ED / NON-EX P ANDES 
BED  WITH  SMALL  LOCAL  SPOUT  AT  10-10 


♦♦EXPERIMENTAL 

DATA** 

RUN  #  06B 

20 

06/08/82  - 

0425 

THERMOCOUPLE  READINGS 

TC 

:#  OUTPUT 

TEMP 

TC# 

OUTPUT 

TEMP 

(mV) 

(DEG  F) 

(mV) 

(DEG  F) 

1 

3.  332 

174.65 

41 

2.  110 

124.98 

4 : 

!  3.331 

174.55 

42 

2.076 

123.54 

!  3 . 336 

174. 76 

43 

1.958 

118.62 

A 

.  .977 

79.23 

44 

2.087 

123.87 

C 

)  ■-> . 

174.55 

45 

2.  112 

124.92 

6 

i  3. 337 

174.81 

46 

2.002 

120.34 

/ 

f  3.343 

175. 05 

47 

.911 

73.51 

£ 

1  o. 

174.71 

48 

.911 

73.49 

s 

*  3.312 

173.78 

49 

.910 

73.46 

1C 

)  3.322 

174.  17 

50 

1 . 430 

96.22 

1 1 

3.319 

174.07 

51 

1.408 

95.  30 

i: 

’  3.316 

173.94 

52 

1.385 

94.33 

l' 

J  3.278 

172.36 

53 

.953 

75.34 

I* 

^  3.282 

172.30 

54 

.948 

75.  11 

1! 

i  3.287 

172.72 

55 

1.202 

86.  16 

l<f 

>  3.286 

172.60 

56 

1.422 

95.85 

1/ 

'  3.230 

170. 36 

57 

1.417 

95.  62 

IE 

3  3.235 

170.55 

68 

2.035 

121.70 

IE 

?  3.238 

170.72 

69 

2.  483 

140. 40 

2C 

)  3 . 232 

170.58 

70 

1.952 

118.20 

23 

l  3. 129 

166.36 

71 

2.  Ill 

124.90 

'“>«■ 

2  3.132 

1 66 .51 

72 

1 . 890 

115.58 

nr 

5  3.127 

166.26 

73 

2.072 

123.26 

\  3.129 

1 66 . 35 

74 

.936 

74.29 

iih 

5  3.130 

166.41 

75 

.940 

74.47 

2d 

i>  3.  131 

166.47 

76 

2.  119 

125.24 

»n- 
^  j 

7  3.128 

166. 34 

77 

2.331 

136. 13 

21 

3  3. 128 

166.31 

78 

1.872 

114.81 

2< 

?  3.110 

165. 39 

79 

2.028 

121 .41 

3< 

5  3.113 

165.70 

-» ! 

l  3.111 

1 65 . 63 

3! 

2  3.108 

1 65 . 43 

•->> 

3  3.077 

164.22 

3J 

H  3.075 

164.14 

3! 

5  3.075 

164. 11 

3i! 

b  3.072 

163.96 

7  3.016 

161.74 

3  3.021 

161.93 

9  3.019 

161.83 

EXPERIMENTAL  DATA** 
RUN#  06B20 
06/08/82  -  0423 


RIGHT-HAND  HEATER 
TEMPERATURE  PROFILE 


COL  1 

COL  2 

COL  3 

COL  4 

AVE 

ROW 

1 

174. 65 

174.55 

174.76 

****** 

174.66 

ROW 

2 

174.55 

174.81 

175.03 

174.71 

174.78 

ROW 

j 

173.78 

174.  17 

174.07 

173.94 

173.99 

ROW 

4 

172.36 

172.50 

172. 72 

172.60 

172.54 

ROW 

5 

170.36 

170. 55 

170. 72 

170.58 

170.55 

NOTES  TEMPERATURES  SHOWN  ARE  IN  DEGREES  FAHRENHEIT 
AND  ARE  ARRANGED  AS  SEEN  FROM  WITHIN  THE  BED 


RIGHT-HAND  HEATER 
TEMP  vs  DEPTH 


TEMP  (DEG  F) 


NOTE:  PLOTTED  TEMPERATURES  ARE  AVERAGES  OF  EACH 
HORIZONTAL  ROW 


*  * CALCULATED  RESULTS** 
RUN  #  06B20 
06/08/82  -  0425 


RIGHT-HAND  HEATER 


AVERAGE  HEATER  TEMPERATURE  = 
AVERAGE  PLATE  SURFACE  TEMPERATURE  « 
AVERAGE  BED  TEMPERATURE  = 
ELECTRICAL  ENERGY  INTO  HEATER  * 
LOSS  FROM  HTR  BACKING  TO  ATM  = 
h (HTR  BACKING- ATM)  *  .52 

LOSS  FROM  HTR  PERIMETER  TO  BED  = 
HEAT  FLUX  THRU  PLATE  = 
h< PLATE-BED)  =  3.39 


173.30  DEG  F 
=  173.28  DEG  F 

89.68  DEG  F 
123.11  BTU/HR 
11.22  BTU/HR 
BTU/HR-FT~2-DEG  F 
13.50  BTU/HR 
98.39  BTU/HR 
BTU/HR-FT'2-DEG  F 


LEFT-HAND  HEATER 


AVERAGE  HEATER  TEMPERATURE  » 
AVERAGE  PLATE  SURFACE  TEMPERATURE  = 
AVERAGE  BED  TEMPERATURE  = 
ELECTRICAL  ENERGY  INTO  HEATER  * 
LOSS  FROM  HTR  BACKING  TO  ATM  « 
h (HTR  BACKING-ATM)  »  .42 

LOSS  FROM  HTR  PERIMETER  TO  BED  - 
HEAT  FLUX  THRU  PLATE  = 
h (PLATE-BED)  =  4.77 


164.87  DEG  F 
164.84  DEG  F 
89.68  DEG  G 
153.54  BTU/HR 
9.20  BTU/HR 
BTU/HR-FT-'2-DEG  F 
19.77  BTU/HR 
124.58  BTU/HR 
BTU/HR-FT''2-DEG  F 


ENERGY  BALANCE 


AVERAGE  AIR  INLET  TEMPERATURE  * 
AIR  OUTLET  TEMPERATURE  * 

AVERAGE  BED  TEMPERATURE  * 
q  < AIR  OUT)  = 
q (LOSS  F/R  WALL)  * 
q (TOTAL  OUT  OF  BED)  = 
q( TOTAL  INTO  BED)  = 

SUPERFICIAL  VELOCITY  « 
SUPERFICIAL  MASS  VELOCITY  = 
PARTICLE  REYNOLDS  NUMBER  = 


73.49  DEG  F 
86. 16  DEG  F 
89.68  DEG  F 
96.97  BTU/HR 
.50  BTU./ HR- 
97.  47  BTU/HR 
256.23  BTU/HR 
.47  FT/SEC 
121.92  LBM/HR-FT  2 
2.  44 
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Fluidized  6.24  in.  configuration. 


♦  ♦EXPERIMENTAL  DATA^ 
RUN  #  06B26 
06/0S/S2  -  0310 


RUN  NUMBER  = 

BED  WIDTH  = 

STATIC  BED  HEIGHT  » 

AIR  FLOW  RATE  » 

AMBIENT  TEMPERATURE  - 
RIGHT  HEATER  VOLTAGE  * 
RIGHT  HEATER  CURRENT  * 
LEFT  HEATER  VOLTAGE  * 
LEFT  HEATER  CURRENT  * 
BED  EXPANSION  - 


06B26 
6.23  INCHES 
12.5  INCHES 
9.23  CFM 
74.30  DEG  F 
60.0  VOLTS 
1.12  AMPS 
60.0  VOLTS 
1.2S  AMPS 
.  1  INCHES 


COMMENTS:  FLUIDIZATION  WITH  2-3  INCH 
DIA  BUBBLES  ERUPTING  AT  LOCATION  3-3 
AND  VICINITY/ACTIVE  BOT  LEFT 


-■  v  "j*  * 


.'.-.•.CAH’, 


♦♦EXPERIMENTAL  DATA# * 
RUN  #  06B26 
06/08/32  -  0510 


THERMOCOUPLE  READINGS 


TC# 

OUTPUT 

TEMP 

TC# 

OUTPUT 

TEMP 

(mV) 

(DEG  F) 

(mV) 

(DEG  F) 

1 

2.  729 

150. 44 

41 

2.  146 

126. 48 

2.728 

150.36 

42 

2.  116 

125.21 

j 

2.724 

150. 13 

43 

1 . 996 

1 20 . 22 

4 

.  961 

78.19 

44 

2.  160 

126.91 

5 

2.855 

1 55 . 53 

45 

2.154 

126. 67 

6 

2.859 

155.67 

46 

2.041 

121 . 97 

7 

2.858 

155. 63 

47 

.  922 

74.00 

8 

2.848 

155.21 

48 

•  r 

73.  98 

9 

2.991 

160. 95 

49 

.922 

73.99 

10 

3 . 000 

161.31 

50 

1.564 

101 . 98 

1 1 

2.  996 

161. 15 

51 

1 . 555 

101.63 

12 

2.988 

1 60 . 35 

52 

1 . 551 

101 . 49 

13 

3.  107 

165. 53  ' 

53 

1.556 

101 . 66 

14 

3.  106 

165. 48 

54 

1.588 

1 03 . 00 

15 

3.  106 

165.49 

55 

1.526 

100. 12 

16 

3.  106 

1 65. 43 

56 

1.551 

1 0 1 . 40 

17 

3.  141 

166.80 

57 

1 . 550 

101.35 

18 

3.  144 

166.92 

68 

1.927 

117.  14 

19 

3.  140 

166.81 

69 

1.958 

118.46 

20 

3.  137 

166.79 

70 

1 . 930 

117.27 

21 

'-nO  0 

JL  •  4.J7 

130.99 

71 

2.368 

135.64 

*->*■? 

2.290 

132.28 

72 

1.916 

116.68 

n*i* 

2.291 

132.28 

73 

2.  039 

121.37 

24 

2.297 

132.51 

74 

1 . 045 

79.  10 

25 

2.340 

1 34 . 30 

75 

1 . 260 

88.51 

26 

2.  336 

134.16 

76 

1 . 839 

113.41 

27 

2 . 335 

134.  10 

77 

1.356 

114.  14 

23 

2.  343 

134.63 

78 

1 . 673 

106. 34 

29 

2.346 

134.56 

79 

1.925 

117.06 

30 

2 . 335 

1 34 . 06 

31 

2 . 339 

134.26 

_■»  ^ 

2.354 

134.79 

._•»  .Ji 

2.  298 

132. 51 

34 

2.285 

131.94 

35 

2.  291 

132. 13 

36 

2.  310 

132.97 

—-7 

a 

129. 73 

33 

2.213 

129.16 

39 

'■"»  nnc 

M  «  Jwaa  J 

129. 44 

40 

2.259 

131.11 
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* EXPERIMENTAL  data** 
RUN#  06B26 
06/08/32  -  0510 


RIGHT-HAND  HEATER 
TEMPERATURE  PROFILE 


COL  1 

COL  2 

COL  3 

COL  4 

AVE 

ROW 

1  !  150.44 

150. 36 

150.  IS 

****** 

150.33 

ROW 

t 

2  !  155.53 

155. 67 

155.63 

155. 21 

155.51 

ROW 

1 

3  !  160.95 

161 . 31 

161. 15 

160.85 

161.07 

ROW 

» 

4  !  165.53 

165.48 

165.49 

165.43 

165.48 

ROW 

1 

5  !  166.80 

166.92 

166. 81 

166.79 

166.83 

NOTE:  TEMPERATURES  SHOWN  ARE  IN  DEGREES  FAHRENHEIT 
AND  ARE  ARRANGED  AS  SEEN  FROM  WITHIN  THE  BED 


RIGHT -KANO  HEATER 
TEMP  vs  DEPTH 


TEMP  (DEG  F ) 


NOTE:  PLOTTED  TEMPERATURES  ARE  AVERAGES  OF  EACH 
HORIZONTAL  ROW 
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**£XPER I MENTAL  DATA** 
RUN#  06B26 
06/08/82  -  0510 


FLUIDIZED  BED 
TEMPERATURE  PROFILE 


*  ^CALCULATED  RESULTS** 
RUN  #  06B26 
06/08/82  -  0510 


RIGHT-HAND  HEATER 


AVERAGE  HEATER  TEMPERATURE  * 
AVERAGE  PLATE  SURFACE  TEMPERATURE  = 
AVERAGE  BED  TEMPERATURE  = 
ELECTRICAL  ENERGY  INTO  HEATER  = 
LOSS  FROM  HTR  BACKING  TO  ATM  * 
h ( HTR  BACKING- ATM)  =  .36 

LOSS  FROM  HTR  PERIMETER  TO  BED  = 
HEAT  FLUX  THRU  PLATE  = 
h< PLATE-BED)  =  9.97 


159.34  DEG  F 

159.79  DEG  F 

101.79  DEG  F 
229.29  BTU/HR 

7.92  BTU/HR 
BTU/HR-FT'2-DEG  F 
20.58  BTU/HR 
200.79  BTU/HR 
BTU/HR-FT'V2-DEG  F 


LEFT-HAND  HEATER 


AVERAGE  HEATER  TEMPERATURE  = 
AVERAGE  PLATE  SURFACE  TEMPERATURE  « 
AVERAGE  BED  TEMPERATURE  * 
ELECTRICAL  ENERGY  INTO  HEATER  = 
LOSS  FROM  HTR  BACKING  TO  ATM  = 
h  <HTR  BACKING-ATM)  »  .07 

LOSS  FROM  HTR  PERIMETER  TO  BED  * 
HEAT  FLUX  THRU  PLATE  * 
h (PLATE-BED)  =  21.26 


132.60  DEG  F 
=  132.54  DEG  F 

101.79  DEG  G 
262.05  BTU/HR 
1 . 54  BTU/HR 
BTU/HR-FT--2-DEG  F 
33.49  BTU/HR 
227.02  BTU/HR 
BTU/HR-FT' 2-DEG  F 


ENERGY  BALANCE 


AVERAGE  AIR  INLET  TEMPERATURE  = 
AIR  OUTLET  TEMPERATURE  = 

AVERAGE  BED  TEMPERATURE  = 
q  < AIR  OUT)  * 
q  (LOSS  F./R  WALL)  * 
q (TOTAL  OUT  OF  BED)  = 
q (TOTAL  INTO  BED)  = 

SUPERFICIAL  VELOCITY  = 
SUPERFICIAL  MASS  VELOCITY  = 
PARTICLE  REYNOLDS  NUMBER  = 


73.99  DEG  F 
100.12  DEG  F 
101.79  DEG  F 
247.04  BTU/HR 
26.35  BTU/HR 
273.39  BTU/HR 
481.88  BTU/HR 
.59  FT /SEC 
150.65  LBM./HR-FT  ' 2 
3 . 0 1 


Non- fluidized  12.0  in.  configuration. 


♦  ♦EXPERIMENTAL  DATA^ 
RUN  #  12B30E 
06/17/02  -  0055 


RUN  NUMBER  ■ 

12B30E 

BED  WIDTH  * 

12.00  INCHES 

STATIC  BED  HEIGHT  = 

12.5  INCHES 

AIR  FLOW  RATE  » 

10.53  CFM 

AMBIENT  TEMPERATURE  - 

78.60  DEG  F 

RIGHT  HEATER  VOLTAGE  - 

46.0  VOLTS 

RIGHT  HEATER  CURRENT  * 

.85  AMPS 

LEFT  HEATER  VOLTAGE  =* 

56.0  VOLTS 

LEFT  HEATER  CURRENT  - 

1.02  AMPS 

BED  EXPANSION 

0.0  INCHES 

COMMENTS s  NON-FLUIDI ZED/NON-EXPANDED 
BED/ NO  BUBBLES  BOTTOM  LEFT 


k)  M  kJ  kO  k)  kJ  k)  kJ  kJ  M  kJ  kJ  k>  kJ  kJ  ki  kJ  kJ  kJ  M  U  <A  Cl  Cl  Cl  Cl  Cl  Cl  Cl  Cl  Cl  Cl  Cl  Cl  Cl  Cl  *-*■  d  Cl  Cl 


OUTPUT 

(mV) 


.  556 
.550 
.  553 
.  140 
.554 
.559 
.565 
.  556 
.541 
.551 
.  545 
.541 
.512 
.515 
.516 
.514 
.468 
.471 
.  468 
.458 
.886 
.886 
•  883 
.886 
.893 
.894 
.891 
.890 
.886 
.888 
.887 
.886 
.869 
.867 
.869 
.866 
.828 
.  3'.*6 
.332 
.831 


**EXPER I MENTAL  DATA** 
RUN  #  12B30E 
06/17/82  -  0055 


THERMOCOUPLE  READINGS 


TEMP 

TC# 

OUTPUT 

TEMP 

(DEG  F) 

(mV) 

(DEG  F) 

183.47 

41 

2.  123 

125.52 

183.20 

42 

2.  145 

126.42 

1 83. 34 

43 

2.041 

122. 10 

89.84 

44 

1.785 

111.20 

183.39 

45 

1.862 

114.45 

1  S3 . 55 

46 

1.795 

111.65 

183.83 

47 

.967 

75.99 

183.41 

48 

.966 

75.92 

182.84 

49 

.966 

75.94 

183.22 

50 

1.438 

96.56 

183.03 

51 

1.429 

96.21 

182.80 

52 

1.403 

95.  1 1 

181.62 

JO 

1.273 

89.41 

181.69 

54 

1.277 

89.56 

181.78 

55 

1.257 

88.55 

181.61 

56 

1.261 

88.  86 

179.84 

57 

1 . 260 

88.80 

179.90 

58 

1.265 

89.01 

179.81 

59 

1.274 

89.40 

179.52 

60 

1.255 

88.87 

1 56 . 60 

61 

1.257 

88.95 

156.63 

62 

1.262 

89.  15 

156.46 

63 

1.239 

88.  13 

156.58 

68 

1.954 

118.29 

156.89 

69 

2.477 

140. 15 

156.94 

70 

2.  173 

127.50 

156.32 

71 

2.527  ' 

142.22 

156.75 

72 

1.974 

119. 13 

156. 59 

73 

2.200 

128.63 

1 56. 65 

74 

1 . 095 

81.29 

156.62 

75 

1.  165 

84.36 

156.51 

76 

1.977 

119.26 

155.85 

77 

2.092 

124. 10 

155.76 

78 

1 . 773 

110.61 

155.81 

79 

1 . 740 

109.20 

155.67 

154.16 

154.47 

154.30 

154.51 

a 


i 

« 

V 

\ 

I 

-1 


♦♦EXPERIMENTAL  DATA** 
RUN#  12B30E 
06/ 17/82  -  0055 


RIGHT-HAND  HEATER 
TEMPERATURE  PROFILE 

COL  1  COL  2  COL  3  COL  4 

ROW  1  !  183.47  183.20  183.34  ****** 

.  183.39  183.55  183.83  183.41 

I 

I 

ROW  3  !  182.84  183.22  183.03  182.80 

I 

I 

ROW  4  !  181.62  181.69  181.78  181.61 

I 

I 

ROW  5  i  179.84  179.90  179.81  179.52 


182.97 

181.68 

179.77 


NOTE:  TEMPERATURES  SHOWN  ARE  IN  DEGREES  FAHRENHEIT 
AND  ARE  ARRANGED  AS  SEEN  FROM  WITHIN  THE  BED 


/-V 

=  \  0 
~  3 
3 

o  7 
— ■  6 

:  5 

CO  4 


RIGHT-HAND  HEATER 
TEMP  vs  DEPTH 


TEMP  (DEG  F) 


NOTE:  PLOTTED  TEMPERATURES  ARE  AVERAGES  OF  EACH 
HORIZONTAL  ROW 


155 


**EXPER I MENTAL  DATA** 
RUN#  12B30E 
06/17/82  -  0055 


FLUIDIZED  BED 
TEMPERATURE  PROFILE 


88.  95 

89.01 

88.36 

89.41 

96.56 

89.  15 

89 . 40 

96.21 

88.  13 

88.87 

88.80 

89.56 

95.  11 

TEMPERA1 

rURES  Sb 

40WN  ARE  ] 

[N  DEGRE 

:es  FAHRENJ 

AND  ARE 

ARRANGE 

•D  AS  SEEI* 

^  FROM  fi 

iBOVE  THE  I 

L  5 


♦  ♦CALCULATED  RESULTS** 
RUN  #  12B30E 
06/17/82  -  0055 


RIGHT-HAND  HEATER 


AVERAGE  HEATER  TEMPERATURE  - 
AVERAGE  PLATE  SURFACE  TEMPERATURE 
AVERAGE  BED  TEMPERATURE  * 
ELECTRICAL  ENERGY  INTO  HEATER  * 
LOSS  FROM  HTR  BACKING  TO  ATM  = 
h  (HTR  BACKING-ATM)  =»  .61 

LOSS  FROM  HTR  PERIMETER  TO  BED  = 
HEAT  FLUX  THRU  PLATE  = 
h (PLATE-BED)  *  3.33 


1 82. 26  DEG  F 
182.23  DEG  F 
90.62  DEG  F 
133.41  BTU/HR 
12.85  BTU/HR 
BTU/HR-FT~2-DEG  F 
14.77  BTU/HR 
105.80  BTU/HR 
BTU/HR-FT"2-DEG  F 


LEFT-HAND  HEATER 


AVERAGE  HEATER  TEMPERATURE  * 
AVERAGE  PLATE  SURFACE  TEMPERATURE 
AVERAGE  BED  TEMPERATURE  = 
ELECTRICAL  ENERGY  INTO  HEATER  = 
LOSS  FROM  HTR  BACKING  TO  ATM  = 
h (HTR  BACKING-ATM)  *  .63 

LOSS  FROM  HTR  PERIMETER  TO  BED  = 
HEAT  FLUX  THRU  PLATE  * 
h (PLATE-BED)  =  7.29 


156.03  DEG  F 
155.98  DEG  F 
90.62  DEG  G 
194.90  BTU/HR 
9.90  BTU/HR 
BTU/HR-FT'2-DEG  F 
19.57  BTU/HR 
165.43  BTU/HR 
BTU/HR-FT  '2-DEG  F 


ENERGY  BALANCE 


AVERAGE  AIR  INLET  TEMPERATURE 
AIR  OUTLET  TEMPERATURE  = 
AVERAGE  BED  TEMPERATURE  = 
q ( AIR  OUT)  * 
q (LOSS  F/R  WALL)  = 
q (TOTAL  OUT  OF  BED)  * 
q (TOTAL  INTO  BED)  ■ 
SUPERFICIAL  VELOCITY  * 
SUPERFICIAL  MASS  VELOCITY  = 
PARTICLE  REYNOLDS  NUMBER  = 


75.95  DEG  F 
88.55  DEG  F 
90.62  DEG  F 
138.71  BTU/HR 
16.49  BTU/HR 
155.20  BTU/HR 
305.56  BTU/HR 
.35  FT/SEC 
91.38  LBM/HR-FT  '  2 
1.83 


Fluidized  12.0  in.  configuration. 


♦♦EXPERIMENTAL  DATA** 

RUN  #  12B48E 
06/17/82  -  0133 

RUN  NUMBER  -  12B48E 

BED  WIDTH  -  12.00  INCHES 

STATIC  BED  HEIGHT  -  12.3  INCHES 

AIR  FLOW  RATE  -  16.37  CFM 

AMBIENT  TEMPERATURE  »  81.30  DEG  F 

RIGHT  HEATER  VOLTAGE  *  46.0  VOLTS 

RIGHT  HEATER  CURRENT  »  .83  AMPS 

LEFT  HEATER  VOLTAGE  »  46.0  VOLTS 

LEFT  HEATER  CURRENT  -  1.02  AMPS 

BED  EXPANSION  *  0.1  INCHES 

COMMENTS:  FLUIDIZATION  WITH  2-3  INCH 
BUBBLES  ERUPTING  AT  LOCATION  4-4 


159 


**EXPER I MENTAL  DATA** 
RUN#  1 2B4SE 
06/17/82  -  0155 


RIGHT-HAND  HEATER 
TEMPERATURE  PROFILE 


COL  1 

COL  2 

COL  3 

COL  4 

AVE 

ROW 

1 

180. 14 

180.21 

180.39 

****** 

180.24 

ROW 

2 

180. 16 

180.41 

1 80 .  63 

180.27 

180. 37 

ROW 

3 

179.05 

179.43 

179.27 

179.04 

179.20 

ROW 

4 

176. 96 

177.  13 

177. 12 

176.84 

177.01 

ROW 

5 

174.67 

174.84 

174.72 

174.35 

174.65 

NOTE:  TEMPERATURES  SHOWN  ARE  IN  DEGREES  FAHRENHEIT 
AND  ARE  ARRANGED  AS  SEEN  FROM  WITHIN  THE  BED 


RIGHT-HAND  HEATER 
TEMP  vs  DEPTH 


NOTE:  PLOTTED  TEMPERATURES  ARE  AVERAGES  OF  EACH 
HORIZONTAL  ROW 
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^EXPERIMENTAL  DATA** 
RUN#  12B48E 
06/17/82  ~  0155 


LEFT-HAND  HEATER 
TEMPERATURE  PROFILE 


COL  1 

COL  2 

COL  3 

COL  4 

Mir 

l-ivt 

ROW 

1 

154.38 

154.49 

154.44 

154.68 

154. 

50 

ROW 

155. SO 

155.81 

155. 73 

155.70 

155. 

76 

ROW 

155. 90 

1 55 . 80 

155.65 

155.58 

155. 

73 

ROW 

4 

154.39 

154.  14 

153. 87 

153.69 

154. 

02 

ROW 

5 

151.65 

151. 55 

151.22 

151.43 

• 

151. 

46 

NOTE:  TEMPERATURES  SHOWN  ARE  IN  DEGREES  FAHRENHEIT 
AND  ARE  ARRANGED  AS  SEEN  FROM  WITHIN  THE  BED 


left-hand  herter 

TEMP  vs  DEPTH 


TEMP  (DEG  F> 


NOTE:  PLOTTED  TEMPERATURES.  ARE  AVERAGES  OF  EACH 
HORIZONTAL  ROW 
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♦♦EXPERIMENTAL  DATA** 
RUN#  12B48E 
06/17/32  -  0155 


FLUIDIZED  BED 
TEMPERATURE  PROFILE 


85.01 

81.79 

81.67 

76.  94 

80.  65 

81.87 

82.01 

82.84 

81.82 

81.90 

81. 79 

91#  *->o 

79.41 

NOTES  TEMPERATURES  SHOWN  ARE  IN  DEGREES  FAHRENHEIT 
AND  ARE  ARRANGED  AS  SEEN  FROM  ABOVE  THE  BED 
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♦♦CALCULATED  RESULTS^ 
RUN  #  12B4SE 
06/17/82  -  0155 


RIGHT-HAND  HEATER 


AVERAGE  HEATER  TEMPERATURE  = 
AVERAGE  PLATE  SURFACE  TEMPERATURE  * 
AVERAGE  BED  TEMPERATURE  = 
ELECTRICAL  ENERGY  INTO  HEATER  = 
LOSS  FROM  HTR  BACKING  TO  ATM  = 
h ( HTR  BACKING-ATM)  =  .51 

LOSS  FROM  HTR  PERIMETER  TO  BED  * 
HEAT  FLUX  THRU  PLATE  = 
h (PLATE-BED)  =  3.17 


178.29  DEG  F 
178.27  DEG  F 
81.46  DEG  F 
133.41  BTU/HR 
10.82  BTU/HR 
BTU/HR-FT-2-DEG  F 
16. 12  BTU/HR 
106.48  BTU/HR 
BTU/HR-FT'V2-DEG  F 


LEFT-HAND  HEATER 


AVERAGE  HEATER  TEMPERATURE  - 
AVERAGE  PLATE  SURFACE  TEMPERATURE 
AVERAGE  BED  TEMPERATURE  = 
ELECTRICAL  ENERGY  INTO  HEATER  ■ 
LOSS  FROM  HTR  BACKING  TO  ATM  » 
h (HTR  BACKING-ATM)  =  .35 

LOSS  FROM  HTR  PERIMETER  TO  BED  * 
HEAT  FLUX  THRU  PLATE  * 
h (PLATE-BED)  «  5.10 


154.29  DEG  F 
154.26  DEG  F 
81.46  DEG  G 
160.10  BTU/HR 
6.44  BTU/HR 
BTU/HR-FT'2-DEG  F 
24.67  BTU/HR 
128.98  BTU/HR 
BTU/HR-FT' 2-DEG  F 


ENERGY  BALANCE 


AVERAGE  AIR  INLET  TEMPERATURE 
AIR  OUTLET  TEMPERATURE  * 
AVERAGE  BED  TEMPERATURE  = 
q (AIR  OUT)  » 
q (LOSS  F/R  WALL)  = 
q (TOTAL  OUT  OF  BED)  = 
q (TOTAL  INTO  BED)  * 
SUPERFICIAL  VELOCITY  = 
SUPERFICIAL  MASS  VELOCITY  = 
PARTICLE  REYNOLDS  NUMBER  * 


75.64  DEG  F 
82.94  DEG  F 
81.46  DEG  F 
126.26  BTU/HR 
0.00  BTU/HR 
126.26  BTU/HR 
.  .  276.25  BTU/HR 
.55  FT/SEC 
143.54  LBM/HR-FT " 2 
2.37 
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5.  A  typical  data  run.  (Bed  thermocouple  probe  *57  reads 
significantly  higher  than  other  bed  thermocouples. 
Reading  exceeds  highest  heater  thermocouple  reading.) 


NOTE:  ALL  VALUES  WERE  RECALCULATED  AFTER  OMITTING 

THERMOCOUPLE  NR.  57  READING.  REVISED  RESULTS 
ARE  PRESENTED  IN  TABLES  AND  CURVES  PREVIOUSLY 
SHOWN. 


♦♦EXPERIMENTAL  DATA*! 
RUN  ft  06B36M 
06/08/82  -  1630 


RUN  NUMBER  -  06B36M 

BED  WIDTH  -  6.25  INCHES 

STATIC  BED  HEIGHT  -  12.5  INCHES 

AIR  FLOW  RATE  -  12.48  CFM 

AMBIENT  TEMPERATURE  -  74.50  DEG  F 

RIGHT  HEATER  VOLTAGE  -  60.0  VOLTS 

RIGHT  HEATER  CURRENT  -  1.12  AMPS 

LEFT  HEATER  VOLTAGE  -  60.0  VOLTS 

LEFT  HEATER  CURRENT  ■  1.28  AMPS 

BED  EXPANSION  -  0.25  INCHES 

COMMENTS*  VIOLENT  FLUIDIZATION  AT  CARRYOVER 
LIMIT 
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**£XPER I MENTAL  DATA** 
RUN  #  06B36M 
06/08/82  -  1630 


THERMOCOUPLE  READINGS 


TC# 

OUTPUT 

TEMP 

TC# 

OUTPUT 

TEMP 

(mV) 

(DEG  F) 

(mV) 

(DEG  F) 

1 

1 . 930 

117.46 

41 

1.773 

110.81 

2 

1.936 

117.68 

42 

1.774 

110.85 

3 

1.951 

118.31 

43 

1.726 

108.83 

4 

.980 

79 . 43 

44 

1.966 

118.82 

5 

2.004 

120.53 

45 

2.  108 

124.75 

6 

2.004 

120.50 

46 

2.007 

120.55 

7 

2.005 

120.57 

47 

.936 

74.62 

8 

2.000 

120.34 

48 

.  935 

74.56 

9 

2.096 

124.34 

49 

.935 

74.57 

10 

2.089 

124.06 

50 

1.631 

1 04 . 85 

11 

2.081 

123.72 

51 

1.592 

103.22 

12 

2.088 

123.96 

52 

1.621 

104.49 

13 

2.  180 

127.75 

53 

1 . 608 

103.89 

14 

2.  185 

127.92 

54 

1.609 

103. 90 

15 

2.  164 

127.07 

55 

1.579 

102.38 

16 

2.  166 

127. 11 

56 

1.629 

104.75 

17 

2.209 

128.88 

57 

2.260 

131. 38 

18 

2.219 

129.27 

68 

1.582 

102.44 

19 

2.214 

129.09 

69 

1.658 

105.70 

20 

2.208 

128.96 

70 

1.942 

117.78 

21 

1.936 

117.53 

71 

2.025 

121.28 

22 

1.948 

118.03 

72 

1 . 56 

101 . 71 

23 

1.947 

117.96 

73 

1.626 

104.32 

24 

1.944 

117.82 

74 

1.313 

90.81 

25 

2.005 

120.38 

75 

1.497 

98.78 

26 

1.987 

119.65 

76 

1.768 

110.39 

27 

1.997 

. 120.04 

77 

1 . 705 

107. 71 

28 

2.017 

120.87 

78 

1.682 

106.72 

29 

2.061 

122.74 

79 

1 . 740 

109.20 

30 

2.053 

122.35 

31 

2.051 

122.30 

*?r> 

2.056 

122.44 

< 

w--> 

2.071 

123.09 

34 

2.067 

122. 88 

35 

2.075 

123.21 

36 

2.079 

1 23. 38 

37 

2.080 

123.42 

33 

2.058 

122.50 

39 

2.055 

122.37 

• 

.40 

2.094 

124. 25 

166 


**EX PER I MENTAL  DATA** 
RUN#  06B36M 
06/08/82  -  1630 


RIGHT-HAND  HEATER 
TEMPERATURE  PROFILE 


COL  1 

COL  2 

COL  3 

COL  4 

AVE 

ROW 

1 

117.46 

117.68 

1 18.31 

****** 

117.82 

ROW 

o 

120. 53 

120.50 

120.57 

120.34 

120.48 

ROW 

3 

124. 34 

124.06 

123.72 

123.96 

124.02 

ROW 

4 

127. 75 

127.92 

127.07 

127.11 

127.46 

ROW 

5 

128.88 

129.27 

129.09 

128.96 

129.05 

NOTE:  TEMPERATURES  SHOWN  ARE  IN  DEGREES  FAHRENHEIT 
AND  ARE  ARRANGED  AS  SEEN  FROM  WITHIN  THE  BED 


R  I  GHT-HAND-  HERTER 
TEMP  vs  DEPTH 


NOTE:  PLOTTED  TEMPERATURES  ARE  AVERAGES  OF  EACH 
HORIZONTAL  ROW 


n 


^EXPERIMENTAL  DATA** 

RUN#  06B36M 

06/08/82  -  1630 

LEFT-HAND  HEATER 
TEMPERATURE  PROFILE 

COL  1  COL  2 

COL  3  COL  4 

AVE 

ROW  1 

117.53  118.03 

117.96  117.32 

117.84 

ROW  2 

120.38  119.65 

120.04  120.87 

120.24 

ROW  4  :  123.09  122.88  123.21  123.38  !  123.14 

I  I 

I  > 

ROW  5  I  123.42  122.50  122.37  124.25  !  123.13 


NOTE:  TEMPERATURES  SHOWN  ARE  IN  DEGREES  FAHRENHEIT 
AND  ARE  ARRANGED  AS  SEEN  FROM  WITHIN  THE  BED 


J 


LEFT-HAND  HEATER 
TEMP  vs  DEPTH 


TEMP  (DEG  F> 


PLOTTED  TEMPERATURES  ARE  AVERAGES  OF  EACH 


*  ^CALCULATED  RESULTS** 
RUN  #  06B36M 
06/08/82  -  1630 


RIGHT-HAND  HEATER 


AVERAGE  HEATER  TEMPERATURE  »  123.77  DEG  F 
AVERAGE  PLATE  SURFACE  TEMPERATURE  =  123.71  DEG  F 
AVERAGE  BED  TEMPERATURE  *  108.07  DEG  F 
ELECTRICAL  ENERGY  INTO  HEATER  =  229.29  BTU/HR 
LOSS  FROM  HTR  BACKING  TO  ATM  =  2.96  BTU/HR 
h (HTR  BACKING-ATM)  «  .19  BTU/HR-FT"2-DEG  F 
LOSS  FROM  HTR  PERIMETER  TO  BED  »  8.13  BTU/HR 
HEAT  FLUX  THRU  PLATE  =  218.20  BTU/HR 
h (PLATE-BED)  *  40.18  BTU/HR-FT^2-DEG  F 


LEFT-HAND  HEATER 


AVERAGE  HEATER  TEMPERATURE  »  121.36  DEG  F 
AVERAGE  PLATE  SURFACE  TEMPERATURE  =  121.29  DEG  F 
AVERAGE  BED  TEMPERATURE  -  108.07  DEG  G 
ELECTRICAL  ENERGY  INTO  HEATER  *  262.05  BTU/HR 
LOSS  FROM  HTR  BACKING  TO  ATM  »  .00  BTU/HR 
h  (HTR  BACKING-ATM)  *  .00  BTU/HR -FT''2-DEG  F 
LOSS  FROM  HTR  PERIMETER  TO  BED  -  4.41  BTU/HR 
HEAT  FLUX  THRU  PLATE  »  257.64  BTU/HR 
h  (PLATE-BED)  =  56.11  BTU/HR-FT''2-DEG  F 


ENERGY  BALANCE 


AVERAGE  AIR  INLET  TEMPERATURE 
AIR  OUTLET  TEMPERATURE  - 
AVERAGE  BED  TEMPERATURE  = 
q (AIR  OUT)  ■ 
q (LOSS  F/R  WALL)  * 
q (TOTAL  OUT  OF  BED)  = 
q (TOTAL  INTO  BED)  = 
SUPERFICIAL  VELOCITY  * 
SUPERFICIAL  MASS  VELOCITY  * 
PARTICLE  REYNOLDS  NUMBER  = 


74.58  DEG  F 
102.38  DEG  F 
108.07  DEG  F 
353.81  BTU/HR 
22.30  BTU/HR 
376.11  BTU/HR 
488.38  BTU/HR 
.30  FT /SEC 
202.78  LBM/HR-FT'2 
4.05 
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APPENDIX  F. 


EQUIPMENT  LISTING 


DEVICE 

MANUFACTURER 

MODEL  # 

Axial  Flow 

Turbo  Compressor 

Spencer  Turbine 

550  SCFM 

Rotometer 

Fisher  §  Porter 

34  SCFM 

30- Inch 

Water  Manometer 

Mer iam 

33KA35 

Computer 

Hewlett  -  Packard 
Company 

HP  85 

Data  Acquisition 
System 

Hewlett  -  Packard 
Company 

HP  3497A 

Dual  Flexible  Disk 
Drive 

Hewlett-Packard 

Company 

HP  82901M 

Plotter 

Hewlett-Packard 

Company 

HP  722SA 

Printer 

Hewlett  -  Packard 
Company 

HP  82905B 

Copper  Constantan 
Thermocouples 

Omega  Engineering 

10  Gauge 

Regulated  DC 

Power  Supply 

60V,  2  A  (Left) 

Lambda 

LK345A 

Regulated  DC 

Power  Supply 

60V,  2  A  (Right) 

Hewlett-Packard 

Company 

6296A 

Computating 

Bridge 

Rosemount 

920A 

Calibration 

Bath 

* 

Rosemount 

913A 

Strip  Heaters 

250  Watts 

Watlo 

250  Watts 

1 7 1 


SERIAL  £ 

A69S13 

2151 A 50182 

2 14  3J1 296  5 

C1821 

1552A02618 

110 

8044 


APPENDIX  G. 


EXPERIMENTAL  UNCERTAINTY  ANALYSIS 

In  order  to  arrive  at  an  estimate  of  the  difference 
between  measured  or  calculated  values  and  the  true  values 
of  various  parameters,  the  following  uncertainty  analysis 
was  accomplished: 

1 .  TEMPERATURE  MEASUREMENTS 

A.  Thermocouple  Readings  -  Since  the  thermocouple  out¬ 
put  voltage  was  measured  automatically  by  the  Data  Acquisition 
System,  human  error  in  reading  the  output  was  avoided;  however, 
the  values  arrived  at  by  the  System  were  still  subject  to 
random  error  as  well  as  inaccuracies  in  the  empirical  corre¬ 
lation  developed  during  the  calibration  of  the  thermocouples. 
Because  of  these  effects,  it  is  believed  that  the  thermo¬ 
couple  readings  are  accurate  to  within  ±°F. 

B.  Ambient  Temperature  Reading  -  This  measurement  was 
accomplished  with  a  mercury-in-glass  thermometer  having  mini¬ 
mum  graduations  of  1°F  and  a  range  of  -20°F  to  +140°F.  It 

is  believed  that  the  thermometer  could  be  read  correctly  to 
within  ±0.5°F.  Allowing  for  the  error  within  the  device, 
the  uncertainty  of  the  thermometer  reading  is  considered 
to  be  ±1°F, 
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2.  PRESSURE  MEASUREMENTS 


A.  Bed  Pressure  Drop  -  The  minimum  manometer  scale  read¬ 
ing  was  0.1  inch  with  a  total  range  of  30  inches.  Readings 
of  0.05  inches  were  possible.  For  static  measurements,  the 
pressure  readings  are  considered  to  be  accurate  within  0.2%; 
however,  during  conditions  under  which  the  water  column  was 
fluctuating,  the  uncertainty  was  increased  to  as  much  as  2% 
of  the  full  scale  reading. 

B.  Venturi  Inlet  Static  Pressure  -  The  minimum  scale 
reading  was  0.1  inch  with  a  total  range  of  60  inches.  Readings 
of  0.05  inches  were  possible.  The  pressure  readings  are  con¬ 
sidered  to  be  accurate  to  within  0.1%  of  full  scale. 

C.  Venturi  Differential  Pressure  -  The  minimum  scale 
reading  was  0.001  inches  over  a  range  of  2  inches.  The 
reading  was  considered  to  be  accurate  to  within  0.1%  of 
full  scale. 

D.  Atmospheric  Pressure  -  Minimum  scale  readings  of 
0.001  in.  Hg  over  a  range  of  approximately  32  inches  were 
observed.  As  a  consequence,  the  uncertainty  of  the  reading 
was  considered  to  be  0.11  of  full  scale. 

3.  DISTANCE/SIZE  MEASUREMENTS 

A.  Bed  Width/Depth  -  Measurements  were  considered 
accurate  to  within  ±0.0625  in.  or  0.5%  of  the  apparatus 


B.  Bed  Probe  Depth  -  Due  to  the  bending  motion  imparted 
to  the  probes  by  the  particle  activity,  the  depth  measure¬ 
ments  were  considered  accurate  to  within  the  0.25  in.  or 
1.5%  of  the  bed  depth. 

C.  Bed  Expansion  -  Due  to  the  rapidly  fluctuating  ex¬ 
panded  bed  surface  level,  the  bed  depth  is  considered  to 
have  an  uncertainty  of  ±0.5  in.  or  5%  of  bed  depth. 

D.  Venturi  Inlet/Throat  Diameter  -  These  measurements 
are  considered  accurate  to  within  ±0.001  in.  or  0.2 %  of 
diameter . 

E.  Particle  Diameter  -  Determined  by  optical  microscope 
examination  conducted  by  Morgan,  [Ref.  2],  D^  =  0.0122  [510 
microns)  ±  1%. 

F.  Plate/Wall  Thickness  -  Measurements  were  considered 
accurate  to  within  0.02  in.  or  51  of  the  wall  thickness. 

4.  VOLTAGE/ CURRENT  MEASUREMENTS 

Minimum  ammeter  scale  graduations  were  0.1  Amp  with  a 
total  range  of  2  Amps.  Because  the  power  supplies  were  not 
calibrated  or  compared  with  a  known  standard  an  uncertainty 
of  ±10%  of  the  full  scale  deflection  is  assigned  to  the 
current  reading.  The  minimum  voltmeter  scale  reading  was 
1  volt  with  a  range  of  60  volts.  Once  again,  due  to  the 
lack  of  calibration,  an  uncertainty  of  5%  is  assigned  of 
full  scale  deflection. 
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5.  THERMAL  CONDUCTIVITIES 

Due  to  the  variation  in  conductivities  found  in  the 
literature,  the  following  uncertainties  are  assigned: 

A.  kc  =  232  Btu/Hr-Ft-°F  ±0.5% 

B.  kmb  =  0.112  Btu/Hr- Ft -°F  ±10% 

C.  ki  =  0.025  Btu/Hr-Ft-°F  ±10% 

D.  k£  =  0.112  But/Hr-Ft-°F  ±10% 

6.  DENSITY  MEASUREMENTS 

Particle  Density  -  Established  by  Morgan  to  be  154.6 
Lbra/Ft^  ±3%. 


APPENDIX  H 


ROTOMETER  CALIBRATION  COMPUTER  PROGRAM  LISTING 


:C>  PRINTER  IS  701.100 
20  S-.634 

ZO  CLEAR  @  DISP  “  ROTOMETER  CALIBRATION- 

40  DISR  "  “ 

50  DISP  "THIS  PROGRAM  CONVERTS  VENTURI  PRESSURE  DATA  INTO  ROTOMETER  PLOW  RATES" 
60  PAUSE 

70  CLEAR  £  DISP  “ENTER  P«t*(IN  Hg>“ 

SO  INPUT  P 

PO  CLEAR  @  DISP  “ENTER  Tin (DEG  F) “ 

100  INPUT  T 

110  CLEAR  6  DISP  "ENTER  Y“ 

ICO  INPUT  Y 

130  ‘  CALCULATE  Patm(PSl> 

140  P1“P*14. 7/29. 92 

150  !  CALCULATE  Tin (DEG  R) 

160  Tl-T*459.69 

170  CLEAR  9  DISP  "PRESS  CONTINUE  TO  START  PROCESSING* 

ISO  PAUSE 

190  CLEAR  9  DISP  "ENTER  FLOWMETER  READING* 

200  INPUT  F 

210  CLEAR  9  DISP  "ENTER  Pig (IN  K20) * 

220  INPUT  12 

230  CLEAR  9  DISP  "ENTER  DELTA  P(1N  R0> ” 

240  INPUT  13 

250  !  CALCULATE  FLOWRATES 
260  CALCULATE  Pig  (PSD 
270  P2«I2*14. 7/33. 91/12 
2S0  !  CALCULATE  Pi  (PSD 
290  P3-P1-P2 

300  !  CALCULATE  DELTA  P(PSD 

310  P4-I3«62.4«. 834/12/144 

320  !  CALCULATE  X 

330  X-P4/P3 

340  CALCULATE  RHOi 

350  D-P3* 144/33. 34 /T1 

360  !  CALCULATE  MASS  FLOW (LBM/SEC) 

370  CLEAR  9  DISP  "ENTER  ASSUMED  REYNOLDS  NUMBER- 
380  INPUT  R 

390  CLEAR  »  DISP  "ENTER  DISCHARGE  COEFFICIENT  ASSOCIATED  WITH  ASSUMED  REYNOLDS  N 
UMBER" 

400  INPUT  C 
410  E- 1.032796 
420  K-C*E 

430  M«K*3. 56327/1 44* YtSQR (2*32. 174*D*P4* 1 44) 

440  !  CALCULATE  VOLUME  FLOW  RATE (FT3/SEC) 

450  Vl-M/D 

460  !  CALCULATE  VOLUME  FLOW  RATE (FT3/MIN) 

470  V2-V1X60 

480  •  CALCULATE  REYNOLDS  NUMBER 

490  R2-V1 / 1 4 . 25309*144*0*4.26/ 12/ . 000000379/32. 174 
500  CLEAR  *  DISP  "ASSUMED  REYNOLDS  NUMBER  -*,R 
510  ! 

520  DISP  “CALCULATED  REYNOLDS  NUMBER  -*,R2 
530  ! 

540  DISP  "DO  YOU  WISH  TO  REPEAT  CALCULATION  USING  A  DIFFERENT  REASSUMED  REYNOLDS 
NUMBER?  (Y/N>- 
550  INPUT  Wl« 

560  IF  W1»-"Y-  THEN  360 
570  PRINT  CHR* ( 12) 


530  PRINT  USING  5°0 

590  IMAGE  2/  .  CSX.  "  *  *R0T0METEF.  CAL  I  BRAT  ION**" 

6«>';  =,FINT  USING  ©10 

=; 0  IMAGE  CSX . " A5ME  HERSCHEL-~YPS  VENTURI" 
c20  PRINT  USING  630 

=30  IMAGE  ZZX. . " A. 26/2. 13  IN  DIA"620  PRINT  USING  630  :  r 

©40  PRINT  USING  650  !  F 

=50  IMAGE  /.32X.DD. X, "PERCENT  READING" 

ooO  PRINT  USING  670 

o 70  IMAGE  ZZX." - " 

680  PRINT  USING  690 

690  IMAGE  2/.31X. " * *£XPER IMENTAL  DATA**- 
700  PRINT  USING  710  ;  P 

710  IMAGE  2/, 15X. "ATMOSPHERIC  PRESSURE  -". 14X , DD. DDDDD. X IN  Hg“ 

720  PRINT  USING  730  s  T 

730  IMAGE  15X" VENTURI  INLET  TEMPERATURE  - ” , 13X , DD. D. X, "DEG  Fu 
740  PRINT  USING  750  »  12 

750  IMAGE  15X, "VENTURI  INLET  STATIC  GAGE  PRESSURE  , 2X. DD. DD, X. " IN  H2Q" 

760  PRINT  USING  770  }  13 

770  IMAGE  15X, "VENTURI  DIFFERENTIAL  PRESSURE  »“ , 3X . D. ODD. X . " IN  RED  OIL" 
780  PRINT  USING  790  ;  S 

790  IMAGE  15X. "SPECIFIC  GRAVITY  OF  RED  OIL  -",16X,D.DDD 
800  PRINT  USING  810 

810  IMAGE  2/  ,  3 1 X , " * (CALCULATED  VALUES**" 

820  PRINT  USING  830  ;  PI 

830  IMAGE  2/. 15X. "ATMOSPHERIC  PRESSURE  . 17X , DD. DDDD. X , "PSI “ 

840  PRINT  USING  850  1  P2 

850  IMAGE  15X. "VENTURI  INLET  GAGE  PRESSURE  10X , D. DDDDD, X. “PSI " 

860  PRINT  USING  870  j  P3 

870  IMAGE  15X. "VENTURI  INLET  ABSOLUTE  PRESSURE  5X, DD. DDDDD, X, "PSI" 

880  PRINT  USING  890  I  P4 

890  IMAGE  13X, "VENTURI  DIFFERENTIAL  PRESSURE  8X, D. DDDDD, X. "PSI " 

900  PRINT  USING  910  ;  D 

910  IMAGE  15X. "VENTURI  INLET  AIR  DENSITY  SX. D. DDDDD, X. "LBM/FT3" 

920  PRINT  USING  930  ?  X 

930  IMAGE  13X. “VENTURI  PRESSURE  DROP  RATIO  (X)  10X.D.  DDDDD 

940  PRINT  USING  950  J  Y 

950  IMAGE  15X. "EXPANSION  FACTOR  <Y>  -".24X.D.DD 
960  PRINT  USING  970  ;  R 

970  IMAGE  15X, "ASSUMED  REYNOLDS  NUMBER  ,21X,DDDD 
980  PRINT  USING  990  l  C 

990  IMAGE  15X, "DISCHARGE  COEFFICIENT  <C>  -".18X.D.DDD 
lOOO  PRINT  USING  1010  j  E 

1010  IMAGE  15X. "VELOCITY  OF  APPROACH  FACTOR  <E>  -".1 OX, D. DDDDD 
1020  PRINT  USING  1030  !  K 

1030  IMAGE  15X. "FLOW  COEFFICIENT  <K>  21 X , D. DDDDD 

1040  PRINT  USING  1030  5  M 

1050  IMAGE  1SX, "AIR  MASS  FLOW  RATE  , 15X. D. DDDDD, X , "LBM/SEC" 

1060  PRINT  USING  1070  J  V2 

1070  IMAGE  1SX, “AIR  VOLUMETRIC  FLOW  RATE  , 14X, DD. DDD. X , "CFM" 

1080  PRINT  USING  1090  |  R2 

1090  IMAGE  15X, "VENTURI  INLET  REYNOLDS  NUMBER  ■"15X.DDDD 

1100  CLEAR  *  DISP  "DO  YOU  WISH  TO  PROCESS  ANOTHER  SET  OF  READINGS?  (Y/N) " 
1110  INPUT  W2* 

1120  IF  W2*-"Y"  THEN  170 
1130  END 


APPENDIX  I. 

THERMOCOUPLE  CALIBRATION  COMPUTER  PROGRAM  LISTING 


Ifi  n  I  C;  P  *'  11 

20  D I  SP  "  D H CAL" 

“0  0 1 SP  “  " 

40  D I  P  “THIS  C R L I BRAT  I  OH  PROGR 
AM  IS  FOR  USE  WITH  THE  HP--^ 
?7  OATH  ACQU I S I T I  OH  SYSTEM" 
50  DISP  “  “ 

SO  DISP  “UP  TO  50  THERMOCOUPLES 
MAY  BE  CALIBRATED  S I MULT  AMI 
OUSLY" 

70  DISP  "  “ 

so  DISP  "UP  TQ  20  DIFFERENT  BAT 
H  TEMPERATURE  SETTINGS  MAY  B 
E  USED" 

SO  DISP  "  * 

100  DISP  "THERMOCOUPLE  NUMBERS  S 
HOULD  AGF ZEE  WITH  CHANNEL  NUN 
BERS" 

110  PAUSE 
120  DISP  “  " 

130  OPTION  BASE  1 

140  MASS  STORAGE  IS  "=0701" 

150  SHORT  R  < 20 , 50 > , M  < 20 , 2  > 

ISO  DIM  T ( 20 >  2> 

170  DISP  "  " 

1 S 0  D ISP  "HO  W  M A  N  Y  THE  P  M 0 C  0 U P L  E S 
ARE  BEING  CALIBRATED*" 

ISO  INPUT  K 
200  DISP  "  " 

210  DISP  "HOW  MANY  BATH  THEMPERA 
TURE  SETTINGS  ARE  BEING  USED 
?  " 

220  INPUT  N 
230  H$=" AF" 

240  DISP  "  " 

250  DISP  "WHAT  IS  FIRST  CHANNEL 
NUMBER?" 

2 SO  DISP  "(TWO  DIGITS)" 

270  INPUT  Y$ 

2 SO  H$C33=Y$ 

2?0  C  =  V A  L  (  Y  $  ) 

300  H*C53="RL" 

31©  DISP  "  " 

320  DISP  "WHAT  IS  LAST  CHANNEL  N 
UMBER?" 

330  DISP  "(TWO  DIGITS)" 

340  INPUT  HtC7J 
350  H T C ? 3  = " V T 1 V R 5 " 

3S0  C*=Hf 
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N 
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i 

% 

3 


V 

\ 

\ 
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370  DISP  “  " 

3S0  DISP  "WHAT  DO  YOU  WISH  TO  CA 
LL  YOUR  DATA  FILE?" 

390  DISP  " <SIX  CHARACTERS  MAX)" 
4O0  INPUT  0# 

410  FOR  1=1  TO  H 
420  DISP  “  “ 

430  DISP  "  TEMPERATURE  SETTING 
#  “  i  I 

440  DISP  “  " 

450  DISP  "WHAT  IS  THE  PLATINUM  B 
RTH  SENSOR  RESISTANCE?" 

460  INPUT  T (  1 ,  1  > 

470  DISP  "  “ 

4S0  DISP  "WHAT  IS  THE  CORRESPOND 
ING  BATH  TEMPERATURE  «X>?“ 
490  INPUT  T<I,2> 

5O0  DISP  "  " 

510  DISP  “TO  START  SAMPLING  PRES 
S  CONTINUE" 

520  PAUSE 

530  CLEAR  709 

540  OUTPUT  70S  ; C$ 

350  FOR  .1=1  TO  K 
560  OUTPUT  709  ."AS" 

=70  ENTER  709  ;  E 
5S0  R < I , J > =BY 1 000 
590  NEXT  J 
600  BEEP 
6ie  » EXT  I 
620  CiSP  "  " 

630  CLEAR  0  CISP  "  SAMPLING 

COMPLETE" 

640  DISP  "  " 

650  DISP  “RESULTS  WILL  BE  PRINTE 
D  MOMENTARILY" 

660  DISP  "  " 

670  DISP  "THERMOCOUPLE  OUTPUT  IS 
IN  M  I  L I VOLTS" 

680  DISP  "  " 

690  DISP  "BATH  TEMP.IS  IN'  CELS  I U 

C*  11 

700  DISP  "  " 

710  PRINT  "DATA  IS  STORED  IN  FIL 
E  NAMED  " 

720  PRINT  USING  730  ;  CHR*(34>,D 
t  .•  CHR-- '  34> 

730  IMAGE  1 1 X jH-6R.fi 
740  CREATE  D$, 1,16500 
750  ASSIGN*  1  TO  Df 
760  E=C- 1 
770  FOR  1=1  TO  K 
780  FOR  J=1  TO  N 
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790  M<J,1)*R<J,I> 

306  M<J,2>*T<J#2> 

S10  PRINT#  1  ;  M(J, l),n<J*2> 

320  NEXT  J 
330  F  =  E+ 1 

340  PRINT  USING  350 
350  IMAGE  2s ,  -tt*t*Xttt**:*XX***:t 
XXXX.XXXXXXXX" 

360  PRINT  USING  370 
870  IMAGE  ttXXXX.XtX.tXX.X.tX.X.XX.XX.XX.t 
*********"  ,/  ,X 
330  PRINT  USING  330  F 
390  IMAGE  s  ,$x,  “THERMOCOUPLE  #■' , 
X,  2D 

300  PRINT  USING  910 

910  IMAGE  €>:,  “  * X  X  X 1 X X  X X  X  X X.  X  X  X  X X  “ 

•  2s ,  X 

320  PRINT  USING  930 
930  IMAGE  "TEMP'S  4X, "PLAT" ,  SX, "T 
C" , 5X, "BATH" 

940  PRINT  USING  350 
950  IMAGE  X>  "SET'S  3X* "RESIST'S  3X 
/ "OUTPUT" , 3X , "TEMP" 

960  PRINT  USING  970 
970  IMAGE  “ **** "  ,  3X,  “ ******:''  .  3>s 
“  *:** XXX  '* ,  3X  ,  » tttt  ”  -  s 
330  FOR  L=1  TO  N 
990  PRINT  USING  1000  >  L,T<L,1>, 
M  SL  ,  1  >  ..  M  <  L  *  2  > 

1  000  I  MAGE  X  ..  DO ,  4X ,  DO  DDDD ,  2X ,  DD 
. DC  j  , 3X • DD . DD 
1010  NEXT  L 

1020  CLEAR  S  DISP  "  THERM 000 

UPLE  #'SF 

1030  .DISP  "READY  FOR  LEAST  SOUAR 
ES  CURVE  FITTING" 

1040  DISP  "  " 

1050  DISP  "ENTER  O  FOR  A  PONER  F 
UNCTION  < Y=AX-'B)  " 

1060  DISP  "  " 

107O  DISP  "ENTER  1  FOR.  AN  EXPONE 
*  NTIAL  FUNCTION  <Y«Re~BX>" 
1030  DISP  -  " 

1090  DISP  "ENTER  2  THRU  10  FOR  A 
POLY  'OMIAL  OF  N  TERMS  < DEG 
PEE  =  iM)" 

1  1  00  D I SP  "  <  Y=C  1  +C2X^C3Xa-2+C4X^3 
+  .  .  .  C \ N > XA N-  1  >)" 

1110  DIM  Y  <  20  >  ,  X  •:  20  >  .  A  <  1 0 ,  1  o  y  ,  E:  < 
10/  10  >  ..  C  <  1  ©  ,  D  •:  1 0  '< 

1120  MAT  Y=2EP<M> 

1130  MAT  X»ZEF<N> 

1140  FOR  J«i  TO  H 
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1130  X<J>=M<‘  J,  1  > 

1160  Y  <  J>=M< J  > 2> 

1170  NEXT  J 
1180  INPUT  Q 

1190  REM  "CALCULATE  LOGARITHMS 
F  X-  AND  Y- VALUES  IF  NECES 
fl  R  Y  M 

1200  IF  Q>  =2  THEN  128© 

1210  FOR  J=1  TO  N 

1220  LET  Y< J>=LOG<M< J, 2> > 

1230  NEXT  J 

1240  IF  Q= 1  THEN  1280 

1250  FOR  J=1  TO  N 

1260  LET  X< J>=LOG<M(J, 1 >> 

1270  NEXT  J 

1280  REM  "CALCULATE  ELEMENTS  OF 
A-MATRIX  AND  O-VECTOR" 

1290  LET  N 1 =Q 

1300  IF  NlJs2  THEN  1320 

1310  LET  N 1=2 

1320  MAT  A=2ER<N 1 iNl) 

1330  MAT  0=2ER<N1> 

1340  FOR  J=1  TO  HI 
1350  FOR  K=1  TO  HI 
1360  IF  J+K>2  THEN  1390 
1370  LET  A  <  J  >  lO  =N 
1380  GOTO  1420 
1398  FOR  L=1  TO  N 
14O0  LET  A< J,K>=R< J, K>+X<L J+fc 
-2> 

1410  NEXT  L 

1420  NEXT  K 

1430  FOR  L= 1  TO  N 

1440  IF  J> 1  THEN  1470 

1450  LET  0 < J >  =0 < J )  +  Y <' L > 

1460  GOTO  1488 

14  70  LET  D  •'  J  >  =0  <  J  >  +Y  <  L  >  *X  <  L  >  (  J- 
1  > 

1480  NEXT  L 
1490  NEXT  J 

15O0  REM  "SOLVE  SIMULTANEOUS  LI 
NEAR  EQUATIONS" 

1510  MAT  Es  =  I N V < A ) 

1520  MAT  C  =  B 1 0 

1530  REM  "PRINT  EQUATION  FOR  CUR 
VE  FIT" 

1540  PRINT  USING  1550 

1550  IMAGE  - 


1560  IF  Q> 1  THEN  167© 

1578  LET  Cl *EXP<C<  1  >  > 

1-580  IF  0s  1  THEN  1630 
1590  PRINT  USING  1600 
1600  IMAGE  5X,  "sgXPOMER  FUNCTION 

*** "  ,  ,  x 
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1610  PRINT  MY  =  " s Cl ; "*X~" ; CC2> 
1620  GOTO  I860 
1630  PRINT  USING  1640 
1640  IMAGE  3X, "***EXPONENTIAL  FU 
NOTION***" ,/,X 

1650  PRINT  "Ys" ; Cl ; “*EXPt “ ; C<2> , 
11 $  X  >  " 

1660  GOTO  I860 
167&  IF  C 1 2 > >  =  0  THEN  1740 
1630  PRINT  USING  1630 
1630  IMAGE  3X> "***PGLYNOMI AL  FUN 
CT I  ON***  '* 

1700  PRINT  USING  1710  ;  Q-l 
1710  IMAGE  10X,  "  < ORDER  =  ”  ,  G ,  "  >  "  , 
/  >  X 

1720  PRINT  " Y  = " :C(1);C(2)/ "X" ; 
1730  GOTO  1730 
1740  PRINT  USING  1750 
1750  IMAGE  3  X  .•  “  *  %  *  POLY  N  0  M I A  L  FUN 
CT I  ON*** “ 

1760  PRINT  USING  1770  j  Q-l 
1770  IMAGE  10X/“< ORDER  =“,D," >"> 

S  j  X 

1730  PRINT  " Ys M  ^  c < 1 > ; "  +  "  ;  C  <  2  >  ; " * 

X  “  ; 

1730  IF  Q=2  THEN  I860 
1300  FOR  Js3  TO  Q 
1310  IF  C<J)>*8  THEN  1848 
1320  PRINT  C<J>  J-l , 

1830  GOTO  1858 

1340  PR  I  NT  "  +  C  <  J >  >  "  *X“-  “  ;  J- 1  ; 

1 850  NEXT  J 

1360  REM  “PRINT  INPUT  VALUES  OF 
X  AND  Y  AND  CALCULATED  VALU 
ES  OFY" 

1370  IF  Q>  =2  THEN  1358 
1380  FOR  J=l  TO  N 
1330  LET  Y < J  >  *EXP < Y <  J } > 

1380  NEXT  J 

1316  IF  Q*1  THEN  1358 

1320  FOR  J=l  TO  N 

1338  LET  X  <  J  >  «EXF' <  X  (  J  >  > 

1348  NEXT  J 
1350  PRINT 

1360  PRINT  .USING  1370 
13  70  I  MAGE  / ,  «;• ,  "  X "  ,  6X  •  "  Y  <  ACT  >  * 
•5X,"Y  <  CALC > " 

1330  PRINT  USING  133© 

1330  I  MAGE  3X ,  "  ***  "  .  5X ..  "  **  *****  " 

,  5X,  " *** ***** " 

2800  LET  S=0 
2O10  FOR  J=1  TO  N 
202©  IF  Q>*2  THEN  2080 
403©  IF  Q=1  THEN  206© 


2650  GOTO  2120 

2O60  LET  Y1=C1*EXP<C<2>*X< J;> 
2070  GOTO  2120 
2080  LET  Y 1  =C'  ( 1  > 

2090  FOR  K=2  TO  Q 

2100  LET  Y 1  =Y  1 +C < K >  *-X <  J  )  <  K-  1  > 

2110  NEXT  K 

2120  LET  S=S+<Y< J>-Y1 )-2 
2130  PRINT  USING  2140  ;  X <  J >  ,  Y C  J 
>  ,  Y1 

2140  I  URGE  X  ,  DO  .  ODD  >  5X ..  DO  .  DO ,  6X , 
on  oonn 
2158  NEXT  J 
2160  PRINT 

2178  PRINT  "SUM  OF  SQUARE  ERRORS 
*“  iS 

21S0  CLEAR  S  DISP  "DO  YOU  DESIRE 
ft  PLOT?" 

2190  DISP  "CENTER  Y  OR  N>° 

2200  INPUT  P* 

2210  IF  P*="N"  THEN  3690 

2220  GGSUB  2250 

2230  PAUSE 

2240  GOTO  3690 

2250  DIM  Yl<4>, ftl<5> ,M1<5> 

2260  Y 1  <  1  >  /  Y 1  <  3  >  *  I NF  0  Y1C2>,Y1< 
4>«-INF  ®  Al<l>,fll<2>,fil<3> 
,  A 1  <  4  ft  1  <  5 >  =0 
2279  FOR  J=1  TO  N 
22S0  IF  Y 1 < 2 >  < X <  J >  THEN  Y1<2>*X< 
J  > 

2290  IF  Y 1 < 1 > > X < J >  THEN  Y 1(1) =  X C 
J> 

2300  I F  Y 1  <  4  >  <Y  <  J  >  THEN  V 1 < 4 )  =  Y < 
J) 

2310.  IF  Y 1  < 3 > > V < J >  THEN  Y1<3>«Y< 
J> 

2320  fil  < 1 >=ftl <1 >+X< J>  tj  ft  1  <  2  >  =  ft  1 
(2>+X<J>*X<J)  8  fll<3>  =  RlC3> 

+  Y  <  J  > 

2330  A1(4)«A1<4>+Y<J)*Y<J)  <2  Al< 

5  >  =ft 1 <  5  >  +X ( J>  t-V  <  J  >. 

2340  NEXT  J 

2350  t11Cl)=ftl(l>/N  8  Ml  C2.)  =  Cftl  <2 
>-fll  <  1  >*ftl  (  1  >/'N>/<N-l  >  8  Ml 
<  3 >  =  ft 1 <  3> /N 

2360  Ml C  4  >  =  <  ft  1 C  4  }  -  ft  1 <3>4fil <3.'.'N) 
'<N-1>  8  Ml C  5  >  =  <  ft  1 C  5 ) - ft  1 C 1 ) 
tft  1  C  3  )  ,'N  >  ✓  ft  1  <  2  >  -ft  1  <  1  ,•  4ft  1  <  1 
>/\V> 

2370  R5=M  1  C 3 > -M 1  (  1  '< *M  1(5)  8  05  =  0 
2380  CL’EftR  8  D5=  1  8  DISP  "AUTO  X 
“SCALING? = Y/N" 
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2350  INPUT  R* 

2400  OH  FHR  GOTO  2386/2410,2426 

2410  xi=Yi<n  e  x2=yi<2>  e  goto 

2470 

2420  CLEAR  6  DISP  "ENTER  MI  HI  HUH 
VALUE  FOR  X— AX i S" 

2430  INPUT  XI 

2440  CLEAR  S  DISP  "ENTER  MAXIMUM 
VALUE  FOR'  X-AXIS"  . 

2450  INPUT  X2 

2460  IF  X 1 >  =X2  THEN  BEEP  &  GOTO 
2420 

2470  CLEAR  0  DISP  "DO  YOU  WISH  V 
ERTICAL  OR  HORIZONTAL  LABEL 
S  ON  THE  X-AXIS?  :  V-'H“ 

24B0  INPUT  R * 

2450  IF  NOT  LEN < Rf >  THEN  BEEP  « 

G 0 TO  2 4 7 0 

2300  S?  =  l  >  IF  u=Cf < RtCl , 13>  =  " V" 
THE N  S 5  =  0  0  G 0 TO  2520 
2510  IF  UP£* <R*C1 , 13)#"H"  THEN  B 
EEP  e  GOTO  2470 

2520  CLEAR  0  DISP  "DESIRED  NUMBE 
R  OF  INTERVALS  ON  X-AXIS?  < 
<■16)" 

2530  INPUT  LI 

2540  IF‘L1<1  OR  LI >16  OR  L1#INT< 
LI)  THEN  BEEP  6  GOTO  2520 
2550  CLEAR  <2  DISP  “DESIRED  NUMBE 
R  OF  SUBDIVISIONS  BETWEEN  L 
ABELED  INTERVALS?" 

2560  INPUT  L3 

2570  S8*6  &  IF  L3=0  THEN  S8  =  l  *5 
L3=  1 

25S0  IF  L3 < 1  OR  L3>L1  THEN  BEEP 
0  GOTO  2550 

2556  CLEAR  6  DISP  "DO  YOU  DESIRE 
AUTOMATIC  SCALING  OF  THE  Y 
-AXIS?  =  Y/N" 

2600  INPUT  R* 

2610  ON  FHR  GOTO  2556 . 2676, 2626 
2620  CLEAR  0  DISP  “ENTER  MINIMUM 
DESIRED  VALUE  FOR  Y-AXIS" 
2630  TNPUT  Y 1 

2640  CLEAR  0  DISP  "ENTER  MAXIMUM 
DESIRED  VALUE  FOR  Y-AXIS'.' 
2650  INPUT  Y2 

2660  IF  Y 1  > Y2  THEN  BEEP  @  GOTO  2 
620  ELSE  268© 

2676  Y 1  =  Y 1  3 >.  6  Y2=Y1<4) 

26S0  CLEAR  6  DISP  "ENTER  DESIRED 
NUMBER* OF  INTERVALS  ON  THE 
Y-AXIS  '.<  =  12)  " 
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2720 

2733 

2740 

2750 


2760 


2520  INPUT  L2 

27O0  IF  L2 < 1  OR  L2> 12  OP  L2#INT< 

l2>  then  beep  e  goto  zssq 

2710  CLEAR  G  OISP  "ENTER  DESIRED 
NUMBER  OF  SUBDIVISIONS  BET 
WEEN  LABELED  INTERVALS'"  . 
2720  INPUT  L4 

2730  S7=0  G  IF  L4-0  THEN  S7=l  G 
L4=  1 

2740  IF  L4 < 1  OR  L4>12  THEN  BEEP 
fi  GOTO  2710 

? 1  =X2-X  1  Z3=INT<200/L1  >  *L 

1  e  Z2= V2-V 1  6  24= I NT ( 1 44/L 

2760  01=21  "23  6  02=22/24  >  i5  =  *2 
♦ <207-23 >*01  G  X0=X 1-48*01 
S  Y0=Y1 -40*02 

2770  GCL-EAR  G  SCALE  X0 , 25, YO , Y2+ 

<  151-24  >*02 

2780  XAXIS  Y 1 , 2 1  /L 1 ,  X 1  .•  X2+D 1 
2750  YAXIS  XI , Z2/L2, Y1 , Y2+D2 
2800  IF  S8  THEN  2960 
2810  W1 =LGT < ABS<  2 1  /L  1 *L3>  > 

2820  W=LGT<ABS<X1+<X1=0>>>+1 
2R30  IF  W>5-<SGN<X1 >*-l >  OR  WIn- 
3  THEN  2920 
2«4Q  IF  S9  THEN  3520 
2850  LOIR  90  :_4  , 

2860  FOR  J=Xl  TO  X2+D1  STEP  ^l^L 

1  *L3 

2870  MOVE  J-4*D1,Y0 
2880  GOSUB  3320  • 

2880  LABEL  V*C1,V3 
2900  NEXT  J 
2910  GOTO  2960 
2920  MOVE  X0,Y@ 

2930  LOIR  0  '  A 

2940  I 1 =X 1  G  GOSUB  3450  G  11=21/ 
Li  G  Zf=V*  G  GOSUB  345© 

2950  LABEL  "XMI  N=  " &ZSC  1 ,  J  03'&"  :  T I 
CS=U&V$C1 ,  103 
2960  IF  S7  THEN  3140 
2970  Wl=LGT<ABSvZ2/L2*L4>> 

2980  M=LGT  <  ABS <  Y 1 +  <  Y 1 =0  > >  >  +  l 
2990  IF  H>3-<SGN<Y1)=-1>  OR  Wls- 
3  THEN  3070 

3000  LOIR  0  . 

3010  FOR  J = Y 1  TO  Y2  +  D2  STEP  Z2/L 
2*L4 

3020  MOVE  X0,J-4*O2 
3030  GOSUB  3320 


X2+01  STEP  Zl/L 


G  1 1=21/ 
3450 

.i0D'&"  :  T I 


OJ  CJ  M  C‘J  CO  OJ  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  OJ  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  ‘  O  to  CO  CO  (0  CO  CO  CO  CO 

^  -f*  £»  4*  4*  4*  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  fO  fO  (0  (0  to  fv)  to  (O  IS)  f  O  *■“*■  »-**-**-*►-**-*►-*•-*©©  '5  O  »3j  Q 

(Ti  <ji  ^  co  io  *-k  ^  co  “O  m  4*  co  f  o  m  o  co  co  “0  *Ti  cn  4^  co  io  >-*  «s)  'o  co  ■■'j  ijj  cn  4^  oi  ro  ^  O  *D  CO  “  I  *."h  Cfl  r- 

(S)ij)i3  ©  'S  O  ©  Q  O  CuMHO  *2)  »S  O  '2*  O  Cl*  G*  O  G*  G>  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  O  o 


LABEL  VfCl,V3 
NEXT  J 
GOTO  3140 
LOIR  30 

I  1 =  Y 1  G  GOSUB  3450 
MOVE  X©+12*01,Y1 
LABEL  "YMIN*“ScVfCl  ,  103 
MOVE  X0+24*D1,Y1 
1 1 =22/L2  G  GOSUB  3450 
LABEL  ,,TICS="StVfCl ,  103 
PENUP  e  LOIR  0 
MOVE  X<D,Y<1> 

FOR  J=1  TO  N  G  MOVE  X<J>-2* 
01,  Y<  .0-4*02  G  LABEL  <2 

NEXT  J  G  BEEP 
FOR  L=0  TO  160 
X 4  =  Y 1 <  D  +L*<  <  Y 1  < 2  >  -Y 1  <  1  >  > / 1 
00  > 

IF  Q>=2  THEN  325G 
IF  Q= 1  THEN  3230 
LET  Y4=C1*X4^C<2> 

GOTO  3230 

LET  Y4=C1*EXP<C<2>*X4> 

GOTO  3290 
LET  Y4=C<1> 

FOP  K=2  TO  Q 

LET  Y  4 = Y  4 + C < K > *  X4-CK-1 > 

next 

PLOT  X4,Y4 
NEXT  L 

BEEP  G  GOTO  3330 
vis* "  e  x=j 

V I = V  A  LI <  X  > 

{F  PO S<V$/"E">  THEN  3400 
G9=LGT<  BBS (X+<  X  =  0  >  > > 

IF  LEN<V*»3  AND  BBS G 3 '•  > 4 - 
<SGN<X>«-1 >  THEN  V=5  6  Vf=“ 

- u  e  RETURN 

IF  LEN < V$ ) <5  THEN  V=LEH< V*> 
G  RETURN 

V*C1 , 53=VfiL$CX>  G  V  =  5 

GRAPH  G  RETURN 

EO=POS<V*, "E-X 

IF  V  f  C  l ,  13="-“  THEN  V*C33=V 

ICE03  G  GOTO  3430 

Vf C23=ViTE03 

V*LEN<V#>  G  IF  V>  5  THEN  PR I 
NT  USING  3560  ;  J  G  V=5  6  V 

w  n 

RETURN 
Vt  =  "  M 

VT=VALT(  ID 


3470  IF  POS<V*« "E">  THEN  3506 
348©  VfCl, lCO  =  VAL*<Il> 

3490  -RETURN 

3500  V*C6,  lCO='.,'tCPO£<V#,  “Eu>3 
3510  RETURN 
3520  LDIR  0  8  L3=-INF 
3530  FOR  J=X1  TO  X2  +  D1  STEP  Zl^L 
1*L3 

3540  GOSUB  3320 

3550  IF  L9> J-<V*4+€>*01  OR  L9>Z5 
+(1-V*3>*D1  THEN  PRINT  USIN 
G  3580  8  GOTO  3610 
3560  IMAGE  “LABEL  DELETED  AT  « , 7 
0 . 40 

3570  MOVE  J+  < 2-V44  >  *D  1 ,  Y 1  - 1  2402 
3580  L3=J+< V44+2) 40 1 
3590  IF  L3>Z5  THEN  MOVE  Z5+C2-V* 
8>*D1 , Y1-124D2 
3600  LABEL  V*C1 , V3 
3610  NEXT  J 
3620  GOTO  2960 
3630  DEF  FNR 

3640  IF  NOT  LEN<Rf >  THEN  11=1  6 
GOTO  3660 

3650  I  1  =P0$ <  "  YN "  ..  IJPCS  (  R$C  1  .  13)5  + 
1 

3660  IF  11*1  THEN  BEEP 
3670  FNR= 1 1 
3680  FN  END 

3630  CLEAR  8  DISP  “DO  YOU  WISH  T 
0  FIT  A  DIFFERENT  TYPE  CURV 
E?< Y/N> " 

3780  INPUT  M t 

3710  IF  W*«"Y"  THEN  3738 

3720  GOTO  3740 

3730  GOTO  1020 

3740 'NEXT  I 

3750  ASSIGN#  1  TO  * 

3760  BEEP 

3770  CLEAR  8  DISP  "  CAL  I  BRAT 

ION  COMPLETE" 

3780  END 
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HEAT  TRANSFER  DATA  ACQUISITION  COMPUTER  PROGRAM  LISTING 

::  mass  storage  is  "sD^oi" 

2  ■'  SHORT  v  (  50  > 

30  DIM  C;3'.7;80>.ft<5>.L(5).Ql*ClC»03 
40  ci*““ 

50  OLE Aft  R  DISF  "  QUICK” 

»C«  DISP  “  ” 

70  DISF'  “THIS  PROGRAM  IS  FOR  USE  IN  SAMPLING  POR  FUTURE  PROCESSING  ” 

30  DISP  "OPERATING  PARAMETERS  AND  TC  OUTPUTS  ARE  STORED  ON  DISK  POR  FUTURE  REDU! 
T I  ON  ” 

«0  PAUSE 

100  CLEAR  3  DISP  "WHAT  IS  THE  DATE?  (MO/DA/YR) " 

110  INPUT  Dl* 

120  CLEAR  8  DISP  “WHAT  IS  THE  TIME?  (MILITARY) “ 

130  INPUT  Tl* 

140  CLEAR  8  DISP  "WHAT  IS  THE  BED  WIDTH?  (INCHES)" 

150  INPUT  X 

160  CLEAR  3  DISP  “WHAT  IS  THE  STATIC  BED  HIEGHT?  (INCHES)” 

170  INPUT  SI* 

180  CLEAR  8  DISP  "WHAT  WAS  THE  HIGHEST  NUMBERED  PROBE  UTILIZED?  (TWO  DIGITS)" 
190  INPUT  L* 

200  L«VAL(L*> 

210  CLEAR  6  DISP  “WHAT  IS  THE  VOLTAGE  TO  TPE  RIGHT-HAND  HEATER?  (VOLTS)" 

220  INPUT  VI 

230  CLEAR  9  DISP  “WHAT  IS  THE  CURRENT  TO  THE  RIGHT-HAND  HEATER?  (AMPS)" 

240  INPUT  II 

250  CLEAR  9  DISP  “WHAT  IS  1>C  VOLTAGE  TO  THE  LEFT-HAND  HEATER?  (VOLTS) “ 

260  INPUT  V2 

270  CLEAR  9  DISP  "WHAT  IS  THE  CURRENT  TO  THE  LEFT-HAND  HEATER?  (AMPS)" 

2B0  INPUT  12 

290  CLEAR  8  DISP  "WHAT  IS  THE  AMBIENT  TEMPERATURE?  (DEG  F) " 

300  INPUT  A1 

310  CLEAR  8  DISP  "WHAT  IS  THE  AIR  FLOW-RATE?  (X> " 

320  INPUT  FI 
330  F2-FH.373 

340  CLEAR  8  DISP  "WHAT  IS  THE  BED  EXPANSION?  (INCHES)" 

350  INPUT  El* 

360  CLEAR  8  DISP  "ENTER  ANY  COMMENTS  YOU  WISH  RECORDED  (MAX  100  CHARACTERS) * 

370  INPUT  01  SC  1 3 

380  CLEAR  8  DISP  “WHAT  DO  YOU  WISH  TO  CALL  YOUR  DATA  FILE?  (SIX  CHARACTERS  MAX) 
390  INPUT  F* 

400  H*-“AF01AL“ 

410  H*C73»"79" 

420  H*C93""VTi VR5TE0* 

430  CLEAR  709 
440  OUTPUT  709  jH* 

450  CREATE  F*, 1.1300 
460  ASSIGN#  1  TO  F* 

470  ASSIGN#  2  TO  "C0EFF2s D700" 

480  PRINT#  1  |  D1«,T1*,X,S1*,L*.L,V1, II, V2, I2,A1,F1,F2,E1*.01* 

490  PRINT  USING  900 

500  IMAGE  "I*****************************" 

510  PRINT  USING  520 

520  IMAGE  3/,3X, "«*EXPER I MENTAL  DATA**" 

530  PRINT  USING  540  J  F* 

340  IMAGE  8X , "RUN  #",X,6A 
550  PRINT  USING  560  |  D1*,T1* 

360  IMAGE  6X, 9A, X, X,4A 
570  PRINT  USING  580 

580  IMAGE  2/.2X, "TC#”,3X, " OUTPUT ",5X, "TEMP" 


PRINT  USING  600 

IMAGE  9X  .  “  imVj  "  .  5>. .  "  :DEG  F> "uISTAOO 
PRINT  USING  c20 

I  MAGE  CX  .  "  — —  "  .  Z X  .  "  — - —  "  .  5)  .  " - " 

C-EAR  >:•  Z  ISP  "PRESS  CONTINUE  WHEN  READY  TQ  BEGIN  SAMPLING" 

PAUSE 

OUTPUT  70®  j " TEC “ 

POP  1-1  TO  7® 

OUTPUT  709  ; "AS" 

ENTER  70®  s  R 

v  a  >  -r*  iooo 

NEXT  I 

OUTPUT  70®  : "TE1TE" 

ENTER  70®  s  S 
BEEP 

clear  <?  disp  "sampling  complete" 

DISP  "RESULTS  WILL  BE  PRINTED  MOMENTARILY" 

READ#  2. 1 
FOR  1-1  TO  7® 

READ#  2  I  C(l> ,C(2> .C(3> .C(4) ,C(5) 

IF  I >L  THEN  810 
GOTO  820 

IF  I <-67  THEN  940 

T  ( I)-C  <  1  >  +C  (2)  tV(  I  >-<-C(3>  tV(  I )  tV  ( I  >  +C  (4>  tV(  I  >  tv  ( I  >  tV(  I )  «C  (5)  tv  ( I  >  tv  ( I  >  tv 


T(n-T(I)  t®/3+32 

PRINT#  1  f  V(I),T(I) 

IF  I #41  THEN  920 
PRINT  USING  870 

IMAGE  3/,2X.  "TC#“.3X. "OUTPUT", 3X, "TEMP" 

PRINT  USING  890 

IMAGE  9X,  " (mV) ",3X, " (DEG  F> " 

PRINT  USING  910 

IMAGE  2X ,  " - " .  3X ,  " - "  ,  SX .  " - " 

PRINT  USING  930  ?  I.V(I),T(I) 

I MAGE  2X . DD , 3 X , DD . DDD , 4 X . DDD . DDDD 
NEXT  I 

ASSIGN#  1  TO  t 
ASSIGN#  2  TO  t 
R(1)-(T(1)+T(2)*T(3) )/3 
R(2>— (T(5> wT (6) *T (7) ♦T !8> >/4 
R(3)-(T(9)+T(10)*T(11) ♦T (12) )  /4 
5  R(4)« (T(13> *T (14) *T (13) *T ( 16) ) /4 
)  R  (5)  «  (T  ( 17) +T  < 18) +T ( 19) +T (20) ) /4 
J  L  ( 1 )  -  <T  (21 )  +T  (22)  *T  (23)  #T  (24)  ) /4 

>  L  (2)  -  <T  (23)  wT  (26)  *J  (27)  «-T  (28)  )  /4 
)  L  (3)  -  <T  (29)  *T  (30)  *T  (31 ) +T (32) ) /4 
)  L  (4)  -  (T  (33)  *T (34) *T (33) +T (36) > /4 
)  L(3)-(T(37)wT(3B)*T(39)wT(40))/4 

>  Gl-(T(41)*T(42)wT(43))/3 
)  62— (T (44) wT (43) *T (46) > /3 
)  P1-(T(47)+T(48>*T(49) )/3 

>  81-0 

>  Xl-0 

>  FOR  1-30  TO  L 

>  IF  1-33  THEN  1130 

>  GOTO  1170 
)  X 1-1 

•  GOTO  1260 

•  IF  1-65  THEN  1190 


;  1 30  GOTO  1210 

; : =0  xi«2 

.2'."'.'  SC'S  1260 

121C  1”  I=oo  THEN  1230 

1221  GOTO  1250 

1230  X 1 =3 

1240  GOTO  1260 

1250  B 1 “B 1 *T ( I } 

1260  NEXT  I 
1270  N1-L-30+1-X1 

1280  HI  -  (R  ( 1 )  -<-R  (2)  i-R  (3)  *R  (4)  fR  <5>  >  /3 

1290  H2«(L(l>-*-L(2)"-L'3>-'-L<4>*L<5>  )  /3 

1300  B1=B1/N1 

1310  El -V 1*11*3. 4121 

1320  E2-V2*I2*3. 4121 

1330  IF  XI >1  THEN  1360 

1340  P2«T<33) 

1350  GOTO  1400 
1360  IF  XI >2  THEN  1390 
1370  P2« (T (55) *T (65) ) /2 
1 380  GOTO  1400 

1390  P2- <T C55) +T <65>*T (66) ) /3 
1400  PRINT  USING  1410  ;  X 

1410  IMAGE  3/,  "BED  WIDTH  -  •• .  DD .  DD .  X ,  “ IN” 

1420  PRINT  USING  1430  j  SI* 

1430  IMAGE  "STATIC  BED  HEIGHT  -  ",5A. -  IN” 

1440  PRINT  USING  1430  s  F2 

1430  IMAGE  "AIR  FLOW  RATE  -  ".DDD.DD."  CFM“ 

1460  PRINT  USING  1470  ;  A1 

1470  IMAGE  "AMBIENT  TEMPERATURE-  “.DD.DD, “  DEG  F1 
1480  PRINT  USING  1490  J  VI 

1490  IMAGE  "RIGHT  HEATER  VOLTAGE-  ",DD.D,"  VOLTS' 
1300  PRINT  USING  1510  ;  II 

1310  IMAGE  "RIGHT  HEATER  CURRENT  -  ",D.DD. "  AMPS' 
1320  PRINT  USING  1330  :  V2 

1330  IMAGE  "LEFT  HEATER  VOLTAGE  -  ”.DD.D. "  VOLTS' 
1340  PRINT  USING  1330  1  12 

1330  IMAGE  "LEFT  HEATER  CURRENT  -  " , D . DD , "  AMPS" 
1560  PRINT  USING  1370  j  El* 

1370  IMAGE  "BED  EXPANSION  -  ”.4A."  INCHES" 

1380  PRINT  USING  1390  ?  01* 

1390  IMAGE  "COMMENTS i  ",!OOA 


6160  END 


APPENDIX  K 


HEAT  TRANSFER  DATA  REDUCTION  COMPUTER  PROGRAM  LISTING 


10  SHORT  v 1 30  > 

00  PRINTER  IS  701.100 
21  PRINT  CHR*  ( 27 >  S< "  LI  1 L " 

30  PRINT  CHR*(27>  VL16D" 

40  DIM  C  (3)  .TOO)  .  R  <3)  ,  L  (3)  , Ol*C 1001 
30  01*-"  " 

60  £9-0 
70  B9-0 

80  CLEAR  $  D ISP  "  PROCESS" 

90  DISP  H  " 

100  DISP  "THIS  PROGRAM  IS  FOR  USE  IN  PROCESSING  DATA  WHICH  HAS  BEEN  PREVIOUSLY  R 
ECORDED. “ 

110  DISP  “  “ 

120  DISP  "HEAT  TRANSFER  COEFFICIENTS  ARE  CALCULATED,  DISPLAYED,  AND  PRINTED." 

130  PAUSE 

140  CLEAR  8  DISP  "WHAT  IS  THE  NAME  OF  THE  DATA  FILE?  <SIX  CHARACTERS  MAX)" 

130  INPUT  F* 

160  CLEAR  8  DISP  “PRESS  CONTINUE  TO  BEGIN  PROCESSING" 

1 70  PAUSE 

100  ASSIGN#  1  TO  F* 

190  READ#  1  j  D1*,T1*.X.S1«.L*,L,V1,I1,V2.I2,A1,F1,F2,E1*,01* 

200  F2-F1*. 32434*. 793133 

210  CLEAR  0  DISP  "DO  YOU  WISH  A  PRINT  OUT  OF  THE  EXPERIMENTAL  DATA?  CY/N) “ 

220  INPUT  R8* 

230  IF  R8*-"Y”  THEN  260 

240  B9-1 

230  GOTO  1110 

260  CLEAR  8  DISP  •  COMMENT  STATEMENT* “ 

270  DISP  "  “ 

280  DISP  01# 

290  DISP  "  * 

300  DISP  "DO  YOU  WISH  TO  MAKE  ANY  CORRECTIONS?  <Y/N) " 

310  INPUT  R9* 

320  IF  R9*-“N"  THEN  360 
330  DISP  "  " 

340  DISP  "ENTER  CORRECTED  COMMENT  STATEMENT" 

330  INPUT  01# 

360  ! 

370  ! 

380  !  C4SPRINT  EXPERIMENT*.  CONDITIONS*** 

390  ! 

400  PRINT  CHR*  < 12) 

410  PRINT  USING  420 

420  IMAGE  3/.30X, "**EX PER I MENTAL  DATA**" 

430  PRINT  USING  440  |  F* 

440  IMAGE  33 X , "RUN  #", X.6A 
430  PRINT  USING  460  |  Dlt.TlS 
460  IMAGE  33X,  GA,  X,  ,  X,  4A 

470  IF  E9«l  THEN  1030 
480  PRINT  USING  490  I  F* 

490  IMAGE  3/.22X, "RUN  NUMBER  «",20X,6A 
500  PRINT  USING  310  |  X 

310  IMAGE  /, 22X, "BED  WIDTH  14X, DD. DD, X, "INCHES" 

520  PRINT  USING  530  j  SI* 

530  IMAGE  /,22X, "STATIC  BED  HEIGHT  -",6X.3A,X. "INCHES’ 

540  PRINT  USING  530  |  F2 

530  IMAGE  / , 22X , "AIR  FLOW  RATE  , I2X , DDD.DD. X, "CFM" 

560  PRINT  USING  570  |  A1 
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AD-A150  785  HEAT  TRANSFER  TO  VERTICAL  FLAT  PLATES  IN  A  RECTANGULAR 
GAS-FLUIDIZED  BED(U)  NAVAL  POSTGRADUATE  SCHOOL  HONTEREV 
CA  D  C  NEILV  JUN  84 


UNCLASSIFIED 


F/G  13/7 


NL 


:?  s»-i  then  1390 
: :3v  wi-o 
i ; 9o  z i *c 

; 200  -Or.  i  » l  TO  40 

::K-  Wl-l  THEN  J-O-i  «•  30T0  1230 

;  200  J  -»  *40 

1000  IP  J'L  THEN  1230 

1240  SOTO  1260 

1250  IP  J.-67  THEN  1040 

12o0  IP  J;79  THEN  Zl-1 

1270  IP  21-1  THEN  1010 

1280  PRINT  USING  1290  I  K. V <K> . T <K> . J. V < J) . T < J) 

1290  IMAGE  13X.DD.3X.DD.DDD,4X.DDD.DD.13X.DD.3X.DD.DDD.4X.DDD.DD 
1 000  GOTO  •  1 030 

1010  PRINT  USING  1320  «  K.V(K>,TOO 

1020  IMAGE  13X.DD,3X.DD.DDD,4X.DDD.DD 

1300  GOTO  1060 

1340  J«J*1  »  Wl-l 

1330  GOTO  1250 

1060  NEXT  K 

1070  ASSIGN#  1  TO  * 

1080  PRINT  CHRs ( 12) 

1390  R  < 1 > • <T  < l ) *T (2) *T (3) > /O 
1400  R<2)-<T<5)*T<6)*TC7>*T<8) > /4 
1410  R<3>-  <T(9)*T  UO)*T<  ll)*T<12)  )  /4 
1420  R(4)-(T<13)*T(14)*T(1S)*T<16) ) /4 
1430  R<3)-<T<17)*T<1S)*T<19)*T<20) ) /4 
1440  L < 1 )  —  (T (21 ) *T <221 *T (23) *T (24) ) /4 
1430  L  (2)  *  <T  <231  *T  (26)  VT  (27)  *T  (28) )  /4 
1460  L (3) - CT (29) *T (30) *T  <31 ) *T  <32> ) /* 

1470  L<4)-<T(33)*T(34)*T(33>*T<36) )/4 
1480  L<3)-(T(37)+T(38)*T<39)*T<40>  >/4 
1490  H1-(R<1)*R(2)*«<3)*R(4)*R<3) ) /3 
1300  H2-  <L  <  1  >  *L  (2)  -M.  (3)  *L  (4)  *L  (3)  ) /5 
1510  IF  B9-1  THEN  3300 
1320  ! 

1330  !  *«*DISPLAV  HEATER  TEMP  SNAPSHOT (R-H> *** 

1340  • 

1330  PRINT  USING  1360 

1360  IMAGE  3/, 30X. " **EXPER I MENTAL  DATA**" 

1370  PRINT  USING  1380  j  F* 

1380  IMAGE  33X, "RUN#",X.6A 
1590  PRINT  USING  1600  |  D1*,T1» 

1600  IMAGE  33 X . BA, X, X.4A 

1610  PRINT  USING  1620 

1620  IMAGE  3/, 32X, "RIGHT-HAND  HEATER" 

1630  PRINT  USING  1640 

1640  IMAGE  31 X, "TEMPERATURE  PROFILE- 

1630  PRINT  USING  1640 

1660  IMAGE  /•23X, "COL  1".3X,"C0L  2",3X,"C0L  3",3X,"C0L  4«,7X."AVE" 

1670  PRINT  USING  1680 

1680  IMAGE  23X,  " -  —  — —  - - » 

1690  PRINT  USING  1700  j  T < 1 )  , T <2> , T (3) ,R <1 > 

1700  IMAGE  17X.-R0W  1 ",  X, " ! " , X.  DDD. DD, 2X, DDD. DD, 2X, DDD. DD. 2X , “****«• ",2X, " ! ",2X. 
ODD. DD 

1710  PRINT  USING  1720 
1720  IMAGE  23X , " ! " , 33X , " ! " 

1730  PRINT  USING  1740  ;  T (3) , T<6) . T (7) . T (B> .R (2) 

1740  IMAGE  1 7X, "ROW  2" , X, " ! ",  X, DDD. DD, 2X , DDD. DD,2X, ODD. DD, 2X , DDD. DD, 2X, " ! " ,2X. DD 
D.DD 
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IMAGE  i  .  CCX  "AMBIENT  TEMPERATURE  •  "  .  SX .  DD .  DD .  X .  “ DEG  F" 

5Sv  PRINT  USING  IPO  :  v: 

1° ’’  IMAGE  ■  ,  CCX ,  "R I  GAT  HEATER  VOLTAGE  ■*' , IX,  DD.  D.  X,  "VOLES'* 

=-r i nt  using  s io  :  : : 

©10  IMAGE  .CCX.  "RIGHT  HEATER  CURRENT  .  ©X. D.  DC. X. " AMPS" 
aCO  PRINT  USING  ©GO  s  VC 

©GO  IMAGE  CCX.  "LEFT  HEATER  VOLTAGE  ■" . oX . DD. D.  X  .  "VOi-TS" 

©40  PRINT  USING  650  :  IC 

oSO  IMAGE  /.CCX. "LEFT  HEATER  CURRENT  . 7X , 0. DD. X . " AMPS" 
o©0  PRINT  USING  ©70  j  El* 

670  IMAGE  /.CCX. "BED  EXPANSION  1 1 X. 4A. X. " INCHES" 

©SO  D9-LEN(01*> 

690  01*lD9*1. 100D-"  “ 

700  IF  D9 :  GO  THEN  740 

710  PRINT  USING  720  j  01SC1.D93 

720  I MAGE  / . CCX . “ COMMENTS « " , X . GOA 

7G0  GOTO  970 

740  K9-C5 

750  IF  014CK9.K93-"  ”  THEN  780 
760  ►.9-K9+1 
770  GOTO  750 

7S0  PRINT  USING  790  j  01*11. K93 

790  I MAGE  / . 22X . " COMMENTS I ” . X , 40A 

BOO  D2-K9+3S 

BIO  IF  D9 >02  THEN  B70 

820  K8-K9+1 

830  IF  K8>D9  THEN  970 

840  PRINT  USING  850  j  01 »CK9*1 . D93 

850  IMAGE  2CX.S0A 

860  GOTO  970 

870  L9-0C 

880  IF  01*CL9,L93«"  “  THEN  910 
890  L9-L9*l 
900  GOTO  880 

910  PRINT  USING  920  I  01*CK9+1.L93 
920  IMAGE  CCX . GOA 
930  LB-L9+1 

940  IF  LB >09  THEN  D9-LB 

950  PRINT  USING  960  >  01«CL9-M ,  D93 

960  IMAGE  CCX. SOA 

970  PRINT  CHR* ( 12) 

980  E9-1 
990  GOTO  360 
1000  ! 

1010  !  ***PRINT  THERMOCOUPLE  READINGS*** 

10C0  ! 

1030  PRINT  USING  1040 

1040  IMAGE  /.30X, "THERMOCOUPLE  READINGS" 

1050  PRINT  USING  1060 

1060  IMAGE  /,13X,"TC*",3X, "OUTPUT", 4X, "TEMP", 14X,"TC#",3X, “OUTPUT", 4X, "TEMP" 
1070  PRINT  USING  1080 

1080  IMAGE  COX, " <aV) ", 4X, " (DEB  F) ", 19X, " <*V> ",4X, " (DEG  F) " 

1090  PRINT  USING  1100 

WOO  IMAGE  13X,  " - ",3X,  " - ",4X," - ",14X," - ",3X." - ".4X," - " 

1110  FOR  1-1  TO  79 
1120  IF  I >L  THEN  1140 
1130  GOTO  1150 
1140  IF  I<-67  THEN  1160 
1150  READ*  1  i  V<I),T<I> 

1160  NEXT  I 
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: *3  :■  print  jsins  ;~oC 
l~sO  IMAGE  ZZX.“  "  .  ZZX  •  “ " 

•  — ~  PRINT  USING  *T3;;  i  T'R>  .T-'iO  .  T  i  ;).T:;  Z R  :  Z  ’ 

:~3C  IMAGE  1  "X .  "ROW  Z'  x. DDD. DD. 2X .  DDD.  0D.2x.0DC.0D.2X.  ODD .  DC .  2* 

D.DD 

17«0  PRINT  USING  1300 
1800  IMAGE  zzx.  •• :  ".—X.  •• :  •• 

1310  PRINT  USING  1820  s  T  ■;  13)  .  T  ( 14)  .  T  { IS)  .  T  ( lo'  .  R  ■  4  > 

1320  IMAGE  1 7X .  "ROW  4" ,  x .  “  :  " .  X.  DDD .  DO.  2X .  ODD.  DO.  2X .  ODD.  00. 2>. .  ODD.  DD.  2X 
D.  00 

13*0  PRINT  USING  1840 
1840  IMAGE  23X ,  "  !  "  , 33X . " ! * 

1830  PRINT  USING  1860  i  T < 17) , T < 1 81 . T ( 19) . T (20> . R (3) 

1860  IMAGE  1 7X . "ROW  3“ . X . " ! ” . X. DDD. DO, 2X. DDD. DD. 2X . DDD- DD, 2X. ODD. DD, 2X 
D.DD 

1870  PRINT  USING  1880 

1880  IMAGE  23X «  " - " 

1890  PRINT  USING  1900 

1900  IMAGE  / , 16X . "NOTE:  TEMPERATURES  SHOWN  ARE  IN  DEGREES  FAHRENHEIT" 
1910  PRINT  USING  1920 

1920  IMAGE  22 X , "AND  ARE  ARRANGED  AS  SEEN  FROM  WITHIN  THE  BED" 

1930  PRINT  USING  1940 

1940  IMAGE  21/. 17 X. "NOTE i  PLOTTED  TEMPERATURES  ARE  AVERAGES  OF  EACH" 

1930  PRINT  USING  1960 

1960  IMAGE  23X. "HORIZONTAL  ROW" 

1970  PRINT  CHR6C12) 

1980  ! 

1990  !  *«*PLOT  TEMP  GRAPH (R-H> «« 

2000  ! 

2010  PRINTER  IS  2 

2020  PRINT  USING  2030  s  F« 

2030  IMAGE  4/.10X.-RUN#  “.X,6A 
2040  PRINT  USING  2030 
2030  IMAGE  3/ 

2060  PRINTER  IS  701,100 
2070  Z9«IP(R(3>) 

2080  ZB«RMD<Z9.3) 

2090  IF  ZB-0  THEN  2120 

2100  Z9-Z9-1 

2110  GOTO  2080 

2120  Z6-Z9-13 

2130  Z7-Z9«M3 

2140  GCLEAR  9  SETGU 

2130  LOCATE  23,118.20.80 

2160  !  FRAME 

2170  SCALE  Z6. Z7, 0, 10 

2180  FXD  0 

2190  LAXES  1,1, Z6. 0,3,1 
2200  Z3-Z7-16 
2210  MOVE  23,-3 
2220  DEG 

2230  LOIR  O  9  LORG  4 
2240  LABEL  "TEMP  (DEG  F) * 

2230  Z4»Z6— 3 

2260  MOVE  Z4, 3 

2270  LDIR  90  9  LORG  4 

2280  LABEL  "VERT  POSITION  (IN)" 

2290  Z3-Z6+14 
2300  MOVE  Z3, 12 
2310  LDIR  0  9  LORG  4 
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232v  L4B£i_  "P.rSHT-HANO  HEATER" 
ztz'.'  ;c»ic-i4 
2340  HOVE  22 . 1 1 

2330  _r  ;s  o  c-  ldpg  4 

23oO  ^ABEl  """EH®  v*  DEPTH" 

2770  3ETUU 
2380  y-9 

2390  POP  1-1  TO  5 
2400  HOVE  R ( I ) . V 
2410  Y-Y-2 
2420  LORG  5 
2430  LABEL 
2440  NEXT  I 
2450  COPY 

2460  I  *  *  *D  I  SPLAY  HEATER  TEMP  SNAPSHOT  (L-H> *»* 

2470  f 

2400  PRINT  USING  2490 

2490  IMAGE  3/.30X. “**EXPERIMENTAL  DATA*** 

2500  PRINT  USING  2510  i  F« 

2510  IMAGE  33X ,  " RUN# " . X , 6A 
2520  PRINT  USING  2530  |  D1*,TI* 

2530  IMAGE  33X. BA. X.4A 

2540  PRINT  USING  25SO 

2550  IMAGE  3/ . 33X, "LEFT-HAND  HEATER- 

2560  PRINT  USING  2570 

2570  IMAGE  31 X. "TEMPERATURE  PROFILE" 

2580  PRINT  USING  2590 

2590  IMAGE  / . 2SX , "COL  1",3X,"C0L  2“,3X.MC0L  3",3X, "COL  4",7X,"AVE" 

2600  PRINT  USING  2610 

2610  IMAGE  23X." - " 

2620  PRINT  USING  2630  >  T (21 ) . T <22> , T (23) , T (24) , L Cl > 

2630  IMAGE  17X.”R0W  1 " , X. “ ! " , X , ODD. DO, 2X. DDD. DO, 2X. DDD.DD.2X. DDD. DO. 2X.  " ! "  .  2X.  DD 
D.DD 

2640  PRINT  USING  2650 
2650  IMAGE  23X, " ! ",33X. " ! “ 

2660  PRINT  USING  2670  j  T <251 , T (25> . T<27> , T (28) . L <2> 

2670  IMAGE  17X. "ROW  2” , X, ” ! " , X. ODD. OD, 2X, DDD. DD, 2X, DDD.DD.2X, DDD.DD.2X. " : " . 2X . DD 
D.DD 

2680  PRINT  USING  2690 
2690  IMAGE  23X. " ! ".33X, " 1 " 

2700  PRINT  USING  2710  j  T <29> , T <30> . T <31 > . T (32) .L (3) 

2710  IMAGE  17X, "ROW  3* , X, " : " , X. DDD. DD, 2X, DDD. DD.2X, DDD.DD. 2X. DDD.DD.2X. " : " . 2X. DD 
D.DD 

2720  PRINT  USING  2730 
2730  IMAGE  23X. " t *.33X. " ! * 

2740  PRINT  USING  2750  |  T <33> , T <34) ,T (35) , T (36) , L <♦) 

2750  IMAGE  17X, "ROW  4" . Xf X, DDD. DD, 2X, DDD. DD, 2X, DDD. DD.2X. DDD.DD, 2X.  "!".  2X.  DD 
D.DD 

2760  PRINT  USING  2770 
2770  IMAGE  23X, " ! ", 33X, " i * 

2780  PRINT  USING  2790  I  T<37),T<38) ,  T  (39) ,  T  <40>  , L  <5> 

2790  IMAGE  17X, "ROW  5" ,  X, "!". X. DDD. DO, 2X, DDD. DD, 2X. DDD.DD, 2X, DDD.DD, 2X, ,2X, DD 
D.DD 

2800  PRINT  USING  2810 

2810  IMAGE  23X," -  -  ■ —  - -  — - " 

2820  PRINT  USING  2830 

2830  IMAGE  /, 16X, "NOTEi  TEMPERATURES  SHOWN  ARE  IN  DEGREES  FAHRENHEIT- 
2840  PRINT  USING  2850 

2850  IMAGE  22 X. "AND  ARE  ARRANGED  AS  SEEN  FROM  WITHIN  THE  StS" 

2860  PRINT  USING  2S70 
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'if  — cviv  I^TVH  *."»  '  w* 


iy  -  ki-  -  rrv  -'*#■ — if  ~w-  lp«-  it-- 


■u-' 1  w»  tvw;w  "-f "  -j""  t*  "'vr  » 


*1 

k' 


I 


rs~0  IMAGE  -I  - ,  17X . "NOTE:  ^LOTTED  TEMPERATURES  ARE  AVERAGES  3”  EACH11 

2330  PRINT  USING  2390 

23°0  IMAGE  23 X. "HORIZONTAL  ROW” 

2«  :>0  PR  I  NT  2HR*  i  1 2 ) 

2«10  ' 

2920  '  »*cL0T  TEMP  GRAPH  »«" 


2®40  PRINTER  IS  2 

2950  PRINT  USING  2960  :  P* 

2®60  IMAGE  4  / ,  1 0  X .  ”  RUN#  “.X.6A 
2®T0  PRINT  USING  2980 
2980  IMAGE  3/ 

2990  PRINTER  IS  701,100 
3000  Z9-IP(L<3>> 

3010  ZB-RMDCZ9.5> 

3020  IF  ZB»0  THEN  3050 

3030  Z9-Z9-1 

3040  GOTO  3010 

3050  Z6-Z9-15 

3060  Z7-Z9*15 

3070  GCLEAR  8  SETGU 

3080  LOCATE  25,118.20.80 

3090  !  FRAME 

3100  SCALE  Z6.Z7.0.10 

3110  FXD  0 

3120  LAXES  1, 1.Z6.0.5. 1 
3130  Z3-Z7-16 
3140  MOVE  ZS.-3 
3150  OEG 

3160  LOIR  0  9  LORG  4 
3170  LABEL  "TEMP  (DEG  F) " 

3180  Z4-Z6-3 

3190  MOVE  Z4.5 

3200  LOIR  90  9  LORG  4 

3210  LABEL  "VERT  POSITION  (IN)" 

3220  Z3-Z6+14 

3230  MOVE  Z3, 12 

3240  LDIR  0  9  LORG  4 

3250  LABEL  "LEFT-HAND  HEATER" 

3260  Z2-Z6-14 
3270  MOVE  Z2. 1 1 
3280  LDIR  0  9  LORG  4 
3290  LABEL  “TEMP  vs  DEPTH" 

3300  SETUU 
3310  Y-9 

3320  FOR  1-1  TO  5 
3330  MOVE  L  ( I  >  .  Y 
3340  Y-Y-2 
3350  LORG  5 
3360  LABEL 
3370  NEXT  I 
3380  COPY 

3390  !  PRINT  PAGE  HEADING 
3400  PRINT  USING  3410 

3410  IMAGE  3/,30X."*«EXPERlHENTAL  DATA*** 

3420  PRINT  USING  3430  !  FS 

3430  IMAGE  35X, "RUN#", X,6A 

3440  PRINT  USING  3450  i  D18.T18 

3450  IMAGE  33X,8A. X, X,4A 

3460  PRINT  USING  3470 
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T*-C-  IMAGE  Z-  .  34X.  "FLUIDIZED  BED" 

~*30  PRINT  USING  3490 

-400  IMAGE  "IX. "TEMPERATURE  PROFILE" 

2500  •  READ  BED  SIZE 

3510  oi*“F*:i.Zj 

3520  IF  Cl*»"12"  THEN  4590 
2530  IF  019--10"  THEN  4220 
2540  IF  Ql»-"08"  THEN  38«0 
2550  IF  P9-1  THEN  2850 
25b0  • 

2570  !  **DISPLAY  6  IN  BED  TEMP  SNAPSHOT** 
2580  ' 


2590 

PRINT 

USING 

3600 

3bOO 

IMAGE 

2/, 27) 

t  '  . 

3b  10 

PRINT 

USING 

3620 

3620 

IMAGE 

27X,  “ 

“ .  25  X ,  “  :  •• 

3430 

PRINT 

USING 

3640 

3640 

IMAGE 

27X,“ 

“ . 25X, “ : " 

3430 

PRINT 

USING 

3660  »  T ( 56  > . T ( 53 ) . T ( 30  > 

3440 

IMAGE 

27X.  “ 

X.DDD.DD,2X.DDD.DD,2X.DDD.DD.2X. " 

3470 

PRINT 

USING 

3680 

3480 

IMAGE 

27X.  " 

" .  25X 

3690 

PRINT 

USING 

3700 

3700 

IMAGE 

27X.  “ 

%23X,":" 

3710 

PRINT 

USING 

3720  »  T(S1> 

3720 

IMAGE 

27X,  " 

”.17X,DDD.DD,2X,“:“ 

3730 

PRINT 

USING 

3740 

3740 

IMAGE 

27X." 

N.25X, “  :  - 

3730 

PRINT 

USING 

3760 

3760 

IMAGE 

27X .  “ 

“ . 25 X . «  i  “ 

3770 

PRINT 

USING 

3780  i  T  (57) .  T  (54)  ,  T  (52) 

3780 

IMAGE 

27X ,  “ 

" , X , DDD. DD. 2X . ODD. DD. 2X , DDD. DD, 2X , " 

3790 

PRINT 

USING 

3800 

3800 

IMAGE 

27X." 

".2Sx, " : " 

3810 

PRINT 

USING 

3830 

3820 

PRINT 

USING 

3840 

3830 

IMAGE 

27X, " ! 

" • 2SX, " i " 

3840 

IMAGE 

27X, 

3850  !  CALCULATE  AVERAGE  BED  TEMP 

38*0  81«  <T  <50)  *T (51 >  *T  <52>  ♦T  (33) *T  (54>+T  (5b)  ♦T  <57 >  )  /7 
3870  GOTO  4900 
3880  ! 

3890  !  t (DISPLAY  8  IN  BED  TEMP  SNAPSHOT** 

3900  IF  B9-1  THEN  4180 
3910  ! 

3920  PRINT  USING  3930 

3930  IMAGE  2/.25X,  ” - " 

3940  PRINT  USING  3950 
3930  IMAGE  25X, * i ",29X, " ! " 

3940  PRINT  USING  3970 
3970  IMAGE  25X, - ! ",29X, « I " 

3980  PRINT  USING  3990  |  T (54) , T (53) , T (50) 

3990  IMAGE  2SX, " ! " , SX . DDD. DD, 2X , DDD. DD,2X, DDD. DD, 2X, " ! ■ 
4000  PRINT  USING  4010 
4010  IMAGE  25X,"i",29X,"S" 

4020  PRINT  USING  4030 
4030  IMAGE  25X,“S".29X,"I" 

4040  PRINT  USING  4050  I  T(51) 

4050  IMAGE  2SX, *  ? ", 21X,DDD.DD,2X, " I  * 

4040  PRINT  USING  4070 
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V*  .  *  .  ,  *  , 


>  -  -p  _* 


y.y;>v 


'•*  •  “  -  *  •**  •**  a**j.**  V***  »**  L*»  **■  .** 

-1.: 


•io-v-  image  25x. 

-030  PR  I  NT  USING  4090 

Image  25>: .  - :  - .  2®>: . . 

AIOv  PRINT  USING  4 1  i  v  ;  T  .  T  {*.4)  .  T -:S2.- 
4  HO  IMAGE  USX,  "  :  " «  SX,  ODD,  DD. 2X .  DDD .  DD . 2x .  ODD.  DD.  2X.  "  :  " 

4 l 20  PRINT  USING  4 I GO 
4 1 30  I  MAGE  25X  .  29X .  "  :  “ 

4140  PRINT  USING  4 1 SO 
4 ISO  IMAGE  2SX . " ! " . 29X ,  “  :  " 

4160  PRINT  USING  4170 

41 BO  !  CALCULATE  AVERAGE  BED  TEMP 

4190  81* (T<50)+T(S1 ) +T (S2) *T  <S3) ♦T  <S4>  +T (36) *T (57) ) /7 
4200  GOTO  4900 
4210  ' 

4220  !  ** (DISPLAY  10  IN  BED  TEMP  SNAPSHOT*** 

4230  IF  B9«l  THEN  4S50 
4240  ! 

42S0  PRINT  USING  4260 

4260  IMAGE  2/.23X." - " 

4270  PRINT  USING  4280 
4280  IMAGE  23X , “ ! “ , 33X , “ ! “ 

4290  PRINT  USING  4300 
4300  IMAGE  23X, " ! “ , 33X. M ! " 

4310  PRINT  USING  4320  t  T (SB) , T (56) . T (S3) . T (30) 

4320  IMAGE  23X. " ! ",  X, DDD. DD. 2X . DDD . DD, 2X , DDD. DD. 2X . DDD. DD. 2X . " ! " 

4330  PRINT  USING  4340 
4340  IMAGE  23X . " ! " . 33X , " ! “ 

4330  PRINT  USING  4360 
4360  IMAGE  23X. " ! ",33X,  " ! " 

4370  PRINT  USING  4380  l  T<39).T<31> 

4380  IMAGE  23X. " i " , X , DDD. DD. 18X, DDD. DD, 2X. " ! " 

4390  PRINT  USING  4400 
4400  IMAGE  23X, " ! " , 33X, " ! " 

4410  PRINT  USING  4420 
4420  IMAGE  23X , " ! " . 33X , " 1 " 

4430  PRINT  USING  4440  J  T <60) , T <S7) , T <34) , T (32) 

4440  IMAGE  23X,  " ! " . X. DDD. DD. 2X, DDD. DD,2X,DDD.DD, 2X, DDD. DD. 2X. " ! " 

4430  PRINT  USING  4460 
4460  IMAGE  23X , " ! " . 33X ,  "  !  “ 

4470  PRINT  USING  4480 
4480  IMAGE  23X, " ! " , 33X, " ! " 

4490  PRINT  USING  4300 

4300  IMAGE  23X.  " - * 

4510  PRINT  USING  4320 

4320  IMAGE  /, 1SX, "NOTE:  TEMPERATURES  SHOWN  ARE  IN  DEGREES  FAHRENHEIT" 

4530  PRINT  USING  4340 

4540  IMAGE  21 X , "AND  ARE  ARRANGED  AS  SEEN  FROM  ABOVE  THE  BED” 

4330  !  CALCULATE  AVERAGE  BED  TEMP 

4360  Bl» <T (30) *T (31 >  *T (32)  *T (33) ♦T  (34)  *T (56) *T (57) +T (38) ♦T (39)  ♦T  (60) ) / 10 

4370  GOTO  4930 
4300  ! 

4390  !  (((DISPLAY  12  IN  BED  TEMP  SNAPSHOT((( 

4600  !  CALCULATE  AVERAGE  BED  TEMP 

4610  Bl»  (T  (30)  *T  (31 )  VT  (32)  *T  (33)  ♦T  (34)  *T  (36)  <*T  (37)  VT  (SB)  *T  (39)  *T  (60)  *T  (61  >  »T  (62) 
♦T(63) )/13 

4620  IF  B9*l  THEN  4930 
4630  ! 

4640  PRINT  USING  4630 

4630  IMAGE  2/.19X,"  - -  - - 
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4Oo0  PRINT  USING  4e70 
4c7C  IMAGE  19X.":“.4IX."  !“ 

4o30  AS.  I  NT  USING  4o90 

A-9o  image  l^x.  “ 

4700  P'lNT  USING  4-10  :  7  16I )  .  T  <58>  .  T  <56)  .  T  (37' .  T  !5C>> 

4-10  IMAGE  19X. X. ODD. DD.2X.DDD. DO. 2X. ODD. DD.2X. DDL. DD.2X. ODD. DD.2X.  •  :• 

4720  PRINT  USING  477>0 

477.0  image  i«x.  :  ».«ix.  * 

4740  PRINT  USING  4750 
4730  IMAGE  19X.  -  I  M.41X.  •• :  ■■ 

4760  PRINT  USING  4770  j  T (62) , T (59) . T (51 > 

4770  IMAGE  19X.”:”.X.DDD.DD.2X.DDD.DD.18X.DDD.DD.2X.*:“ 

4780  PRINT  USING  4790 
4790  IMAGE  19X.  ”  !  “,«1X.  "  ! 

4800  PRINT  USING  4810 
4810  IMAGE  19X. " ! ",41X, * ! “ 

4820  PRINT  USING  4830  «  T (63) , T (60) . T (57) .T (54> , T (52) 

4830  IMAGE  19X. “ I  *.X,DDD.DD,2X.DDD.DD.2X,DDD.DD,2X,DDD.DD,2X,DDD.D0.2X. "1  " 

4840  PRINT  USING  4850 
4850  IMAGE  19X. " ! ",41X, * ! * 

4860  PRINT  USING  4870 
4870  IMAGE  19X, * ! " , 41X, * ! “ 

4880  PRINT  USING  4890 

4890  IMAGE  19X.  “  - - — - -  -  ■  . . ■ . . . .  . . — " 

4900  IF  09-1  THEN  4960 
4910  PRINT  USING  4920 

4920  IMAGE  / , 1SX , "NOTE:  TEMPERATURES  SHOWN  ARE  IN  DEGREES  FAHRENHEIT* 

4930  PRINT  USING  4940 

4940  IMAGE  21 X. "AND  ARE  ARRANGED  AS  SEEN  FROM  ABOVE  THE  BED* 

4950  ! 

4960  !  * • (CALCULATE  RESULTS** 

4970  ! 

4980  !  CALCULATE  AVERAGE  AIR  INLET  TEMP 
4990  PI  —  <T (47) —T <48)  —T (49) ) /3 
5000  !  CALCULATE  AIR  OUTLET  TEMP 
3010  P2-TC33) 

5020  !  CALCULATE  RH  HEATER  GUARD  TEMP 
3030  G1-<T<41)*T<42)4T<43) )/3 
5040  !  PRINT  PAGE  HEADING 
5050  PRINT  CHR«<12) 

5060  PRINT  USING  5070 

5070  IMAGE  3/. SOX, ** (CALCULATED  RESULTS*** 

5080  PRINT  USING  5090  |  F* 

5090  IMAGE  33X.*RUN  #* . X. 6A 
5100  PRINT  USING  5110  1  D19.T1* 

5110  IMAGE  53X.8A, X, *-*,X,4A 

5120  PRINT  USING  5130 

5130  IMAGE  2/,32X, "RIGHT-HAND  HEATER* 

5140  !  PRINT  AVERAGE  HEATER  TEMPERATURE 
5150  PRINT  USING  5160  1  HI 

5160  IMAGE  2/, 15X, "AVERAGE  HEATER  TEMPERATURE  ■* , 12X, DDD.DD, X, "DEG  F* 

5170  •  CALCULATE  AND  DISPLAY  ELECTRICAL  ENERGY  INTO  HEATER 
5180  El— V1«I 1 *3. 4121 

5190  !  CALCULATE  AND  DISPLAY  LOSS  OUT  OF  BACK  OP  HEATER 
5200  S3— <H1— 81) «. 1703163 

5210  •  CALCULATE  CONVECTION  HEAT  TRANSFER  COEFFICIENT  FROM  GUARD  TO  ATMOSPHERE 
5220  Cl -83/ ( .347222* (31 -A 1 ) ) 

5230  !  CALCULATE  LOSS  FROM  SIDES  OF  HEATER  TO  ATMOSHFERE 
5240  S3-C1*. 069444* <T (72) -Al> 

5250  !  CALCULATE  LOSS  FROM  TOP  STRIP  TO  ATMOSPHERE 


-IrO  U2»C1*.  0208333*  iTioSi  -A1 1 

- -  CALCULATE  »C3S  -ROM  BOTTOM  STRI®  TO  ATMOSPHERE 
323':  ST«C  1  *  .  020S3T3*  -  ~  1  — i.  1 ) 

52*0  •  CALCULATE  TOT AL  LOSSES  TC  ATMOSPHERE 
3300  n4«G~'-S”'*-LiO'»S7 

3310  1  CALCULATE  HEAT  FLUX  THRU  RH  PLATE 
5020  P3.E1-H4 
5330  VP«0 

5340  '  CALCULATE  SURFACE  TEMPERATURE  QP  PLATE 
5350  P5»H1-F5*. 0208333/ ( (232-. 032* (Hl-70) > *.347222) 

5360  '  CALCULATE  HEAT  TRANSFER  COEFFICIENT  FOR  PLATE  TO  BED 
5370  P7-P3/ (. 347222* (PS-Bl) > 

S3B0  !  CALCULATE  LOSS  FROM  SIDE  STRIPS  TO  BED 
5390  C3-P7*. 069444* <T (73) -Bl) 

5400  1  CALCULATE  LOSS  FROM  TOP  STRIP  TO  BED 
5410  S9-P7*. 0200333* <T(69)-Bl ) 

5420  !  CALCULATE  LOSS  FROM  BOTTOM  STRIP  TO  BED 
5430  S5-P7* . 020S333*  <T (70) -Bl ) 

5440  !  CALCULATE  TOTAL  LOSSES  INTO  BED 
5450  H3»C3*S9*S5 

5460  !  CALCULATE  NEW  HEAT  FLUX  THRU  PLATE 
5470  X9a»£l  -H4— H3 
5400  V9»V9*1 

5490  !  COMPARE  NEW  PLATE  HEAT  FLUX  WITH  OLD  AND  REPEAT  CALCULATION  IF  DIFFERENCE 
EXCEEDS  0.01 

5500  X8-ABS (P3-X9) 

5510  IF  X8>*. 01  THEN  P3-X9  »  GOTO  5340 

5520  •  PRINT  RESULTS 

5530  PRINT  USING  5540  1  P3 

5540  IMAGE  13X. "AVERAGE  PLATE  SURFACE  TEMPERATURE  •“ , 5X. DDD. DD, X, "DEG  F" 

5550  PRINT  USING  5560  s  Bl 

5560  IMAGE  13X. "AVERAGE  BED  TEMPERATURE  ■“ , 13X.DDD. DD, X. "DEG  F" 

5570  PRINT  USING  5500  I  El 

5500  IMAGE  15X. “ELECTRICAL  ENERGY  INTO  HEATER  »“ . SX , DDD. DD, X , "BTU/HR" 

5590  PRINT  USING  3600  1  H4 

5600  IMAGE  13X. "LOSS  FROM  HTR  BACKING  TO  ATM  9X. DDD. DD. X, “BTU/HR" 

5610  PRINT  USING  5620  »  Cl 

3620  IMAGE  13X. "h  <HTR  BACKING-ATM)  , SX. DDD. DD, X, "BTU/HR— FT'“2— DEG  F« 

5650  PRINT  USING  3640  J  H3 

5640  IMAGE  13X, "LOSS  FROM  HTR  PERIMETER  TO  BED  7X . DDD. DD, X . "BTU/HR" 

5650  PRINT  USING  3660  j  P3 

3660  IMAGE  15X, "HEAT  FLUX  THRU  PLATE  17X, DDD. DD. X, "BTU/HR" 

3670  PRINT  USING  3680  )  P7 

5680  IMAGE  13X,"h (PLATE-BED)  14X, DDD. DD, X, “BTU/HR-FT-2-DEG  F" 

5690  ! 

5700  !  CALCULATE  AND  DISPLAY  RESULTS  FOR  LEFT-HAND  HEATER 

3710  * 

3720  PRINT  USING  3730 

5730  IMAGE  2/, 33X, "LEFT-HAND  HEATER" 

5740  f  CALCULATE  LH  HEATER  GUARD  TEMPERATURE  AVE 
5730  G2-<T<44>VT(43>*T(46)>/3 

5760  !  CALCULATE  ELECTRICAL  ENERGY  INTO  HEATER 
5770  E2-V2* 12*3. 4121 

5780  •  CALCULATE  LOSS  FROM  BACK  OF  HEATER  TO  ATMOSPHERE 
5790  94«<H2-G2>  *. 1703163 

5800  !  CALCULATE  CONVECTION  HEAT  TRANSFER  COEFFICIENT  FROM  GUARD  TO  ATMOSHPERE 
3810  C2-G4/ (.347222* CG2— Al) ) 

3820  ’  CALCULATE  LOSS  FROM  SIDES  QF  HEATER  TO  ATMOSPHERE 
3830  S4-C2*. 069444* (T(7B>—Al> 

5840  !  CALCULATE  LOSS  FROM  TOP  STRIP  TO  ATMOSPHERE 
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5250  U3“C2»  .  02033”*  IT  r~s)  -A1  > 

33oV  ’  CALCULATE  ,_0S5  FROM  BOTTOM  STRIP  TO  ATMOSPHERE 
3S7v  3B-U3 

3330  ‘  CALCULATE  TOTAL  LOSSES  TO  ATMOSPHERE 
3S°C>  ne“&4-S4-u3-58 

3 POO  1  CALCULATE  HEAT  FLUX  THRU  LH  PLATE 
3® 10  P4-E0-HO 

3920  •  CALCULATE  SURFACE  TEMPERATURE  OF  PLATE 
3930  P6-H2-P4* . 0203333/ ( (232-. 022*  <H2-70> >  * . 347222 • 

3®40  1  CALCULATE  HEAT  TRANSFER  COEFFICIENT  FOR  PLATE  TO  BED 
5930  P3-P4/ ( . 347222* (P6-B1 ) > 

3960  !  CALCULATE  LOSS  FROM  SIDE  STRIPS  TO  BED 
3970  C4«PB*. 069444* <T (79) -Bl) 

3900  !  CALCULATE  LOSS  FROM  TOP  STRIP  TO  BED 
3990  U1«P8*. 0208333* <T (77) -B1) 

6000  !  CALCULATE  LOSS  FROM  BOTTOM  STRIP  TO  BED 
6010  S6-U1 

6020  !  CALCULATE  TOTAL  LOSSES  TO  BED 
6030  H3«C4-HJ1*S6 

6040  !  CALCULATE  NEW  HEAT  FLUX  THRU  PLATE 
6030  X7-E2-H6-H3 

6060  '  COMPARE  NEW  HEAT  FLUX  TO  OLD  AND  REPEAT  CALCULATION  IF  DIFFERENCE  EXCEEDS 

0.01 

6070  X6-ABS<P4-X7> 

6080  IF  X6>-.01  THEN  P4«X7  ®  SOTO  5920 

6090  !  PRINT  RESULTS 

6100  PRINT  USING  6110  :  H2 

6110  IMAGE  2/. 1SX, "AVERAGE  HEATER  TEMPERATURE  12X, DDD.DD, X, "DEG  F" 

6120  PRINT  USING  6130  J  P6 

6130  IMAGE  1SX, "AVERAGE  PLATE  SURFACE  TEMPERATURE  -* , 5X, DDD.DD, X. "DEG  F” 

6140  PRINT  USING  6130  :  B1 

6130  IMAGE  13X, "AVERAGE  BED  TEMPERATURE  1SX, DDD.DD, X. "DEG  G" 

6160  PRINT  USING  6170  ;  E2 

6170  IMAGE  13X, "ELECTRICAL  ENERGY  INTO  HEATER  •" , BX. ODD. DD. X. "BTU/HR" 

6180  PRINT  USING  6190  t  H6 

6190  IMAGE  13X. "LOSS  FROM  HTR  BACKING  TO  ATM  9X. DDD. DD, X, "BTU/HR" 

6200  PRINT  USING  6210  j  C2 

6210  IMAGE  13X,  "h  <HTR  BACKING-ATM)  BX, DDD.DD, X, "BTU/HR-FT'~2-DE6  F" 

6220  PRINT  USING  6230  1  H3 

6230  IMAGE  15X, "LOSS  FROM  HTR  PERIMETER  TO  BED  7X, DDD. DD, X, "BTU/HR" 

6240  PRINT  USING  6230  J  P4 

6230  IMAGE  1SX. "HEAT  FLUX  THRU  PLATE  17X, DDD.DD, X, “BTU/HR7 
6260  PRINT  USING  6270  1  PB 

6270  IMAGE  13X, “h (PLATE -BED)  , 14X, DDD. DD, X, "BTU/HR-FT~2-DEG  F" 

6280  ! 

6290  •  CALCULATE  ENERGY  BALANCE  VALUES 
6300  ! 

6310  !  CALCULATE  AIR  DENSITY 
6320  X3- 1.323 *29. 92/  (P2*-439.  69) 

6330  •  CALCULATE  ENERGY  CARRIED  OUT  OF  BED  BY  AIR  STREAM 
6340  A3-F2*X3*60*.241* (P2-P1 > 

6330  !  CALCULATE  LOSS  THRU  FRONT  AND  BACK  FACES 
6360  L6-.448»X«<T(73)-T<74>> 

6370  !  CALCULATE  TOTAL  ENERGY  LEAVING  BED 
6380  A6-A3+L6 

6390  !  CALCULATE  TOTAL  ENERGY  INTO  BED  FROM  RH  HEATER 
6400  B3"€ 1 — H4 

6410  !  CALCULATE  TOTAL  ENERGY  INTO  BED  FROM  LH  HEATER 
6420  B4*ME2“H6 

6430  !  CALCULATE  TOTAL  ENERGY  INTO  BED  FROM  ALL  SOURCES 
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*440  A~»*3-*-&4 

*430  ca_cjlate  convection  heat  transfer  coefficient  f^b  faces-atmosfhefe 
j-iov  ;.e»i_o.  :c**x»  i'T  :?«• -a;  ,  ) 

=4-C'  '  CALCULATE  CONVECTION  HEAT  TRANSFER  COEFFICIENT  F /g  FACES-BED 

»480  5’«Ls-  •  C4»x*  .si-T  (T3)  ,  . 

oAFO  '  CALC'JuATE  SUPERFICIAL  velocity 

qSOO  U*FC/  Co* X /  1  44 )  /oO 

6310  '  CALCULATE  SUPERFICIAL  MASS  VELOCITY 
6500  G-U*X3*3600 

o330  •  CALCULATE  PARTICLE  REYNOLDS  NUMBER 

o340  RH3/3600*. 0009166667/. 000000396/32.  174 

6350  '  PRINT  RESULTS 

6560  PRINT  USING  6370 

6370  IMAGE  C/. 34 X. "ENERGY  BALANCE" 

6580  PRINT  USING  6390  ;  PI 

6390  IMAGE  2/, 13X. "AVERAGE  AIR  INLET  TEMPERATURE  . 9X , DDD. DD. X . “DEG  F" 
6600  PRINT  USING  6610  s  P2 

6610  IMAGE  13X , "AIR  OUTLET  TEMPERATURE  . 16X. DDD. DD. X, "DEG  F" 

6620  PRINT  USING  6630  J  B1 

6630  IMAGE  13X, "AVERAGE  BED  TEMPERATURE  -" . 13X . DDD. DD. X, "DEG  F" 

6640  PRINT  USING  6630  i  AS 

6630  IMAGE  15X, "q (AIR  OUT)  •  " , 27X.DDD. DD. X, "BTU/HR" 

6660  PRINT  USING  6670  I  L6 

6670  IMAGE  13X. "q (LOSS  F/R  WALL)  21 X, DDD. DD, X, "BTU/HR" 

6680  PRINT  USING  6690  )  A6 

6690  IMAGE  13X."q (TOTAL  OUT  OF  BED)  1 BX. DDD. DD. X, "BTU/HR" 

6700  PRINT  USING  6710  l  A7 

6710  IMAGE  15X, "q (TOTAL  INTO  BED)  20X , DDD. DD. X , "BTU/HR" 

6720  PRINT  USING  6730  s  U 

6730  IMAGE  13X, "SUPERFICIAL  VELOCITY  , 17X, DDD. DD. X, "FT/SEC" 

6740  PRINT  USING  6730  S  G 

6730  IMAGE  15X, "SUPERFICIAL  MASS  VELOCITY  •" . 7X, DDD. DD, X, "LBM/HR-FT"2" 

6760  PRINT  USING  6770  J  R 

6770  IMAGE  13X, "PARTICLE  REYNOLDS  NUMBER  -" . 20X, DDD. DD 
6780  END 
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